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THE  BUCKLING  THEORY  OF  TURBULENT  FLOW 


Tne  research  overviewed  in  this  final  report  had  as  primary  objective  the 
establishment  of  a  purely  theoretical  foundation  for  the  existence  of 
turbulent  flows  and  engineering  formulas  that  summarize  the  transport 
|  potential  of  such  flows.  From  the  early  research  -  proposal  stage,  this  work 

reflected  my  view  that  "tradition"  has  given  us  a  fluid  mechanics  research 
methodology  that  is  generally  divorced  from  thermodynamics,  and,  with  special 
i  regard  to  turbulence  research,  a  research  methodology  that  is  divorced  from 

theory.  That  even  the  analytical  work  in  20th  century  turbulence  research  is 
rooted  in  empiricism  is  best  illustrated  by  how  the  empirical  notion  of  "eddy" 
I  is  never  questioned  and  conveniently  hidden  behind  the  time-averaged 

terminology  inherited  from  Osborne  Reynolds. 

The  present  research  project  was  conducted  over  the  three  year  period 
[  that  coincided  with  the  academic  years  1980-81,  81-82,  and  82-83.  A  one-year 

extension  was  added  to  the  original  contract  period  in  order  to  facilitate  the 
completion  of  ongoing  thesis  work  and  my  own  writings  and  preparations  for 
the  1984  sabbatical  leave  and  move  to  Duke  University.  Although  the  actual 
work  is  already  documented  in  detail  in  the  three  annual  reports  [1-3],  I 
welcome  this  opportunity  to  review  it  one  more  time  and  to  stress  its  position 
relative  to  pre-1980  thinking. 

The  present  work  -  the  buckling  theory  of  turbulent  flow  -  departs  from 
fluids  mechanics  tradition  in  two  important  respects.  First,  this  new  theory 

j  focuses  on  the  equilibrium  of  a  finite-size  flow  system  (e.g.,  jet,  wake, 

i 

■  boundary  layer),  as  opposed  to  the  equilibrium  of  the  infinitesimally  small 

|  fluid  packet  used  as  subject  in  the  Navi er-Stokes  equations.  Second,  it 

j  relies  on  the  thermodynamics  inspired  idea  that  the  flow  system  possess 

'  conservative  mechanical  properties  (e.g.  elasticity  in  bending)  if  it  is 
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inviscid. 

The  analytical  conclusion  of  the  buckling  theory  is  that  an  inviscid 
stream  of  finite  transversal  dimension  D  can  buckle  into  a  sinusoidal  shape 
whose  longitudinal  wavelength  Xg  scales  with  D.  In  this  manner  the  buckling 
theory  unveils  two  obvious  but  previously  unexplained  scales  of  the  stream, 
the  buckling  (meandering)  wavelength  Xg  and  the  time  scale  associated  with  the 
growth  of  the  deformation  (i.e.,  the  time  of  eddy  formation)  tg  “  Xg/(V/2), 
where  V  is  the  longitudinal  velocity  scale  of  the  flow.  This  new  theory  is 
therefore  a  theory  of  large  scale  orderly  structure  in  turbulent  flow. 

The  buckling  theory  -  its  enunciation,  experimental  verification  and 
application  to  heat  transfer  engineering  -  has  been  published  extensively 
through  three  equally  important  channels  of  peer  review:  journal  articles, 
listed  chronologically  as  Refs.  [4-15],  graduate  these  [16-20]  and  books 
[21-23].  The  permanent  contact  with  the  research  community  is  greatly 
responsible  for  the  diversity  of  the  topics  studied  from  the  point  of  view  of 
stream  buckling,  and  for  the  balance  between  the  theoretical  and  experimental 
segments  of  this  study. 

The  strategy  that  served  as  guide  for  the  definition  of  individual 
research  topics  during  the  four-year  period  consisted  of  identifying: 

I.  The  boundaries  of  the  new  theoretical  viewpoint,  i.e.,  how 
far  the  theory  extends  into  the  territory  that  is  not 
explained  by  any  existing  theory. 

II.  The  intersections  between  the  buckling  theory  and  existing 
theories,  i.e.,  pinpointing  the  phenomena  that  are  explained 
by  two  theories. 

If  the  domain  labeled  B  in  Fig.  1  represents  the  territory  discovered  until 
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today  from  the  point  of  view  of  8uckling  theory,  and  if  the  domain  S 
represents  the  things  explained  by  an  existing  theory  (e.g.,  hydrodynamic 
Stability  theory),  then  the  phenomenon  plotted  as  point  I  is  an  example  of 
work  done  on  front  (I).  Specific  examples  of  advances  of  type  I  are  the 
explanation  of  the  observed  sinusoidal  deformation  during  the  transition  to 
turbulence*  (Ref.  [22],  chapter  6),  the  Colburn  analogy  between  momentum  and 
heat  transfer  in  turbulent  flow  (Ref.  [22],  chapter  7)  and  the  "entrainment 
hypothesis"  or  the  linear  time-averaged  growth  of  all  turbulent  mixing  regions 
(Ref.  [22],  chapter  8). 

The  intersection  between  the  new  theory  and  an  existing  theory  is 
represented  by  points  labeled  II  in  Fig.  1.  One  discovery  of  type  II  is  the 
notion  that  the  transition  to  turbulence  in  flows  of  many  geometries  and 
origins  is  associated  universally  with  a  “local  Reynolds  number"  of  order 
0(102)  (Ref.  [22],  chapter  6;  the  "local  Reynolds  number"  is  based  on  the 
local  longitudinal  velocity  scale  and  the  local  transversal  length  scale). 
Although  this  apparently  universal  criterion  of  transition  to  turbulence  is 
consistent  with  both  theories  [10,11],  it  was  first  stated  in  the  realm  of 
buckling  theory  [4], 

As  a  second  example  of  an  advance  of  type  II  -  this  time  an  observation 
explained  by  both  S  and  B,  but  historically  stated  first  in  S  -  I  use  this 
opportunity  to  communicate  that  the  classical  Richardson  number  criterion 
regarding  the  instability  of  a  density-stratified  flow  follows  in  very  few 
steps  from  the  main  result  of  the  buckling  theory.  Consider  the  inviscid  jet 
of  transversal  length  scale  D  and  longitudinal  velocity  V  sketched  in  Fig.  2. 
According  to  the  buckling  theory,  since  the  stream  is  inviscid  it  has  the 

*  The  theory  of  hydrodynamic  stability  cannot  explain  this  observation: 
recall  that  in  stability  theory  the  "di sturbance"  shape  is  arbitrary  and 
assumed  given. 
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property  to  buckle  with  a  longitudinal  wavelength  Ag:  as  demonstrated  in 
chapter  4  of  Ref.  [21],  the  subsequent  and  certain  instability  of  the 
infinitesimally  buckled  stream  can  be  thought  of  as  the  effect  of  the  lateral 
net  force  FUD  that  acts  on  a  control  volume  (elbow)  of  height  0  and  length 


Ab/2, 


fud  ~  PV20a  , 


where  the  angle  a  is  infinitesimally  small. 

The  lateral  instability  and  eventual  roll-up  (eddy  formation)  may  not 
occur  if  the  flow  proceeds  through  a  stably  stratified  environment 

p(y)  =  po  -  by  ,  (2) 

where  b  =  -  dp/dy  >  0.  In  such  cases  the  finite-size  control  volume  will  be 
pushed  down  by  a  restoring  force  of  order 

Fdown  '  *B°9Ap  (3) 

where  Ap  ”  baAg  is  the  density  defect  scale  of  the  control  volume. 

The  condition  for  stability  can  then  be  written  sequentially  as 

Fdown  >  Fup  (4) 

AgDgbaAg  >  pV^Da  (5) 


-4-2  >  (£)  (6) 

4)  lB 

The  left-hand-side  in  this  inequality  is  the  Richardson  number.  According  to 
the  buckling  theory  the  right-hand-side  is  always  a  constant  of  order  0(1). 

In  conclusion,  the  stability  criterion  via  buckling  theory  is 

Ri  >  0(1)  (7) 
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The  same  criterion  is  a  classical  result  of  hydrodynamic  stability  theory, 
however,  the  analysis  summarized  as  eqs.  (1)  -  (7)  is  orders  of  magnitude 
simpler  than  the  corresponding  hydrodynamic  stability  analysis. 

In  conclusion,  in  Fig.  1  the  Richardson  number  criterion  (7)  is  covered 
by  theory  S  and  B:  this,  however,  does  not  mean  that  theory  S  alone  (or 
theory  B  alone)  explains  everything  about  transition  in  stably  stratified 
shear  flow.  For  example,  theory  S  still  does  not  explain  the  regular 
sinusoidal  deformation  that  serves  as  initial  condition  for  hydrodynamic 
i nstabil ity. 

The  work  started  on  fronts  I  and  II  is  far  from  over.  The  success  of 
buckling  theory  and  the  fact  that  none  of  the  claims  made  so  far  on  behalf  of 
this  theory  have  been  proven  false  in  the  post-1980  literature,  are  two  good 
reasons  why  this  theory  will  continue  to  attract  attention.  Considering  the 
naked  simplicity  of  buckling  theory,  and  recalling  the  thermodynamics  orgins 
of  the  essential  idea,  I  invite  my  colleagues  to  join  in  this  research  and  to 
give  new  life  to  the  words  spoken  once  by  Josiah  Willard  Gibbs: 

"One  of  the  principal  objects  of  theoretical  research  in 
any  department  of  knowledge  is  to  find  the  point  of  view 
from  which  the  subject  appears  in  its  greatest  simplicity." 


Adrian  Bejan 
Duke  University 
April  16,  1985 
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ABSTRACT 

This  letter  outlines  the  analogy  which  exists  between  inviscid  jets 
and  elastic  columns  in  axial  compression.  It  is  shown  that  straight 
inviscid  jet  columns  possess  the  property  of  sinusoidal  infinitesimal 
buckling.  The  buckling  wavelength  scales  with  the  transversal  dimen¬ 
sion  of  the  jet.  The  repeated  buckling  and  breakup  of  the  jet  column 
is  responsible  for  the  observed  whiplash  motion  of  turbulent  jets. 

The  buckling  theory  predicts  correctly  the  natural  frequency  of  the 
whiplash  motion  and  the  Reynolds  number  for  the  laminar-turbulent 
transition  in  free  jet  flow. 


Introduction 

One  aspect  of  jet  turbulence  receiving  increased  attention  is  the  large 
scale  "orderly"  structure  with  a  length  scale  of  the  same  order  as  the  jet  dia¬ 
meter.  Crow  and  Champagne  [l]  showed  that  as  the  Reynolds  number  increases  from 
2  3 

10  to  10  ,  the  whiplash  motion  of  jets  evolves  from  a  sinusoid  to  a  helix  and 
eventually  to  a  train  of  axisymmetric  waves.  Similar  observations  have  been  re¬ 
ported  by  Reynolds  [2].  The  characteristic  snake-like  shape  of  a  turbulent 
round  jet  is  shown  very  clearly  in  Fig.  1. 

The  purpose  of  this  letter  is  to  offer  a  theoretical  explanation  for  the 
observed  large-scale  periodic  structure  of  turbulent  jets.  The  explanation  is 
founded  on  a  very  interesting  analogy  which  exists  between  jet  flows  and  slender 
elastic  columns  in  axial  compression  [4].  It  may  be  recalled  that  Euler's 
theory  of  infinitesimal  buckling  (indifferent  equilibrium)  in  axially  compressed 
columns  rests  on  only  two  premises  [5,6]: 
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(i)  the  slender  column  is  straight  and  in  axial  compression; 

(ii)  if  subjected  to  a  separate  bending  test,  the  column  develops  in  its 
cross-section  a  resistive  bending  moment  which  is  proportional  to 
the  induced  curvature. 

Slender  elastic  columns,  of  course,  meet  these  two  conditions.  However,  it  is 
shown  in  the  next  section  that  exactly  the  same  conditions  are  met  by  inviscid 
streams  discharging  freely  into  larger  reservoirs.  Consequently,  the  column 
(control  volume)  occupied  by  an  inviscid  stream  buckles  sinusoidally,  and  the 
stream  mixes  periodically  with  the  stagnant  ambient.  From  this  result,  we  con¬ 
clude  that  the  natural  property  of  inviscid  streams  is  to  follow  a  sinusoidal 
(meandering)  path  as  they  travel  through  a  stagnant  ambient.  This  natural 
property  of  inviscid  flow  is  the  basis  for  the  unexplained  turbulent  behavior 
of  fluids. 


The  Static  Equilibrium  of  the  Jet  Envelope 

Consider  a  straight  inviscid  jet  of  density  p  ,  uniform  velocity  V  and 
cross-sectional  area  A  ,  as  shown  in  Fig.  2.  The  static  pressure  inside  the 
jet  and  in  the  ambient  fluid  is  Pq  .  Imagine  now  a  stationary  envelope 
(control  surface)  which  surrounds  a  certain  length  of  the  jet.  This  envelope 
and  two  transversal  end-cuts  define  a  stationary  control  volume. 

In  the  spirit  of  the  thermodynamics  of  open  flow  systems  [7],  the  only 
forces  which  act  on  this  control  volume  are  the  inlet  and  outlet  compressive 
forces 

C  =  p  A  V2  (1) 

Forces  C  are  shown  schematically  in  the  lower  half  of  Fig.  2  where  the  control 
volume  is  symbolized  by  the  solid  line.  At  this  point  we  conclude  that  the 
fluid-filled  column  represented  by  the  control  volume  satisfies  condition  (i) 
for  buckling. 

When  the  control  volume  is  slightly  curved,  each  face  of  the  transversal 
cut  is  exposed  not  only  to  a  compressive  force  but  also  to  a  bending  moment. 
Consider  the  separate  bending  test  in  which  the  je-  is  held  (forced  to  flow)  in 
a  slightly  curved  duct.  The  radius  of  curvature  of  the  duct,  R  ,  is  infinite¬ 
ly  greater  than  the  transversal  dimension  of  the  jet,  D.  Bernoulli's  equation 
for  a  streamline  [8]  dictates 

J  P  V2  +  Pq  =  j  p  V2  (z)  +  P(z)  (2) 


’  '  •  *  i  *  .  **  ‘  i  „  *•  ..  '  ^  „  1  „  ■>  ,  •  ,  '  ,  .  •  ‘  „ 
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where  the  right-hand-side  of  the  equation  corresponds  to  the  curved  duct. 
Coordinate  z  is  measured  radially  from  the  jet  centerline  towards  the  center 
of  curvature.  Radial  equilibrium  of  the  jet  fluid  in  the  curved  duct  requires 
also  [9] 

pv"  (z)  9P  .  . 

R  ~  3z  U 

CO 

In  the  limit  of  vanishingly  small  curvature,  D/R  -*  0  ,  equations  (2)  and  (3) 


v  (z)  =  V(1  +  — ) 


P(z)  =  P  - 
o 


The  bending  moment  acting  over  the  cross-section  is 

M  =  j  J  <pv2  +  P)  z  dA  =  2 


where  .1  =  Jj  z  dA  is  the  area  moment  of  inertia  of  the  jet  cross-section. 

A 

In  conclusion,  the  cross-sectional  bending  moment  M  of  a  nearly  straight 
jet  column  is  proportional  to  the  curvature  of  the  column,  1/R  .  This  means 

that  the  inviscid  jet  column  satisfies  condition  (ii)  for  the  infinitesimal 
buckling  of  a  straight  column. 

Consider  now  the  static  equilibrium  of  the  straight  jet  envelope  shown  in 
Fig.  2.  Clearly,  the  axial  compressive  forces  C  balance  each  other.  In 
addition,  the  excentricitv  bending  moment  CY  —  however  small  —  must  be 
balanced  at  all  times  by  the  cross-sectional  bending  moment  [10] 

CY  +  M  -  M(x)  =  0  (7) 

o 

2  2 

Noting  that  1/R^  =  -d  Y/dx  ,  the  rotational  equilibrium  condition  (7)  consti¬ 
tutes  a  differential  equation  for  the  equilibrium  centerline  Y(x)  , 

2 

pv2l  i-g-  +  pV2AY  +  M  =  0  (8) 

,  ^  o 

dx 


The  solution  satisfying  the  nozzle  conditions  Y  =  Y*  =  0  at  x  =  0  is 


(9) 
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where  the  amplitude  K  is  indeterminate  and  infinitely  small  compared  with  the 
transversal  length  scale  t'l/h  or  D  . 

We  arrive  at  the  important  conclusion  that  the  equilibrium  centerline  of 
the  jet  envelope  is  a  sinusoid  of  infinitely  small  amplitude.  The  wavelength 
of  this  trajectory  scales  only  with  the  transversal  dimension  of  the  jet, 

X  =  2it/l/A  (10) 

For  example,  \/D  =  tt/2  for  a  round  jet  of  diameter  D  ,  and  X/D  =  ir//J  for 
a  flat  (two-dimensional)  jet  of  thickness  D  .  An  analysis  of  the  meandering 
contour  of  the  round  jet  of  Fig.  1  yields  X/D  =  1.2  which  is  in  good  agreement 
with  the  theoretical  value  of  1.57  (it  is  very  likely  that  the  24%  discrepancy 
between  the  two  values  is  caused  by  the  fact  that  the  contour  shown  in  Fig.  1 
overestimates  the  real  diameter  of  the  turbulent  jet) . 


Formation  of  Large-Scale  Vortices 


The  analogy  between  the  infinitesimal  buckling  of  elastic  columns  and  in- 
viscid  jet  envelopes  terminates  with  the  equilibrium  centerline  given  by  equa¬ 
tion  (9) .  Whereas  in  slender  elastic  columns  the  equilibrium  of  the  small- 
amplitude  sinusoid  is  indifferent  [10],  in  an  inviscid  jet  the  equilibrium  is 
unstable.  The  unstable  equilibrium  of  the  inviscid  column  is  described  very 
well  by  the  classical  theory  of  hydrodynamic  instability  [  11 ] .  The  slightest 
deviation  of  the  jet  from  its  rectilinear  shape  leads  to  the  formation  of 
lateral  "lift"  forces  which  consistently  tend  to  amplify  the  deformation.  As  a 
result,  the  inviscid  jet  breaks  up  periodically  as  its  overextended  elbows 
penetrate  and  mix  with  the  stagnant  ambient.  The  degenerated  elbow  region 
becomes  a  "large  scale  turbulent  structure"  [l]  which  continues  to  move  down¬ 
stream  with  a  speed  of  order  V/2  . 

The  repeated  buckling  and  breakup  phenomenon  accounts  for  the  observed 
whiplash  and  fluctuating  motion  of  turbulent  jets.  The  period  of  this  fluctua¬ 
tion  scales  with  the  buckling  time  t  =  X/V  ,  i.e.  ,  with  the  time  of  fluid 

B 

travel  between  successive  elbows  (breaks)  in  the  jet  column.  The  natural  fre¬ 
quency  of  jet  fluctuation  can  be  expressed  in  dimensionless  form  as  a  Strouhal 
number 


which,  as  shown  in  the  preceeding  section,  has  a  value  of  order  0.5.  This  order 
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of  magnitude  estimate  agrees  very  well  with  experimental  measurements  of  the 
natural  frequency  of  turbulent  jets  exposed  to  a  range  of  external  excitation 
frequencies.  For  example,  Bechert  and  Pfizenmaier  [ 12 ]  reported  St  =  0.5  for 
maximum  amplification  of  broadband  jet  noise.  Most  recently,  Acton  [13]  con¬ 
ducted  a  computational  simulation  of  a  round  turbulent  jet  by  modeling  the  shear 
layer  as  a  succession  of  discrete  vortex  ring  elements.  Acton  showed  that  the 
natural  periodicity  of  the  round  jet  was  St  =  0.47  ,  and  that  the  jet  was  most 
sensitive  when  forced  at  approximately  the  same  frequency  (St  =  0.5). 


The  Transition  to  Turbulence 


The  buckling  property  and  fluctuating  nature  of  inviscid  jets  provide  a 
theoretical  basis  for  predicting  the  transition  to  turbulence.  The  transition 
from  laminar  jet  flow  to  turbulent  flow  occurs  when  the  stationary  ambient  is 
no  longer  capable  of  viscously  communicating  with  the  jet.  The  viscous  commu¬ 
nication  time  between  the  jet-ambient  interface  and  the  jet  centerline  follows 
from  Stokes'  first  problem  £  14 ] 


t 


v 


(12) 


The  jet  is  free  to  buckle,  i.e. ,  to  get  out  of  hand,  if  the  viscous  diffusion 
time  tv  is  longer  than  the  buckling  and  breakup  time  tB  .  Defining  the 
buckling  number  N  as  the  ratio  of  these  two  characteristic  times  of  the  flow 
configuration. 


N 


B 


(13) 


we  have  a  criterion  which  predicts  the  fluctuating  (turbulent)  behavior  of  the 

jet.  Using  equation  (12)  and  t  =  \/V  ,  the  buckling  number  can  be  written 

B 

also  as 


„  =  2™ 

B  V  16 


(14) 


For  a  round  jet,  the  transition  criterion 


mb  >  1 


becomes 


VD 

v 


> 


25 


(15) 


It  should  be  noted  that  this  order-of-magnitude  estimate  of  the  transition 
Reynolds  number  agrees  very  well  with  experimental  observations.  For  example, 
Viilu  [15]  found  the  value  of  11.2  for  the  Reynolds  number  of  breakdown  of  the 
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steady  laminar  jet.  Viilu's  observations  were  later  confirmed  by  Reynolds  [2] 
who  reported  the  range  10  <  Re  <  30  for  the  transition  Reynolds  number. 

The  buckling  number  criterion  of  transition  to  turbulence,  equation  (13) , 

explains  also  why  the  observed  transition  Re  is  a  number  considerably  greater 

than  unity.  For  a  dimensionless  group  to  truly  delineate  the  transition  from 

one  mechanism  to  another,  it  must  have  a  value  of  order  one  which,  after  all, 

reflects  the  balance  between  the  competing  mechanisms.  The  Reynolds  number  is 

not  the  correct  dimensionless  group  to  describe  transition  to  turbulence.  The 

Reynolds  number  is  the  experimental  (measurable)  reflection  of  the  N  transi- 

B 

tion  criterion. 


Concluding  Remarks 

This  letter  unveiled  a  theoretical  basis  for  predicting  the  large  scale 
structure  and  fluctuating  behavior  of  inviscid  jets.  The  analogy  between  jet 
envelopes  and  elastic  columns  in  axial  compression  showed  that  the  natural 
tendency  of  inviscid  jets  is  to  buckle  over  a  precise  wavelength  which  scales 
only  with  the  jet  diameter.  The  repeated  buckling  and  breakup  of  the  jet  column 
is  responsible  for  the  sinusoidal,  river-like,  path  of  turbulent  jets  and  also 
for  their  natural  whiplash  motion.  The  predicted  natural  frequency  of  the  jet 
agrees  very  well  with  measurements  from  harmonic  excitation  experiments. 

Finally,  we  learned  that  the  transition  from  laminar  to  turbulent  jet  flow 
occurs  when  the  buckling  and  breakup  time  is  shorter  than  the  time  of  viscous 
diffusion  across  the  jet.  This  last  conclusion  is  supported  strongly  by  experi¬ 
mental  observations. 

The  buckling  property  described  in  this  letter  emerges  as  the  fundamental 
property  serving  as  origin  for  turbulent  motion  in  inviscid  flows. 

Acknowledgement .  This  research  is  being  supported  by  the  Office  of  Naval 
Research. 
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Suleiman  and  Munson1  recently  presented  a  sequence 
of  interesting  experiments  concerning  the  buckling  of 
thin  viscous  fluid  layers.  In  their  paper  the  authors 
referred  to  Taylor’s  observations  on  the  buckling  of 
slender  viscous  filaments. 2  His  early  observations  and 
the  Suleiman  and  Munson  experiments  point  toward  an 
important  analogy  between  the  buckling  of  elastic  col¬ 
umns  and  the  buckling  of  viscous  filaments  in  longitudi¬ 
nal  compression. 

The  object  of  this  Comment  is  two-fold.  First,  it 
brings  to  the  readers’  attention  the  fact  that  a  buckling 
theory  of  thin  viscous  layers  already  exists.  Second, 
it  points  out  that  the  buckling  of  slender  columns  in  ax¬ 
ial  compression  is  not  a  property  only  of  elastic  solids 
and  highly  viscous  fluids,  but  also  a  property  of  irtvis- 
cid  columns  (streams). 

The  buckling  of  a  thin  viscous  layer  was  considered 
theoretically  by  Buckmaster,  Nachman,  and  Ting. 3 
They  showed  that  viscous  layers  in  longitudinal  com¬ 
pression  satisfy  two  basic  requirements:  (i)  the  vis¬ 
cous  layer  is  in  axial  compression  and  the  compressive 
force  is  proportional  to  the  relative  velocity  between 
the  two  ends  of  the  layer:  (ii)  if  curved,  the  viscous 
layer  develops  in  its  cross  section  a  bending  moment 
which  is  proportional  to  the  time- rate  of  change  in  the 
local  curvature.  The  momentum  equations  integrated 
over  the  viscous  layer  led  Buckmaster,  et  rtf.3  to  a 
global  equation  for  the  evolution  of  the  layer  centerline. 
This  equation  was  solved  by  assuming  various  initial 
disturbances  as  the  starting  shape  in  the  evolution  of 
the  viscous  layer.  In  a  subsequent  paper,  Buckmaster 
and  Nachman1  have  extended  this  theory  to  the  case 
where  surface  tension  effects  play  an  important  role. 


Of  fundamental  interest  are  the  similarities  between 
the  buckling  of  a  viscous  layer  (the  viscida  problem) 
and  the  buckling  of  a  slender  elastic  column  (the  elas- 
tica  problem).  A  close  examination  of  the  Euler  theory 
for  slender  elastic  columns  reveals  that,  for  infinitesi¬ 
mal  buckling,  elastic  columns  must  satisfy  only  two 
conditions  (Ref.  5,  p.  184):  (iii)  the  column  is  straight 
and  in  axial  compression;  (iv)  if  subjected  to  a  separate 
bending  test,  the  column  develops  a  resistive  bending 
moment  in  its  cross  section  which  is  proportional  to 
the  local  curvature.  The  infinitesimal  budding  of  the 
straight  column  follows  from  invoking  static  equilibrium 
for  the  system  sketched  in  Fig.  1.  For  a  straight  or 
nearly  straight  column  there  are  two  equilibrium  condi¬ 
tions  to  consider.  The  first  is  the  obvious  balance  of 
compressive  forces  C  in  the  longitudinal  direction. 

The  second  is  one  of  rotational  equilibrium  which  states 
that  the  eccentricity  bending  moment  Cy,  however 
small,  must  be  balanced  at  all  times  by  the  cross-sec¬ 
tional  bending  moment  M.  In  the  case  of  a  column  con¬ 
taining  elastic  material,  a  special  bending  test  of  pre¬ 
scribed  curvature  (Ref.  5,  p.  37)  combined  with  knowl¬ 
edge  of  the  elastic  properties  of  individual  fibers  in 
the  slender  column  leads  to  the  notion  of  a  cross-sec¬ 
tion  bending  moment  proportional  to  the  local  curvature 
(iv).  Combining  this  notion  with  the  rotational  equilib¬ 
rium  condition  yields  the  equilibrium  shape  of  the  straight 
column,  namely,  a  sinusoid  of  infinitely  small  amplitude. 
The  infinitesimal  amplitude  is  indeterminate,  hence,  the 
equilibrium  of  the  nearly  straight  elastic  column  is  indiffer¬ 
ent.  If  present,  the  slightest  lateral  force  is  able  to  push  the 
column  away  from  the  straight  equilibrium  shape  into  a 
nearly  straight  (sinusoidal)  equilibrium  shape  (Ref.  5, 
p.  184). 
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FIG.  1.  (a)  Translational  and  rotational  equilibrium  of  a 
slender  column  ( control  volume)  in  axial  compression,  (b) 
The  natural  buckling  wavelength  of  circular  and  two-dimen¬ 
sional  inviscid  jets. 


In  this  discussion  it  is  important  to  include  the  case 
of  inviscid  fluid  columns:  Do  such  columns  conform  to 
buckling  of  type  (i,  ii)  or  to  buckling  of  type  (iii,  iv)? 
Thermodynamic  reasoning  alone  suggests  that  inviscid 
columns  should  buckle  according  to  model  (iii,  iv)  be¬ 
cause,  like  elastic  solids  and  unlike  highly  viscous  flu¬ 
ids,  inviscid  fluids  are  free  of  entropy  generation. 8 

Referring  again  to  Fig.  1,  we  note  that  a  column  con¬ 
taining  inviscid  fluid  is  none  other  than  the  imaginary 
control  surface  drawn  around  an  inviscid  jet  flowing 
through  an  inviscid  fluid  at  rest.  Since  the  column  is 
stationary,  its  equilibrium  is  described  by  the  two 
static  conditions  (force  and  moment)  discussed  in  con¬ 
nection  with  the  Euler  buckling  of  elastic  columns. 

The  remaining  problem  is  to  determine  what  special 
forms  C  and  M  take  in  the  case  of  inviscid  stream  col¬ 
umns.  For  the  compressive  force  on  the  control  vol¬ 
ume,  C,  we  know  from  the  thermodynamics  of  open 
flow  systems  that7 

C=pAV\  (1) 

where  p,A,  V  are  the  density,  cross-sectional  area, 
and  velocity  of  the  stream.  In  order  to  determine  the 
cross-sectional  bending  moment  M,  we  conduct  a  spe¬ 
cial  bending  experiment  where  the  inviscid  stream  is 
held  in  a  duct  of  known  radius  of  curvature,  Rm.  In  the 
limit  of  vanishingly  small  curvature,  the  inviscid  flow 
equations  show  that  the  resistive  bending  moment  in  the 
cross  section  is8 

M  =  p\*l/R„,  (2) 

where  /  is  the  area  moment  of  inertia  of  the  cross  sec¬ 
tion  (Ref.  5,  p.  37).  Equations  (1)  and  (2)  show  that 
inviscid  fluid  columns  obey  conditions  (iii,  iv)  for  in¬ 


finitesimal  buckling.  The  corresponding  equilibrium 
shape  of  the  column  is  a  sinusoid  of  infinitely  small 
amplitude,8  whose  wavelength  is 

X  =  2  t(l/A)ul.  (3) 

Inviscid  streams,  like  elastic  rods  and  viscous  lay¬ 
ers,  possess  the  natural  property  of  buckling.  The 
wavelength  of  the  buckled  shape  scales  only  with  the 
transverse  dimension  of  stream;  for  example,  Eq.  (3) 
for  a  jet  of  round  cross  section,  yields  \/D  =  it.  2.  Un¬ 
like  buckled  elastic  columns,  whose  equilibrium  is  in¬ 
different  (Ref.  5,  p.  184),  buckled  inviscid  columns 
are  unstable.  The  post-buckling  evolution  of  the  invis¬ 
cid  column  is  well  understood,  forming  the  subject  of 
the  classical  theory  of  hydrodynamic  stability. 9  How¬ 
ever,  the  buckling  property  is  to  be  recognized  as  re¬ 
sponsible  for  the  wave-like  “disturbance"  assumed 
routinely  ( empirically )  as  a  starting  point  in  any  hydro- 
dynamic  stability  analysis. 

In  essence,  the  buckling  property  of  inviscid  streams 
guarantees  that  such  streams  cannot  flow  straight 
through  another  fluid  or  through  a  flexible  duct.  This 
new  property  of  inviscid  fluids  stands  at  the  very  root 
of  the  phenomenon  of  turbulence .  For  example,  the  in¬ 
teraction  between  the  elbows  of  the  buckled  stream  and 
the  stagnant  medium  leads  to  the  periodic  formation  of 
large  eddies.  These  large-scale  turbulent  structures 
are  responsible  for  the  river- like  shape  of  turbulent 
jets,  wakes,  and  plumes:  It  is  widely  observed  that  the 
wavelength  of  this  large-scale  meandering  path  scales 
with  the  diameter  of  the  stream,  as  predicted  in  Eq. 

(3).  The  natural  buckling  of  inviscid  streams  is  par¬ 
ticularly  visible  in  the  case  of  streams  flowing  within 
flexible  boundaries,  like  rivers  and  capillary  rivulets 
at  high  Reynolds  numbers.  The  buckling  property 
serves  as  the  theoretical  basis  for  predicting  other 
turbulence  parameters,  for  example,  the  natural  fre¬ 
quency  of  turbulent  jets  (the  Strouhal  number)  and  the 
critical  Reynolds  number  for  transition  to  turbulence.8 
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This  paper  discusses  experimental  observations  of  the  meandering  fall  of  light-weight  tissue  paper  ribbons 
The  photographs  show  that  the  ribbons  assume  a  sinusoidal  shape  with  a  unique  wavelength  which  scales 
with  the  thickness  of  the  airstream  entrained  by  the  ribbon 


The  objective  of  this  paper  is  to  present  a  series  of 
interesting  experimental  observations  concerning  the 
meandering  motion  executed  by  highly  flexible  ribbons 
falling  through  the  air.  The  experiment  consisted  of 
dropping  a  length  of  light-weight  toilet  tissue  paper 
through  the  air  and  photographing  its  shape  as  it  falls 
to  the  ground.  The  reader  may  take  note  of  the  fact  that 
this  falling- ribbon  phenomenon  occurs  naturally  when 
excited  sports  fans  launch  rolls  of  tissue  paper  from  the 
stands  onto  the  playing  field.  Another  natural  phenome¬ 
non  related  to  the  falling-ribbon  experiments  described 
in  this  paper  is  the  "waving  of  flags’’1  and  the  "vibra¬ 
tion”  of  tape  drives  used  in  the  computer  technology.2 
The  classical  prespective  in  the  study  of  flag  waving 
falls  in  the  realm  of  hydrodynamic  stability  theory, 
where  one  questions  the  stability  of  the  flexible  solid 
surface.  The  starting  point  in  the  stability  study  is  the 
assumption  of  an  initial  deformation  of  arbitrary  wave¬ 
length. 

In  the  present  experiments,  ribbons  of  various  lengths 
were  dropped  from  heights  in  the  range  3-1  m,  through 
the  quiescent  air  of  the  laboratory.  The  time  of  free 
fall  was  measured  with  a  digital  stopwatch;  it  was  found 
that  the  ribbon  reached  its  terminal  velocity  very  quick¬ 
ly,  therefore,  the  free-fall  velocity  U  could  be  deter¬ 
mined  by  dividing  the  total  travel  by  the  measured  time 
of  free-fall.  In  order  to  force  the  ribbon  to  fall  "head 
first,”  one  end  was  loaded  with  a  lead  refill  for  a  mech¬ 
anical  pencil. 

This  simple  experiment  was  repeated  many  times  and, 
in  all  cases,  the  photographs  showed  that  the  falling 
ribbon  acquires  a  sinuous  shape;  the  wavelength  of  this 
shape  was  the  same  for  all  the  cases  involving  a  ribbon 
of  fixed  length.  Figures  1(a)  and  1(b)  show  very  clearly 
the  characteristic  sinuous  shape  observed  in  these  ex¬ 
periments.  A  longer  ribbon  (Fig.  lib)]  exhibits  a  rela¬ 
tively  longer  wavelength. 

Another  important  observation  is  the  fact  that  the  sin¬ 
uous  shape  travels  as  a  solid  body  downward,  at  a  speed 
of  order  U,  2,  where  U  is  the  speed  of  the  tissue  paper 
itself.  The  wave  speed  was  measured  photographically, 
as  shown  in  Fig.  2.  This  photograph  was  obtained  with 
the  shutter  open  in  complete  darkness,  while  lighting 
the  falling  ribbon  with  a  strobe  light  three  times,  at 
precise  time  intervals  (I  =0.025  sec).  The  fact  that  the 
sinuous  shape  travels  at  half-speed  is  strong  evidence 
that  the  sinuous  shape  is  produced  not  by  the  tissue  pa¬ 
per,  but  by  the  airstream  whose  centerline  moves  at 
top  speed  U  through  an  ambient  at  rest. 


The  key  measurement  facilitated  by  the  falhng-nbbon 
experiment  is  that  of  the  meander  wavelength.  Measur¬ 
ing  the  distance  between  the  elbows  of  the  sinuous  shape 
of  Figs.  1(a)  and  1(b),  and  averaging  these  measure¬ 
ments  over  the  sinuous  portions  of  each  ribbon,  yields 
the  wavelengths  listed  under  hB  m  Table  I.  The  rela¬ 
tively  small  standard  deviations  of  these  measurements 
indicate  that  the  elbow-to-elbow  distance  does  not  vary 
appreciably  along  the  wavy  portion  of  the  ribbon. 

The  effective  thickness  Ds  of  the  air  stream  entrained 
by  the  ribbon  can  be  calculated  based  on  the  following 
energy-conservation  argument.  During  its  steady  fall  at 
terminal  velocity  U ,  the  ribbon  weight  if7  performs  the 
mechanical  work  VFAL  on  its  ambient;  AT  is  the  linear 
increment  in  downward  travel,  equal  to  f/Af,  where  At 


FIG.  1.  I al  Fall  of  a  1.25  m  ribbon  of  tissue  paper.  >hl  Fall 
of  a  1,83  m  ribbon  of  tissue  paper. 
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TABLE  I.  Laboratory  measurements  of  the  meandering  wavelength  of  the  inviscid  air  stream 
driven  by  a  falling  ribbon. 


Experiment 

Ribbon  dimensions 
length  x  width 

Lx  ft  lcm:) 

Total  weight 
tV(g) 

Terminal 

speed 

C/(  m/s) 

Meander 
wavelength 
\B  (cm) 

Effective 

air-stream 

thickness 

fl,(cm) 

Fig.  1(a) 

125x11.4 

3.21 

2.  93  t  0. 15 

15.5  ±1.3 

5.5  i  0,  56 

Fig.  1(b) 

183x11.4 

4.  97 

3.  04  t  0.09 

25  ±1.5 

7.  S6  ±  0.46 

is  the  time  increment.  The  work  done  by  the  weight  is 
first  converted  into  the  kinetic  energy  imparted  to  the 
airpacket  pierced  by  the  tip  of  the  ribbon  during  the 
time  &t  (at  the  same  time,  the  ribbon-air  train  sheds  a 
moving  air  packet  of  the  same  size:  the  kinetic  energy 
of  this  air  packet  is  eventually  dissipated  in  the  wake). 
Equating  the  two  energy  increments,  we  write 

W±L  =(pD,±Lb)(U\  2)  or  £>,  =  2(«'  pi,2h) ,  (1) 

where  b  is  the  ribbon  width  and  p  is  the  air  density. 

The  results  of  this  calcu  .tion  are  listed  in  Table  I. 
Clearly,  the  meander  wavelength  scales  with  the  air 
stream  thickness. 

A  possible  explanation  for  the  above  observations 
may  be  offered  based  on  the  buckling  theorv  of  inviscid 
streams.114  One  key  result  oi  this  theorv  is  the  univer¬ 
sal  proportionality  which  must  exist  bet  .wen  stream 


thickness  (D)  and  "buckling”  wavelength^a).  For  a 
two-dimensional  stream  one  finds  0=(v  3<  *).\8,  which 
agrees  in  an  order  of  magnitude  sense  with  the  meas¬ 
urements  listed  in  Table  I.  Additional  evidence  sup¬ 
porting  this  explanation  is  the  fact  that,  from  Fig.  1(a) 
to  1(b),  both  >s  and  Dt  increase. 
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Abstract  ■  This  paper  advances  a  theoretical 
explanation  for  the  laceral  periodicity  and  geo¬ 
metric  similarity  of  meanders  observed  in  rivers  of 
many  sizes.  Invoking  the  static  equilibrium  of  a 
- traight  river  bed,  it  is  shown  analytically  that 
the  equilibrium  shape  of  the  bed  is  a  unique 
sinusoid.  The  theoretical  wavelength  of  the 
sinusoidal  shape  is  proportional  to  the  width  of 
the  river,  in  agreement  with  visual  observations 
of  rivers  of  all  sizes. 

1.  The  Geometric  Similarity  of  River  Meanders 

The  most  basic  feature  of  the  meander  phenomen¬ 
on  is  the  proportionality  between  meander  wave¬ 
length  X  and  stream  width  W  (Leopold  and 
Wolman  1960).  As  concluded  in  a  recent  review 
article  (Callander  1978),  significant  theoretical 
progress  has  been  made  in  Che  direction  of  ex¬ 
plaining  the  occurrence  of  meanders  and  accounting 
for  measured  meander  parameters.  The  theories 
developed  over  the  past  thirty  years  have  in 
common  the  thinking  framework  offered  by  the 
Theory  of  Hydrodynamic  Stability:  according  to 
this  approach,  one  analyzes  the  stability  or  in¬ 
stability  of  sinusoidal  waves  (disturbances) 
superimposed  on  Che  straight  river  flow.  As 
summarized  by  Parker  (1976),  consideration  of 
gravity  waves  led  Werner  (1951)  to  the  conclusion 
that  X/W-2F,  where  F  is  the  flow  Froude  number. 
H.3nsen  (1967)  posed  the  linearized  stability 
problem  associated  with  a  meandering  (sinusoidal) 
disturbance  and  found  X/H  »  7  F"/S,  where  H  is 
the  river  depth  and  S  is  the  dip  angle  of  the 
river  bed.  Anderson  (1967)  analyzed  transverse 
oscillations  and  obtained 

\/W  =  constant  x  (F  H/W)4.  A  number  of  improved 
models  of  hydrodynamic  stability  have  been  re¬ 
ported  more  recently:  for  example,  Callander 
(1969),  Sukegawa  (1970),  Hayashi  (1974),  Enge- 
lund  and  Skovgaard  (1973)  and  Parker  (1976). 

The  chief  contribution  of  the  theoretical 
research  reviewed  above  and,  in  greater  detail, 
in  Callander  (1978),  is  to  have  shown  that  the 
bed  of  a  channel  with  straight  banks  is  unstable: 
since  the  bed  is  composed  of  moving  sediment, 
the  amplitude  of  a  certain  class  of  disturbances 
is  likely  to  grow  the  fastest.  In  addition, 
this  research  predicted  correctly  the  direction 
of  migration  and  the  downstream  wave  speed  of 
meanders.  These  conclusions  are  important  and 
will  be  adopted  without  debate  in  the  theory 
developed  in  this  paper.  However,  the  chief 
limitation  of  the  existing  theories  is  also 
important:  note  should  be  made  of  the  fact  that 
the  ratio  '/W  predicted  by  hydrodynamic 
stability  considerations  (Table  1)  does  not  have 

Copyright  1982  by  the  American  Geophysical  Union. 

Paper  number  2L1032. 

0094-82 76 /82/002L- 1032$ 3. 00 


a  characteristic  (unique)  range  of  values, 
contrary  to  the  statistical  evidence  compiled 
over  natural  streams  of  widely  varying  sizes 
(Leopold  and  Wolman  1960).  This  limitation  is 
accepted  from  the  start  by  every  stability 
analyst  who  postulates  the  existence  of  distur¬ 
bances  of  every  wavelenarh  and  then  tries  to 
identify  which  of  these  disturbances  will 
develop  the  most  rapidly. 

Inasmuch  as  the  constancy  of  the  ratio  \/W 
appears  to  be  an  intrinsic  property  of  all 
streams,  the  stability  theories  explain  the  evolu¬ 
tion  (behavior)  of  exist ing  meanders  but  not  their 
origin. 

The  theory  constructed  in  this  paper  addresses 
the  question  which  has  not  yet  been  addressed, 
namely,  why  does  a  straight  (undisturbed)  river 
choose  a  sinuous  shape  of  precise  wavelength?  Be¬ 
fore  presenting  the  theoretical  answer  to  this 
fundamental  question,  the  author  finds  it  neces¬ 
sary  to  review  a  class  of  very  useful  experiments 
which  shed  light  on  the  natural  properties  of  the 
(river)-(f lexible  bed)  system. 

2.  The  Stream  Plate  Experiments 

The  original  stream  plate  experiment  was 
proposed  by  Tanner  (1960)  as  a  means  of  visual¬ 
izing  the  natural  tendency  of  streams  to  meander. 
The  same  technique  w  s  used  in  a  comprehensive 
study  by  Gorycki  (1973a).  The  experiment  consists 
of  a  smooth  plane  surface  which  supports  a  water 
jet  flow  issuing  from  a  nozzle  tangent  to  the 
surface.  Thus,  the  stream  plate  experiment  is  the 
laboratory  version  of  water  tricklings  commonly 
observed  on  shower  walls  and  car  windshields.  The 
water  columns  generated  in  this  fashion  meander  in 
much  the  same  way  as  rivers  do.  However,  there  is 
one  important  advantage  to  this  experiment,  name¬ 
ly,  the  opportunity  to  observe  the  incipient  phase 
of  the  meandering  process  in  a  straight  jet. 

Of  interest  here  are  Tanner's  and  Gorycki' s 
observations,  and  the  clever  mechanical  analog 
visualized  by  Gorycki  to  simulate  river  meander¬ 
ing.  First,  the  stream  plate  visualization  of  the 
meander  formation  process  is  strong  evidence  that 
the  meander  is  a  property  of  the  stream;  this 
property  is  independent  of  the  effects  of  sedi¬ 
ment.  In  comparing  various  stream  plate  experi¬ 
ments,  Gorycki  (1973b)  argues  further  that  plate 
toughness  is  not  necessary  for  meander  formation. 

The  experiments  also  showed  (Tanner  1962)  that 
meandering  Is  independent  of  the  secondary  flow 
or  disturbances  which  may  be  present  in  the 
nozzle.  This  conclusion  is  strengthened  by  the 
statement  made  by  Schumm  and  Khan  (1972)  who  ob¬ 
served  meanders  made  in  the  laboratory  under 
straight  entrance  conditions:  "...a  perturbation 
or  disturbance  of  the  flow  may  not  be  an  essen¬ 
tial  cause  of  meandering"  (Gorycki  1973a,  p.  178). 

Gorycki  presented  also  a  mechanical  analog  of 
the  river  meandering  mechanism:  he  held  a  slender 
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TABLE  1.  Comparison  ot'  the  Present  Theory  with  Meander  Wavelength  Observations 


l—  - 

Reference 

A/W 

Remarks 

present  the. ry 

2 . 11 

incipient  buckling 
(absolutely  straight  stream) 

Gorycki  (1973a) 

2-3 

stream  plate  simulations 

i 

} 

Leopold,  Wolman  and 

Miller  (1964), 

Leopold  and  Wolman  (1970), 

Dury  (1964) 

?-3 

field  studies  in  straight 
rivers 

i 

Schumm  and  Khan  (1972) 

Table  3 

3.24  r  0.64 

laboratory  channel 
during  development 

Keller  (1972) 

3-5 

channel  experiments 
during  development 

l 

i 

i 

5-7 

during  late  stages  of 
development 

! 

Leopold  and  Wolman  (1960) 

Table  1 

6.5  -  11 

natural  (long  history) 

j _ L 


piece  of  elastic  (a  slender  cylindrical  column  of 
rubber)  between  two  parallel  pieces  of  glass.  He 
then  moved  the  glass  pieces  relative  to  one  an¬ 
other,  in  the  direction  parallel  to  the  cylinder 
axis.  As  a  result,  the  piece  of  elastic  assumed 
a  shape  which  resembles  very  closely  that  of  a 
meandering  scream  in  a  stream  plate  experiment. 

It  is  shown  in  the  next  section  that  the 
mechanism  responsible  for  elastic  meandering  is 
also  responsible  for  river  meandering.  In  Mechan¬ 
ical  Engineering,  which  is  the  present  author's 
education,  the  sinusoidal  shape  of  the  elastic 
column  has  been  explained  as  the  buckling  property 
of  slender  spaces  in  longitudinal  compression. 

The  buckling  of  slender  columns  represents  an 
important  and  voluminous  chapter  in  the  centuries- 
old  discipline  of  Strength  of  Materials  (Den 
Hartog,  1961).  The  theoretical  basis  for  this 
chapter  was  established  by  Euler,  who  pointed  out 
that  in  order  for  a  slender  space  to  buckle  into 
an  equilibrium  sinusoidal  shape  it  must  satisfy 
only  two  conditions  (Love  1927) : 

(i)  the  slender  space  must  be  in  a  state  of 
axial  (longitudinal)  compression. 

(ii)  the  material  which  fills  the  space  must 
be  such  that  if  the  space  is  subjected 
to  a  separate  bending  test  of  prescribed 
curvature,  then  the  space  develops  in  its 
cross-section  a  resistive  bending  moment 
(couple)  which  is  proportional  to  the 
local  curvature. 

The  piece  of  elastic  described  by  Corycki 
satisfies  the  above  conditions.  But  it  is  impor¬ 
tant  to  keep  in  mind  that  conditions  (i)  and  (ii) 
do  not  refer  to  a  specific  material  such  as  an 
elastic  solid:  they  refer  to  a  space  (column)  of 
finite  sice.  The  objective  of  the  following 
analysis  is  to  prove  that  the  river  and  its  bed, 
as  a  sLender  space,  also  satisfy  conditions  (i) 
and  (ii)  necessary  for  sinusoidal  buckling. 


3.  The  Static  Equilibrium  of  a  Straight  Kiver 

Consider  a  straight,  inviscid,  river  flow  of 
uniform  velocity  V,  density  c,  width  W  and 
depth  H.  The  bed  of  the  river  is  horizontal. 

In  the  stationary  frame  of  reference  of  the  bed, 
the  river  cross-section  exhibits  a  uniform  com¬ 
pressive  stress  £>v02  +  P0(z),  where  P0  is 
the  excess  pressure 

P„(z)  “  pg  (H  -  z)  (1) 

and  g  is  the  gravitational  acceleration.  The 
resultant  of  this  compressive  stress,  integrated 
over  the  river  cross-section,  is 

C  =  pv  8 :  WH  (1  +  )  .  (2) 

Therefore,  in  the  frame  of  reference  of  the  bed, 
the  straight  river  is  a  slender  column  in  longi¬ 
tudinal  compression.  This  means  that  condition 
(i)  is  satisfied.  In  order  to  see  that  the  bed, 
as  a  duct,  is  in  a  state  of  longitudinal  compres¬ 
sion,  the  reader  should  think  of  a  piece  of  garden 
hose  through  which  the  flowrate  is  high  (turbu¬ 
lent).  The  hose  is  pushed  axially  by  the  reactive 
compressive  force  associated  with  the  stream 
leaving  through  the  open  end.  The  same  hose  is 
pushed  axially  in  the  opposite  direction  by  the 
impactive  force  associated  with  the  stream  enter¬ 
ing  it  from  the  faucet. 

As  shown  in  Fig.  1,  the  static  equilibrium  of 
the  bed  requires  two  statements,  one  for  transla¬ 
tional  equilibrium  (obvious,  CSC)  and  the  other 
for  rotational  equilibrium.  The  rotational  equil¬ 
ibrium  condition  is  made  necessary  by  the  fact 
that  no  straight  river  is  ever  "mathematically" 
straight;  in  other  words,  all  straight  rivers  are 
subjected  to  an  infinitesimally  small  couple  CY 
due  to  the  imperfect  colinearity  of  the  axial 
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forces  C.  The  rotational  equilibrium  condition 
is  (Den  Hartog  1961) 

CY  -  M  +  M,.  -  0  (J) 

where  >1  is  the  ooupie  acting  over  the  river 
cross-section. 

It  is  easy  to  show  that  cue  river  column  also 
satisfies  condition  (ii),  or  chat  a  nec  couple  M 
is  present  whenever  the  column  Y (X)  is  locally 
curved.  Consider  a  separate  bending  test  in 
which  the  river  radius  of  curvature  R  is  infin¬ 
itely  greater  chan  the  river  width  U  (Fig.  2). 
The  new  velocity  V  and  pressure  distribution  P 
in  the  cross-section  car.  be  determined  immediate¬ 
ly  from  Bernoulli's  equation  (Prandtl  1969) 


Fig.  2.  Cross-section  through  the  river. 


proves  that  rivers  also  obey  condition  (ii)  for 
sinusoidal  buckling. 

4.  The  Natural  Sinusoidal  Shape  of  Rivers 


|  pV2  +  P 


combined  with  a  local  force  balance  in  the  radial 
direction  y  (Prandtl  1969), 

=  (5) 

R»  3y 

In  the  limit  of  vanishingly  small  curvature, 

W/R^  -*■  0,  we  obtain 

V  =  V0  (1  -  )  (6) 


In  the  limit  of  infinitesimally  snail 
deviations  from  the  rectilinear  shape,  the  river 
curvature  1/R^  is  equal  to  -d2Y/dX-.  Based  on 
this  approximation  and  equation  (9),  the  static 
equilibrium  condition  (3)  becomes 


CY  +  |  pV02  HW3Y"  +  M0  -  0 


The  general  equilibirum  shape  of  the  river  bed 
follows  from  equation  (10), 


™  0  y 

V(X)  -  —  [cos(24  )  -  1], 


P  -  cV,2  |  +  jg(H  — z)  ■ 


A  related  result  is  that  the  free  surface  z^(y) 
acquires  a  slight  tilt 


Zf(y)  -  h  +i5- 


where  the  wavelength  \  has  a  precise  value  given 
by 

A  .  - H - _  (12) 

[2  +  gH/V.i]1* 


For  shallow  rivers  and  for  stream  plate  simula¬ 
tions  (gH/V q 2  «  l)  we  find 


Due  to  the  slight  tilt,  the  center  of  mass  of  the 
cross-section  shifts  from  y  =  0  to 
V 3  -  V32W2/(12gHR00),  where  y.«  W. 

According  to  the  standard  methods  of  mechanical 
engineering,  the  net  bending  moment  M  about  the 
vertical  line  passing  through  the  center  of  mass 
of  the  cross-section  is 


y-W/2  z-zf(y) 


fr 

M  ”  (pV2+P)  (y-yo)  dzdy 


pV,-HW3 


v«-W/2  z-0 


This  result  demonstrates  that  the  river  cross- 
section  experiences  a  bending  moment  as  soon  as 
the  river  trajectory  has  curvature.  Equation  (9) 


h  ■  —  =  2.22,  universal  constant 

W  rr 


In  conclusion,  the  natural  (equilibrium)  shape 
of  the  river  bed  is  a  sinusoid  whose  wavelength 
is  a  precise  multiple  of  the  river  width.  The 
amplitude  of  this  shape  is  unknown  (infinitely 
small),  because  the  analysis  leading  to  equation 
(12)  invoked  the  static  equilibrium  of  a  straight 
river.  It  is  well  known,  however,  that  the 
equilibrium  of  a  nearly  straight  river  bed  is 
unstable  (cf.  Hydrodynamic  Stability  Theory, 
Section  1)  and  that  the  highly  regular,  sinuous, 
shape  determined  here  is  destined  to  grow  in 
amplitude. 

5.  Conclusion 


^  V 


Fig.  1.  Static  equilibrium  of  a  straight  river 
bed  (view  from  above). 


The  fundamental  contribution  of  this 
theoretical  argumenc  is  the  prediction  of  a 
universal  proportionality  between  meander  wave¬ 
length  and  river  width,  equations  (12,  13).  This 
prediction  is  supported  strongly  by  observations 
of  meander  formation  in  straight  rivers  and  in 
stream  plate  experiments.  Attention  is  drawn 
first  to  the  stream  plate  experiments  (Gorycki 
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;y/3.i,  p.  L'9,  Figs.  3,  8,  9,  10,  13)  which 
conclusively  show  "chat  the  spacing  between  the 
sinuous  curves  or  point  bars  in  straight  streams 
is  approximately  two  to  three  tines  the  stream 
width",  as  in  equation  (13).  Similar  values  of 
■  /W  in  straight  rivers  were  reported  by  numerous 
field  studies;  for  example,  in  Leopold.  Woiman  and 
Miller  (196-.),  Leopold  and  Woiman  (1970,  Fig.  7.S) 
and  Dury  (1969,  Figs.  26,  28).  A  ratio  3/W  be¬ 
tween  2  and  3  appears  to  be  a  universal  feature  of 
all  straight  streams.  Keller  (1972,  p.  1534) 
showed  that  only  in  the  late  stages  of  meander 
development  \/W  reaches  values  in  the  range  of 
5-7,  whereas  during  development  the  ratio  \/W  is 
in  che  range  3-5. 

Table  1  shows  a  summary  of  experimental  obser¬ 
vations  next  to  the  constant  ratio  >/W  predicted 
by  the  presenc  theory  of  river  buckling.  The 
present  theory  agrees  with  observations  made  in 
straight  or  nearly  straight  rivers,  which  conform 
to  che  type  of  system  selected  here  for  analysis 
in  Section  2.  Table  1  shows  also  that,  as  time 
passes,  che  ratio  •  /W  increases  from  the  initial 
theoretical  value  (2.22)  to  the  natural  (long- 
history)  value  of  approximately  10.  Thus,  the 
present  theory  offers  a  concise  explanation  for 
the  origin  of  meanders  in  straight  rivers. 
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The  nonaxisymmetric  (buckling)  flow  regime  of  fast  capillary  jets 
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This  paper  reports  an  experimental  study  of  the  nonaxisymmetric  flow  of  a  fast  liquid  jet 
discharging  into  the  atmosphere.  The  nonaxisymmetric  shape  of  the  jet  was  photographed  and 
subjected  to  a  wavelength  analysis.  The  results  of  the  wavelength  analysis  demonstrate  that  the  jet 
shape  is  governed  by  a  narrow  band  of  wavelengths  associated  with  a  characteristic  value /lm.. 
which  scales  with  the  jet  diameter  D.  It  is  shown  that  the  experimental  observations  are  in 
agreement  with  predictions  based  on  hydrodynamic  stability  theory  and  buckling  theory. 


I.  INTRODUCTION 

The  problem  of  capillary  jet  flow  and  breakup  has  a 
long  history  beginning  with  the  qualitative  studies  of  Bi- 
done1  and  Savart,2  which  were  extended  by  Savart,  Plateau, 
and  Rayleigh,  and  summarized  later  by  Rayleigh.3  These 
studies  focused  on  the  low-speed  regime  where  the  jet  forms 
radially  symmetric,  regularly  shaped,  drops  of  measurable 
frequency.  Rayleigh  studied  the  symmetric  (varicose)  break¬ 
up  theoretically,  by  imposing  hypothetical  infinitesimal  dis- 
trubances  on  the  jet  and  examining  the  stability  (or  instabil¬ 
ity)  of  each  disturbance  in  time.  Rayleigh’s  theory  was  sum¬ 
marized  and  extended  in  several  directions  by 
Chandrasekhar.4 

Much  of  the  post-Rayleigh  work  focused  almost  exclu¬ 
sively  on  the  axisymmetric  Ivaricose)  regime,  although  there 
have  been  a  number  of  instances  in  which  a  nonaxisymme- 
tnc  breakup  mode  was  observed.  Crane,  Birch,  and  McCor¬ 
mick3  employed  an  electronically  driven  vibrator  to  study 
the  dispersion  curve  of  low-speed  jets  (the  dispersion  curve  is 
the  graphical  representation  of  the  response  of  the  capillary 
jet  to  a  continuous  set  of  imposed  disturbance  frequencies). 
Their  results  agreed  only  qualitatively  with  Rayleigh’s.  A 
similar  experiment  was  described  by  Donnelly  and  Glaber- 
son,6  who  studied  the  response  of  a  capillary  jet  to  external 
disturbances  generated  by  a  loudspeaker  and  audio  oscilla¬ 
tor  Donnelly  and  Glaberson,  like  Crane  etal.,  noted  parcels 
of  liquid  between  the  large  drops  predicted  by  Rayleigh’s 
theory.  They  termed  these  parcels  “ligaments”  and  account¬ 
ed  for  their  appearance  by  arguing  that  the  ligaments  were 
due  to  higher-order  harmonics  present  in  the  disturbing  fre¬ 
quency.  Donnelly  and  Glaberson  found  excellent  agreement 
between  Rayleigh’s  linearized  theory  and  their  experimental 
results,  despite  the  fact  that  Rayleigh's  theory  does  not  pre¬ 
dict  the  existence  of  ligaments. 

During  the  past  fifteen  years  we  have  witnessed  a  large 
volume  of  research  aimed  at  explaining  and  predicting  the 
formation  of  ligaments  in  the  process  of  varicose  break-up; 
examples  of  this  research  effort  are  the  theoretical  work  of 
Yuen,7  Nayfeh,®  and  Lafrance.’  Experimentally,  the  liga¬ 
ment  and  satellite  drop  formation  mechanism  was  investi¬ 
gated  by  Goedde  and  Yuen,10  Rutland  and  Jameson"  and, 
in  a  comprehensive  three-paper  study,  by  Chaudhary  and 
Redekopp,'2  and  Chaudhary  and  Maxworthy  13 14 

Relative  to  the  wealth  of  information  on  the  varicose 


regime,  the  nonaxisymmetric  breakup  is  practically  un¬ 
known.  Photographs  of  the  meandering  path  of  fast  capillary 
jets  appeared  as  early  as  1931  in  the  writings  of  Weber15  and 
Haenlein16:  these  photographs  were  reproduced  later  in  a 
famous  textbook  by  Prandtl,17  who  referred  to  the  crests  of 
nonaxisymmetric  shape  as  “wavy  bulges.”  The  subject  of 
nonaxisymmetric  breakup  resurfaced  only  recently  in  the 
literature,  triggered  by  the  need  for  improved  fire-fighting 
equipment.  Hoyt,  Taylor,  and  Runge18  reported  an  experi¬ 
mental  study  of  the  breakup  of  fast  water  jets  and  on  the 
effect  of  adding  drag-reducing  polymer  to  the  water  solu¬ 
tion.  The  authors  refer  to  the  meandering  section  of  the  jet  as 
an  “unstable  wave  region."  Greater  photographic  resolution 
of  the  meandering  breakup  regime  was  achieved  in  a  subse¬ 
quent  descriptive  study  by  Hoyt  and  Taylor. 19 

In  summary,  much  of  the  existing  work  on  the  breakup 
of  capillary  jets  has  dealt  with  the  axisymmetric  (varicose) 
regime.  The  work  on  the  nonaxisymmetric  regime  is  sketchy 
and,  in  all  cases,  qualitative.  The  object  of  this  paper  is  to 
report  a  quantitative  study  of  the  nonaxisymmetric  flow  re¬ 
gime  of  fast  capillary  jets.  For  the  first  time,  the  photo¬ 
graphed  shape  of  such  jets  is  subjected  to  a  rigorous  wave¬ 
length  analysis  which  shows  conclusively  that  the 
nonaxisymmetric  shape  is  governed  by  a  characteristic, 
meander-type,  wavelength  which  scales  with  the  jet  diame¬ 
ter. 


FIG  I  Sc  lematic  of  enpenmental  apparatus  and  nozzle  design. 
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II.  EXPERIMENT 

The  breakup  modes  of  a  capillary  jet  issuing  into  the 
surrounding  atmosphere  were  studied  in  the  laboratory  us¬ 
ing  the  apparatus  shown  in  Fig.  1.  The  fluid  reservoir  con¬ 
sisted  of  a  1.83  m  tall  Plexiglas  cylinder  with  an  internal 
diameter  of  14  cm.  The  reservoir  had  a  number  of  fluid 
drainage  ports  distributed  equidistantly  over  its  height.  The 
nozzle  adaptor  included  a  “rounded”  internal  duct  design, 
and  was  located  1 S  cm  from  the  bottom  of  the  cylinder  in 
order  to  avoid  the  flow  distortion  caused  by  the  bottom.  The 
reservoir  was  safely  pressurized  to  2  atm  (30  psig)  while  the 
cylinder  was  full,  yielding  a  jet  velocity  range  of  0-25  m/sec. 
The  range  of  low  jet  velocities  was  produced  without  pres¬ 
surization  by  controlling  the  height  of  the  reservoir  column 
via  an  appropriate  drainage  port.  Figure  1  shows  also  the 
two  nozzles  employed  in  this  study.  The  3  mm  nozzle  was 
made  from  a  plastic  compound  which  was  cast  in  a  preci¬ 
sion-made  mold  and  later  machined  to  final  dimensions.  The 

1.1  mm  nozzle  was  machined  directly  from  a  Plexiglas  rod. 

The  jet  flow  was  recorded  photographically  using  the 
set  up  shown  schematically  in  Fig.  2.  The  photographic 
equipment  consisted  of  a  Hasselblad  500EL/M  view  camera 
fitted  with  extension  tube  No.  21  for  detailed  closeup  shots,  a 
Sunpak  model  320  photoflash  and  a  45  cm  X  75  cm  section  of 
translucent  glass  for  diffusing  the  light  from  the  flash.  As 
shown  in  Fig.  2,  the  jet  was  positioned  between  light  source 
and  camera.  The  proper  combinations  of  f-stop  and  shutter 
speed  (in  total  darkness  at  “full”  flash  power),  determined 
after  a  number  of  trials,  are  reported  here  in  Table  I. 

To  provide  a  reasonable  range  of  fluid  properties,  this 
study  was  based  on  three  different  fluids; 

I.  distilled  water, 

II.  glycerol  in  water  solution,  30%  by  volume, 

III.  glycerol  in  water  solution,  70%  by  volume. 

The  physical  properties  of  the  three  fluids  are  reported  in 
Table  II.  For  each  of  the  three  fluids  and  the  two  nozzles, 
four  different  jet  velocities  ranging  from  2  m/sec  to  20  m/sec 
were  observed.  Thus,  a  total  of  24  jets  were  observed  and 
recorded. 


FIG.  2.  Photographic  arrangement,  and  the  coordinates  of  the  photo- 
graphed  jet  segment. 


TABLE  I.  Proper  exposure  settings  for  photography  in  total  darkness  at 
“full"  flash  power. 


Film  type 

f-stop 

Shutter  speed 

Kodak  EPR  120 

4 

1/30  sec 

Kodak  PXP  120 

56 

1/30  sec 

III.  THE  CHARACTERISTIC  WAVELENGTH 

The  domain  covered  by  the  present  study  is  shown  on 
the  Weber  number-Reynolds  number  chart  of  Fig.  3.  The 
following  definitions  apply, 

We  =  pV1D  /a,  (1) 

R  e=pVD/p  (2) 

where  p,  V,  D,  a,  and  p  are,  respectively,  the  jet  density, 
velocity,  diameter,  surface  tension  (in  contact  with  air),  and 
viscosity.  As  illustrated  in  Fig.  2,  (a)  and  (b)  denote  the  ex¬ 
tremities  of  the  photographed  portion  of  the  jet.  The  com¬ 
plete  photographic  record  is  available  in  a  thesis  written  by 
Stockman.20  Due  to  space  limitations,  in  the  present  paper 
we  analyze  only  a  representative  sample  of  this  record. 

Three  photographs  of  the  fast  capillary  jet  flow  are 
shown  in  Figs.  4(a),  5(a),  and  6(a).  Each  photograph  corre¬ 
sponds  to  one  of  the  three  different  liquids  used  in  this  study. 
It  is  useful  to  take  a  close  look  at  the  shape  (contour)  of  the 
photographed  jets  in  order  to  recognize  the  large-scale 
meandering  path  followed  by  the  jet.  As  the  jet  fluid  viscosity 
increases  from  Figs.  4(a)  to  Fig.  6(a),  the  sharpness  of  the 
meandering  path  is  enhanced  to  the  point  where,  in  Fig.  6(a), 
the  sinusoidal  contour  of  the  jet  is  illustrated  with  amazing 
clarity. 

The  central  objective  of  our  study  was  to  document  in 
quantitative  terms  the  meandering  shape  of  fast  capillary 
jets.  To  meet  this  objective,  the  jet  contours  were  subjected  to 
a  wavelength  analysis.  In  each  case,  the  jet  contour  was  pro¬ 
jected  (enlarged)  on  a  screen  and  traced  by  hand  on  paper. 
This  operation  produced  two  waveforms,  one  for  the  upper 
edge  of  the  jet  column  and  another  for  the  lower  edge.  The 
waveforms  were  then  digitized  and  fed  into  a  computer  pro¬ 
gram  which  calculated  their  Fourier  transforms  and  deter¬ 
mined  the  respective  power  spectra  and  cross  correlation 
functions. 

In  order  to  learn  how  the  characteristic  wavelength  var- 


TABLE  II.  Physical  properties  of  the  working  fluids. 


Fluid 

Density 

Ig/cm'l 

Viscosity 

IcSl 

Surface 

tension 

lg/s:l 

I. 

Distilled  water 
(18*0 

1 

1 

73 

II 

Glycerol-water, 

30%  by  volume  l20  'Cl 

1.18 

17 

68 

III 

Glycerol-water, 

70%  by  volume  l20  ’Cl 

1  24 

333.6 

64.5 
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FIG.  3.  Weber  number-Reynolds  number  domain  covered  by  the  present  experiments. 


ies  with  position  along  the  jet,  each  photographed  contour 
was  divided  into  a  number  of  segments  [for  example,  three 
segments  for  the  jet  column  shown  in  Fig.  4(a)].  Each  seg¬ 
ment  was  analyzed,  and  the  results  are  presented  in  Figs.  4-6 
as  power  spectra  with  A  /D  on  the  abscissa  and  P  *  on  the 
ordinate.  P  •  is  defined  by 

P*=  -^lim-Lf  |*(/)|2<*,  (3) 

L)  r— »  i  jo 

*(/)-  P  x(t)e-^dt,  (4) 

J  — ■  • 

where  Jt(r )  is  the  contour  waveform,  a>  —  2nff=  V /A,  and  T 
is  the  sample  length. 

Figures  4(b)— 4(d)  demonstrate  that,  regardless  of  longi¬ 
tudinal  position  along  the  jet  axis,  the  upper  and  lower  wave¬ 
forms  have  a  single  (narrow)  band  of  wavelengths  which 
dominate  the  power  spectrum.  In  this  study  we  refer  to  the 
predominant  wavelength  (/l  corresponding  to  maximum  P  •) 
as  the  meander  wavelength  A^ .  Figures  4(b)-4(d)  show  also 
that  the  upper  and  lower  waveforms  have  the  same  meander 
wavelength  and,  in  all  cases,  the  meander  wavelength  scales 
with  the  jet  diameter.  Furthermore,  the  upper  and  lower 
waveforms  are  in  phase:  this  conclusion  follows  from  the 
cross  correlation  function 

1  1  CT 

r*  =  —  lim  —  x{t  |y(r  +  r]dt,  (5) 

D  r-*®  /  Jo 

where  x{t )  and y\t )  are  the  two  waveforms,  and  r  is  the  prede¬ 
termined  phase  shift  (lag)  between  the  two  waveforms.  As 
shown  in  Fig.  4(e)  the  cross  correlation  function  r*  reaches 
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its  peak  value  at  zero  lag,  which  indicates  that  the  upper  and 
lower  waves  are  in  phase.  In  conclusion,  the  flow  regime 
documented  in  this  study  is  not  axisymmetric  (varicose),  but 
one  which  is  characterized  by  a  large-scale  sinuous  shape  of 
wavelength 

Similar  conclusions  regarding  the  existence  of  a  charac¬ 
teristic  meander  wavelength  emerge  from  the  analysis  of 
Figs.  3(a)  and  6(a).  Due  to  space  limitations,  only  two  sam¬ 
ples  are  reported  here  as  Figs.  5(b)  and  6(b),  while  the  com¬ 
plete  record  of  the  wavelength  analysis  may  be  found  in  Ref. 
20. 

Figure  7  shows  a  summary  of  the  measurements 
yielded  by  the  present  study.  The  plotted  Am„  represents  the 
average  over  the  given  x  segment,  however,  in  reality  the 
meander  wavelength  is  continuous  in  x. 

IV.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

A  theoretical  interpretation  of  the  present  results  is  pos¬ 
sible,  based  on  both  the  theory  of  hydrodynamic  stability 
and  the  buckling  theory  of  fluid  columns.  Batchelor  and 
Gill2 1  considered  the  linear  stability  problem  associated  with 
an  inviscid  round  jet  discharging  into  a  quiscent  fluid.  They 
showed  that  sufficiently  far  downstream  from  the  nozzle,  the 
jet  is  least  stable  to  a  temporal  nonaxisymmetric  ("sin¬ 
uous"21)  distrurbance  whose  axial  wavelength  is  larger  than 
several  times  the  jet  diameter.  Similar  conclusions  were 
reached  by  Mattingly  and  Chang,22  who  treated  the  linear 
stability  of  spatial  disturbances  imposed  on  the  same  jet  con¬ 
figuration.  In  addition,  Mattingly  and  Chang22  studied  the 

M  G.  Stockman  and  A.  Be)an  1 508 


4HA- 


v.. 


FIG.  4.  Meandering  jet  of  fluid  l,0«  1.1  mm.  V  =  1 7  m/sec,  a  =  0  mm.  6=165  mm.  |a|  photograph,  (b)  spectral  density  of  segment  51-71  mm  down¬ 
stream  from  nozzle,  (c)  spectral  density  of  segment  71-91  mm  downstream  from  nozzle,  (d|  spectral  density  of  segment  91-111  mm  downstream  from 
nozzle,  let  cross-correlation  function  of  segment  51-71  mm  downstream  from  nozzle.  The  upper  and  lower  waveforms  are  labeled  A  and  O,  respectively. 


natural  instability  of  the  jet  in  the  laboratory  and  reported 
excellent  agreement  between  experimental  measurements 
and  theoretical  stability  predictions.  The  same  problem  and 
conclusions  were  discussed  in  a  most  recent  study  by  Lopez 
and  Kurzweg.23 

The  connection  between  the  nonaxisymmetric  shape 
and  fast  capillary  jets  and  the  least  stable  disturbance  pre- 
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dieted  in  Refs.  21-23  was  recognized  by  Hoyt  and  Taylor.24 
Based  on  photographs  similar  to  the  ones  obtained  in  the 
present  study,  Hoyt  and  Taylor24  were  able  to  identify  a  visi¬ 
ble  axial  wavelength  of  what  is  clearly  a  nonaxisymmetric  jet 
shape.  In  Fig.  10  of  their  study,  Hoyt  and  Taylor24  report 
that  the  visible  axial  wavelength  increases  in  the  down¬ 
stream  direction,  much  in  the  same  manner  as  /im„  of  fluids 
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FIG.  5.  Meandering  jet  of  fluid  II,  D  =  3  mm.  V  =  11.2  m/s,  j  =  146  mm, 
b  =  324  mm.  (a)  photograph,  lb)  spectral  density  of  segment  168-183  mm 
downstream  from  nozzle;  the  upper  and  lower  waveforms  are  labeled  A 
and  O,  respectively. 


I  and  II  considered  in  the  present  study  (Fig.  7).  Inspired  by 
the  existing  theoretical  work  of  BatcheL:  and  Gill21  and 
Mattingly  and  Chang,22  Hoyt  and  Taylor24  interpreted  their 
two-dimensional  photographic  record  as  a  helical  (three-di¬ 
mensional)  instability  with  long  axial  wavelength,  as  predict¬ 
ed  by  stability  theory. 

In  a  more  recent  experimental  report.  Freeman  and 
Tavlarides23  showed  that  when  a  liquid  jet  is  suspended  in  a 


FIG.  6.  Meandering  jet  of  fluid  III.  D  =  3  mm,  V  =  15.2  m/sec,  a  =  140 
mm.  b  =  305  mm  lai  photograph,  ibi  spectral  density  of  segment  267-278 
mm  downstream  from  nozzle,  the  upper  and  lower  waveforms  are  labeled 
A  and  O.  respectively 
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FIG.  7.  The  measured  meandering  or  buckling  wavelength  versus  longitu¬ 
dinal  position  along  the  jet. 


cofluent  stream,  the  jet  develops  a  nonaxisymmetric  shape 
when  the  relative  velocity  between  jet  and  stream  reaches  a 
high  enough  value.  Figures  2  and  3(a)  published  by  Freeman 
and  Tavlarides23  show  a  sinuous  contour  whose  ratio  (axial 
wavelength)/(diameter)  appear  to  be  nearly  identical  to  the 
ratio  visible  in  Fig.  6(a)  of  the  present  study.  Specifically, 
averaging  over  three  complete  wavelengths  visible  in  Fig.  2 
of  Ref.  25,  we  obtain  Am.,/Z)~1.67.  Also,  averaging  over 
four  complete  wavelengths  visible  in  Fig.  3(a)  of  Ref.  25  we 
estimate  Am„/Z)~1.47.  Note  that  these  two  values  of 
kmtx/D,  1.67 and  1.47, fall  right  in  the  middle  of  the  narrow 
band  of  characteristic  wavelengths  revealed  by  the  power 
spectrum  of  Fig.  6(b)  in  the  present  study. 

To  summarize,  classical  hydrodynamic  stability  argu¬ 
ments  predict  correctly  the  instability  of  the  jet  column  to 
nonaxisymmetric  disturbances,  as  well  as  the  scale  of  the 
axial  wavelength  of  such  disturbances.  However,  there  is  one 
additional  result  which  now  has  been  documented  by  three 
independent  experiments  (Refs.  24  and  25  and  the  present 
study)  which  is  not  predicted  by  existing  hydrodynamic  sta¬ 
bility  analyses.  This  additional  result  is  the  tendency  of  the 
nonaxisymmetric  •frave  to  show  a  Am„/2>  value  which  ap¬ 
proaches  1.5  in  a  region  close  enough  to  the  nozzle  where,  as 
discussed  by  Hoyt  and  Taylor,24  the  nonaxisymmetric  dis¬ 
turbance  has  not  had  time  to  be  amplified  due  to  the  form 
drag  interaction  between  the  liquid  jet  and  the  ambient  air. 

Insight  into  the  origins  of  this  additional  feature  is  of¬ 
fered  by  the  buckling  of  fluid  columns. 23-31  It  is  worth  noting 
that  as  a  theoretical  viewpoint  in  fluid  mechanics,  the  buck¬ 
ling  theory  is  much  newer  than  hydrodynamic  stability  the¬ 
ory.  The  novelty  of  the  concept  of  fluid  column  buckling 
may  indeed  be  responsible  for  the  early  interpretation  of 
nonaxisymmetric  disturbances  in  fast  capillary  jets  as  heli¬ 
cal.  In  fact,  the  experimental  record  available  for  this  inter¬ 
pretation  is  exclusively  two-dimensional  and,  as  such,  the 
same  record  can  be  interpreted  as  evidence  of  local  buckling 
in  a  plane  determined  randomly  by  the  presence  of  random 
disturbances  at  the  air  interface. 

This  alternative  interpretation  is  recommended  strong- 
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For  ihe  same  reason,  the  fact  that  for  fluids  I  and  II  AmilI 
increases  with  x  is  a  reflection  of  the  thickening  of  the  air 
stream  entrained  by  the  liquid  jet.  Thus,  sufficiently  far 
downstream  the  liquid  jet  meanders  according  to  the  buck¬ 
led  shape  of  the  surrounding  (thicker)  air  stream. 

The  connection  between  the  nonaxisymmetric  flow  of 
fast  jets  and  fluid  buckling  requires  further  study.  Some  re- 
seachers  have  already  expressed  the  view  that  fluid  buckling 
may  serve  as  origin  for  the  turbulent  motion  of  fluids.30-53 
Along  the  same  lines,  it  is  interesting  to  note  Lopez  and 
Kurzweg’s  statement  that  the  nonaxisymmetric  instability 
of  jet  flow  may  actually  account  for  “the  breakdown  pheno¬ 
menon  in  boundary  layer  flow,"  hence,  for  the  well-docu¬ 
mented  bursting  process.34 


FIG.  8  .The  theoretical  (Ref.  30.3 1)  shape  of  a  buckled  in  viscid  jet  vis-a-vis  a 
closeup  of  Fig.  6(aj. 


ly  by  the  first  observations  of  fast  capillary  jet  flow,  in  parti¬ 
cular  by  Weber’s15  Figs.  1(b),  1(d),  and  1(e).  In  these  early 
photographs  the  lateral  deformation  of  the  liquid  column 
has  a  local  (nonperiodic)  character,  however,  the  wavelength 
of  this  deformation  is  always  a  characteristic  multiple  of  the 
jet  thickness.  The  same  effect  is  visible  in  Fig.  6(a)  of  the 
present  study,  where  entire  sections  of  the  jet  appear  to  be 
undisturbed  while  other  segments  show  the  characteristic 
wave  A^, . 

Theoretically,  it  has  been  shown  that  a  column  of  vis¬ 
cous  fluid  can  buckle  in  a  way  similar  to  rods  in  axial  com¬ 
pression,^ 16,27  however,  the  wavelength  of  the  buckled  shape 
depends  solely  on  the  wavelength  of  the  initial  lateral  distur¬ 
bance.  On  the  other  hand,  for  a  column  of  inviscid  fluid,  it  is 
found  that  the  buckling  wavelength  always  scales  with  the 
column  diameter,30,31  A  /D  =  tt/2  =  1.57.  This  prediction 
can  then  serve  as  explanation  for  the  observed  axial  wave¬ 
length  of  incipient  nonaxisymmetric  deformations  in  fast 
capillary  columns. 

Figure  8  compares  graphically  the  theoretical  buckling 
shape  of  a  round  jet30,3 1  with  a  closeup  view  of  the  jet  photo¬ 
graphed  in  Fig.  6(a).  The  similarity  of  the  two  waveforms  is 
striking. 

A  strong  indication  that  the  buckling  theory  of  inviscid 
jets  accounts  correctly  for  the  meandering  wavelength  docu¬ 
mented  in  this  study  is  that  the  wave  of  Fig.  6(a)  agrees  well 
with  the  wave  photographed  by  Freeman  and  Tavlarides.25 
Note  that  for  fluid  III  the  Reynolds  number  is  in  the  range 
10-100  and  the  Weber  number  in  the  range  103— 104;  these 
ranges  are  quite  different  from  Ref.  25  where  Re  — 103  and 
We  —  14.  The  fact  that  despite  such  discrepancies  the  photo¬ 
graphed  wavelengths  agree  with  each  other  and  with  the 
buckling  wavelength  rrD  /2  supports  the  buckling  theory 
very  strongly.  The  key  result  of  the  buckling  theory  is  that 
the  buckling  wavelength  must  depend  only  on  D,  in  other 
words,  it  must  be  independent  of  V and  physical  properties. 
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Theory  of  Instantaneous  Sinuous  Structure  in  Turbulent  Buoyant  Plumes 


A.  Bejan,  Boulder.  Colorado 


Abstract.  This  paper  reports  a  theory  which  explains  the  flicker- 
ins  motion  of  turbulent  plumes  as  well  as  their  large-scale  sinuous 
structure.  The  theory  is  based  on  the  fact  that  the  inviscid  region 
of  the  plume  (the  "plume  column")  possesses  elastic  properties 
analogous  to  those  of  elastic  rods  subjected  to  longitudinal 
compression.  It  is  shown  that  the  straight  plume  column  is  not 
stable  and  buckles.  The  distance  between  two  consecutive  elbows 
is  proportional  to  the  local  plume  diameter,  in  other  words,  the 
shapes  of  all  buckled  piumes  are  geometrically  similar.  A  buckled 
plume  collapses  periodically  due  to  the  interaction  of  its  lateral 
elbows  with  the  stagnant  ambient.  This  interaction  is  responsible 
for  the  intermittent  formation  of  large-scale  buoyant  eddies  on 
the  periphery  of  the  turbulent  plume. 

Theorie  der  veranderlichen  getvundenen  Struktur  in  turbulenten 
Auftriebsfahnen 

Zusammenfassung.  Diese  Arbeit  versucht,  die  flattemde  Bewe- 
gung  turbulenter  Auftriebsfahnen  und  ihre  groBriiumige  gcwun- 
dene  Struktur  aufzuklaren.  Man  geht  davon  aus,  daB  der  nicht- 
ziihe  Bereich  der  Fahne  (die  ..Fahnensaule")  elastische  Eigen- 
schaften  besitzt  ilhnlich  jener  elastischer  Stabe.  die  longitudinalcn 
Verdichtungen  unterworfen  sind.  Die  gerade  Fahnensaule  ist 
nicht  stabil  und  verbiegt  sich.  Der  Abstand  zwischen  zwei  Krtim- 
mungen  ist  dem  drtlichen  Fahnendurchmesser  proportional,  d.  h. 
die  Formen  aller  eekrummten  Fahnen  sind  geomeinsch  ahnlich. 
Die  gekriimmte  Fahne  bricht  periodiseh  zusammen  infoige  der 
Wechselwirkung  ihrer  zeillichen  Kriimmungen  und  der  ruhenden 
L  mgebung.  Dadurch  entstehen  intermittierend  groBriiumige  Auf- 
triebswirbel  am  L'mfang  turbulenter  Auftriebsfahnen. 


1  Introduction 

Turbulent  plumes  represent  an  extremely  frequent  natural 
phenomenon,  in  the  atmosphere  as  well  as  in  the  hydro¬ 
sphere.  Furthermore,  in  heat  transfer  and  environmental 
engineering,  the  turbulent  plume  constitutes  one  of  the 
most  effective  mixing  mechanisms  known  to  man.  The 
state-  of  knowledge  on  turbulent  plumes  has  been  sum¬ 
marized  by  Turner  [1|  and.  as  part  of  a  comprehensive 
monograph  on  buoyancy-driven  (lows,  by  Turner  [2], 
Another  extensive  overview  of  the  field  was  published 
most  recently  by  Fischer.  List.  Koh.  Imberger  and  Brooks 
[  >].  A  critical  examination  of  these  review  works  leads  to 


the  conclusion  that  the  favored  approach  in  turbulent 
plume  research  is  one  where  the  "complications"  of  the 
turbulent  flow  field  can  be  smoothed  out  by  the  classic 
method  of  time-averaging  (Townsend.  [4]).  Consequently, 
most  of  the  theoretical  work  on  buoyant  jets  relies  on  the 
a.xtsy mmteric  flow  model  which  lends  itself  very  well  to 
boundary  layer-type  analysis. 

Although  mathematically  attractive,  the  time-averaged 
turbulent  plume  concept  does  gross  injustice  to  the  phys¬ 
ics  of  the  phenomenon.  A  turbulent  exhaust  plume  does 
not  rise  straight  up  into  the  air,  with  a  fixed  shape 
resembling  that  of  an  inverted  cone  or  funnel.  A  real 
plume  executes  a  periodic  lateral  movement:  the  plume  is 
not  straight,  rather,  it  has  recognizable  bends  separated  by 
distances  of  the  same  order  as  the  plume  diameter.  The 
smoke  rising  from  a  barbecue  or  camp  fire  "snakes"  into 
the  air.  This  large  scale  motion,  which  ts  so  effectively 
blacked  out  by  time-averaging,  can  be  seen  clearly  in 
Fig.  I. 

The  object  of  this  paper  is  to  present  a  theory  which 
explains  for  the  first  time  the  instantaneous  structure  and 
the  "flickering"  nature  of  turbulent  piumes. 


2  The  Concept  of  Elastic  Plume  Column 

Consider  a  fluid  of  density  o  rising  vertically  through  a 
heavier  fluid  of  density  n  +  Jo.  We  model  the  starting  sec¬ 
tion  of  this  buoyant  stream  as  inviscid.  This  means  that  for 
a  certain  height  at  the  base  of  the  plume  we  are  neglecting 
the  shear  interaction  between  the  plume  fluid  and  the 
stagnant  ambient.  This  is  permissible  for  a  height  of  the 
order  of  4-6  times  the  plume  diameter,  where,  as  shown 
by  Crow  and  Champagne  [5],  the  plume-ambient  shear 
layer  is  thin  compared  with  the  plume  diameter  Having 
neglected  the  interaction  at  the  plume-ambient  interface, 
it  is  helpful  to  imagine  the  plume  fluid  as  being  sur¬ 
rounded  by  a  thin  flexible  sleeve  (control  surface)  wlt-ch 
is  stationary  with  respect  to  the  stagnant  ambient 
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Fig.  I.  Large-scale  periodicstrueturein  the  plume 
above  a  natural  gas  well  on  fire  -  Douglas  Pass. 
Colorado  (AP  Laserphoto.  reprinted  from  the 
Longmont  Daily  Times-Call. Oct.  lSs-19.  1 9X0) 


Imagine  now  a  cut  perpendicular  to  the  plume  axis.  If 
the  column  is  straight,  in  the  cross-section  generated  bv 
the  cut  we  distinguish  oniv  the  compressive  force  (impulse 

[6]i 

C  =  f f  o  1  d  4  .  (2.1) 

t 

where  L  and  A  are  the  plume  velocity  and  cross-sectional 
area.  If  the  plume  column  is  slightly  curved,  we  distin¬ 
guish  also  a  cross-  actional  bending  moment  M.  This 
bending  moment  is  due  primarily  to  the  fact  that  the 
plume  fluid  travels  relatively  faster  through  regions  of  the 
cut  which  are  located  on  the  inside  bank  of  the  bend.  It  is 
easy  to  show  that  the  cross-sectional  bending  moment  M 
is  present  as  soon  as  the  column  acquires  curvature. 
Consider  for  this  purpose  a  bending  test  in  which  the 
uniform  stream  (o.  V.  .4)  passes  through  a  slightly  curved 
duct:  the  radius  of  curvature  of  the  duct  centerline  Rr_.  is 
infinitely  greater  than  the  plume  diameter.  D.  The  Bernoulli 
equation  applied  along  a  streamline  requires  i  Prandtl  [7j). 


where  the  left  hand  side  applies  to  the  straight  section 
(/’  =  ())  while  the  right  hand  side  corresponds  to  the  curved 
section.  Equilibrium  in  the  radial  direction  requires  also 
(Prandtl  [8]). 


cjtR_  dP 
R-c  dr 


(2.3) 


where  i  is  the  radial  coordinate  measured  in  the  plane  of 
curvature  away  from  the  plume  centerline,  toward  the 
center  of  curvature  Combining  I  qs.  (2.2)  and  (2.3).  we 
obtain 


I  I  i  - 


R . 


•  o 


(2.4) 


The  cross-sectional  bending  moment  M  can  now  be  calcu¬ 
lated  by  writing 

3/  =  j)  (o C  +  P)  v  d.-t  (2.5) 

a 

which,  using  Eqs.  (2.4)  and  (2.2).  yields 
ol-  ... 

3/  =  ^—  ffrd.4.  (2.6) 

''X  .1 

Note  that  the  area  integral  appearing  in  Eq.  (2.6)  is  the 
area  moment  of  inertia  (/  =  nD't 64)  employed  routinely 
in  the  study  of  elastic  beam  flexure  (Den  Hartog  [9]).  The 
group  o  V"  plays  the  role  of  “modulus  of  elasticity":  this  is 
why  it  becomes  increasingly  difficult  to  manually  bend  a 
hose  as  the  turbulent  (inviscid)  flow  rate  through  it 
increases. 

At  this  point  we  conclude  that  the  plume  column  is  in 
axial  compression  and  that  the  cross-sectional  bending 
moment  is  proportional  to  the  curvature  of  the  centerline 
(I  /?*).  With  such  properties,  the  straight  plume  column 
becomes  analogous  to  an  elastic  column  subjected  to  axial 
compression  (Den  Hartog  [10]).  This  analogy  guarantees 
that  if  the  plume  is  tall  enough,  it  will  buckle  like  an 
elastic  rod  compressed  between  the  ends.  In  the  next 
section  we  develop  the  ’'buckling  theory"  of  turbulent 
plumes  viewed  as  elastic  columns. 

3  The  Natural  Buckling  of  a  Straight  Plume 
Rising  from  Rest 

Consider  the  stationary  control  volume  defined  bv  two  cut' 
normal  to  the  plume  column,  at  r  and  r  dr  d  ig  2).  The 
balance  of  vertical  forces  acting  on  this  volume  dictates 
df 

—  =  U>  4  q  .  <  5. 1 ) 

dr 
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Fig.  2.  Static  equilibrium  of  a  segment  of  vertical  plume  column 


The  conservation  of  mass  through  the  control  volume 
requires  also 

V(z)  A  (:)  =  Q,  constant,  (3.2) 

where  Q  is  the  volumetric  flowrate  through  the  plume 
column.  Combining  these  conservation  statements  we 
obtain 


r<--)  =  (2^0--)  ,  (3.3) 

where  we  made  the  additional  assumption  that  the  plume 
fluid  originates  from  rest.  K(0)  =  0. 

As  in  the  buckling  theory  of  a  vertical  flagpole  (Timo¬ 
shenko  and  Gere  [11];  Den  Hartog  [12]),  the  rotational 
equilibrium  of  the  plume  element  dc  requires 

C  dY  +  \f  -  (M  +  dM)  =  0  (3.4) 

or,  using  Eq.  (2.1)  and  (2.6). 


n  A  V* 


dY_ 

dr 


_d_ 

dr 


Hi-/ 


=  0. 


(3.5) 


Note  that  in  the  limit  of  vanishingly  small  curvature. 
\/Rx  -  -  d-)'/dr:.  Taking  into  account  the  r-dependence 
of  V  and  A.  Eq.  (3.3),  the  equation  of  flexure  (3.5) 
becomes 

d3K  4 7r(2 g  Ao/o)vi  d  Y 

17+ - (36) 

The  first  important  conclusion  of  this  analysis  is  that  the 
plume  column  develops  a  “structure”  whose  characteristic 
vertical  dimension  scales  with 


Q; 


32  ,t-  cf  J  o/o 


The  vertical  length  scale  can  also  be  written  as 


constant 


(3.7) 


(3.7a) 


where  D( :)  is  the  plume  diameter  at  level  r  (note  that 
D( :)  —  r-1'1).  Introducing  the  dimensionless  coordinate 
v  =  r/r0.  the  equation  of  flexure  becomes 


d1)' 

dr:l 


dT 

d.v 


0 . 


<3  s> 


The  general  solution  to  this  equation  is  expressible  in 
terms  of  Bessel  functions  (Watson  [13]) 

=  A\  xh2J2/ 3(4/5  .v5'*)  +  K,  xv2J-2! 5(4/5  .vs,i)  (3.9) 
d.v 

where  A,  and  K..  are  arbitrary  constants.  Noting  that 
d)' 

hm  —  =  A,  .v  +  K, .  (3.10) 

r-0  d.v 

we  set  A,  =  0.  which  accounts  for  the  fact  that  in  the  very 
beginning  the  buoyant  fluid  rises  vertically.  The  shape  of 
the  plume  column.  K(.v).  is  obtained  by  integrating  (3.9). 

>’(.v)  =  A  f  J zo(m)  dm  (3.11) 

o 

vv  here 

2  =  4/5  .V3'4.  (3.12) 


The  integral  appearing  in  (3.1 1)  was  calculated  with  sixth¬ 
digit  accuracy,  by  first  expanding  J2,i (m)  in  a  power 
series.  The  principal  values  of  this  integral  are  reported 
here  in  Table  I.  because  they  are  not  available  in  the 
literature. 


Table  I 


J’  my*  dm 

0 

0.5 

0.172638 

1 

0.492870 

1.5 

0.852135 

*) 

1.168617 

2.5 

1.380308 

3 

1.452918 

3.5 

1.384381 

4 

1.203166 

4.5 

0.960384 

5 

0.717376 

5.5 

0.531599 

6 

0.443935 

6.5 

0.470241 

7 

0.598886 

7.5 

0.794658 

8 

1.007952 

8.5 

1.186949 

9 

1.289889 

9.5 

1.294589 

10 

1.203059 

10.5 

1.040267 

II 

0.847522 

11.5 

0.672151 

12 

0.556040 

12.5 

0.525758 

13 

0.586530 

MO 
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Fig.  3.  The  buckled  shape  of  the  plume  column  centerline 


The  magnified  shape  of  the  buckled  plume  column  is 
shown  in  Fig.  3  as  the  function  >'(.v)/A'.  It  should  be 
remembered  that  Eq.  (3.11)  is  the  result  of  a  small 
amplitude  analysis  (  df/dr  <1  1).  hence,  the  amplitude  K 
is  negligible  when  compared  with  the  vertical  length 
scale  The  static  equilibrium  of  the  straight  plume 
column  is  indifferent  [10];  however,  given  the  slightest 
lateral  disturbance,  the  plume  buckles  in  accordance  with 
the  vertical  periodicity  shown  in  Fig.  3.  The  buckled 
shape  is  a  succession  of  elbows  whose  spatial  frequency 
increases  with  height.  The  first  four  elbows  are  located  at 

v.  *  2.878.  a  .  =  5.101,  v,  =  7.101,  v,  =  8.967. 

According  to  the  classical  arguments  of  Hydrodynamic 
Stability  Theory  [14].  the  sinuous  shape  will  be  amplified 
due  to  the  dynamic  interaction  of  the  elbows  with  the 
stagnant  ambient.  It  is  this  interaction  which  leads  to  the 
formation  of  large  turbulent  eddies  visible  in  most  atmo¬ 
spheric  plumes  (Fig.  I). 


4  Formation  and  Evolution  of  Large  F'.ddies 

T  he  mechanism  of  elbow  turbulence  formation  is  shown 
schematically  in  Fig.  4  Speaking  only  qualitatively,  in  the 
first  phase  of  the  process  the  segment  of  plume  column 
located  in  the  vicinity  of  a  natural  elbow  is  pushed  to  the 
side  hy  the  horizontal  resultant  of  the  two  axial  compres¬ 
sive  forces.  The  lateral  movement  of  the  elbow  is  eventu- 
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ally  blocked  by  the  pressure  build-up  associated  with 
either  the  stagnation  of  some  plume  fluid  into  the  still 
ambient,  or  with  accelerating  some  of  the  ambient  fluid. 
This  pressure  build-up  is  consistent  with  the  fact  that  in  a 
column  of  finite  curvature  the  pressure  increases  in  the 
radial  direction.  Eq.  (2.3).  At  the  end  of  the  first  phase  of 
the  process,  the  elbow  has  degenerated  into  a  buoyant 
eddy  of  mixed  fluid:  this  eddy  rotates  as  shown  in  the 
figure,  and  rises  slower  that  the  unmixed  plume  fluid 
found  at  the  same  altitude. 

The  second  phase  of  the  mechanism  is  triggered  by  the 
effect  of  high  pressure  nodules  already  present  on  both 
sides  of  the  plume  column:  the  plume  buckles  in  a  new 
mode,  feeding  its  stream  through  the  structures  generated 
by  the  old  set  of  elbows.  At  the  end  of  the  second  phase, 
the  new  column  gives  birth  to  a  new  set  of  large  eddies 
which  continue  to  rise.  In  general,  the  eddies  produced  by 
the  upper  elbows  interfere  with  older  eddies  rising  from 
lower  levels.  This  interaction  is  partly  responsible  for 
lateral  growth  of  the  turbulent  plume,  as  v  increases. 

The  most  regular  feature  in  the  evolution  of  the 
turbulent  plume  column  is  the  root  section.  0  <  v  <  a,  . 
which  precedes  the  formation  of  the  first  elbow.  The  root 
of  the  plume  swings  back  and  forth,  not  necessarily  in  the 
same  plane,  with  a  characteristic  frequency.  /,  =  1/(2  Jr). 
The  half  period  At  of  this  flickering  motion  can  be 
calculated  by  integrating  the  first  of  equations  (3.3):  At  is 
the  time  interval  needed  by  plume  fluid  to  rise  from  :  =  0 
to  the  first  elbow.  r,=.v,r0.  Combining  this  result  with 
Eqs.  (3.7)  and  (3.3).  we  discover  that  the  Strouhal  number 
based  on  the  plume  velocity  and  diameter  at  the  first 
elbow',  r,  and  Z>, .  is  a  universal  constant  for  all  inviscid 
(turbulent)  plumes 

St,  =  ^i  =  0.534.  14.1) 

1 1 

Another  interesting  result  of  the  buckling  theory  is  that 
the  large  elbow  structure  (the  distance  between  elbows) 
scales  only  with  the  plume  diameter,  which  means  that 
the  buckled  shapes  of  all  inviscid  plumes  are  geometrically 
similar.  This  prediction  agrees  with  the  observed  structure 


v  v 


Fig.  4.  Mechanism  of  formation  of  large-scale  eddies  in  a  turbu¬ 
lent  plume  rising  from  rest 
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Fig.  5.  Universal  geometry'  of  buckled  plume  columns 


of  atmospheric  plumes.  According  to  Eq.  (3.7a).  the 
plume  diameters  corresponding  to  the  first  four  elbows 
are. 

—  =  ~4r  =  3.071,  2.662.  2.45.  2.332.  (4.2) 

-0  ^ » 

Figure  5  shows  a  scale  drawing  of  a  round  plume  and  the 
relative  location  of  the  first  four  elbows.  This  universal 
geometry  agrees  very  well  with  visual  observations  of 
turbulent  plumes  rising  above  barbecues  and  camp  fires. 
It  also  agrees  with  the  plume  photographed  in  Fig.  1. 

5  l  aminar  Plumes  vs.  Turbulent  Plumes:  The  Transition 

The  instantaneous  sinuous  structure  and  the  eddy  forma¬ 
tion  mechanism  described  so  far,  apply  only  to  large 
Reynolds  number  (turbulent)  plumes  which  can  be  mod¬ 
eled  as  inviscid.  The  present  theory  does  not  apply  to 
laminar  plumes.  It  is  important  to  note,  however,  that  the 
universality  of  the  sinuous  structure  (Fig.  5)  povides  a 
theoretical  explanation  for  the  phenomenon  of  transition 
from  laminar  to  turbulent  plume  flow. 

According  to  F.q.  (4.1).  the  base  of  the  plume  fluctua¬ 
tes  within  a  characteristic  time  interval  which  is  propor¬ 
tional  to  the  diameter  and  inversely  proportional  to  the 
velocity.  This  local  characteristic  time. 


:ai 

represents  the  "natural  heartbeat"  of  the  plume  as  an 
inviscid  stream.  If  the  plume  flow  is  highly  viscous 
(laminar),  the  plume  column  has  a  different  characteristic 
time:  the  viscous  local  time  is  associated  with  the  travel  of 
information  by  viscous  diffusion  from  the  plume-ambient 
interface  to  the  plume  centerline  (total  distance  D'2\.  The 
viscous  communication  time  is  given  by  the  solution  to 
Stokes'  first  problem  [15], 


—  (  1'  t Visions) 


1.  or  ( visioi 


The  plume  flow  remains  laminar  if  the  ambient  can 
communicate  viscously  with  the  stream  faster  than  the 
stream  can  fluctuate.  If  the  column  can  fluctuate  faster 
than  the  viscous  diffusion  time.  then,  of  course,  the  flow 
opts  for  the  inviscid  (turbulent)  regime.  At  transition,  the 
viscid  and  inviscid  time  scales  are  of  the  same  order  of 
magnitude:  combining  Eqs.  (5.1)  and  r 5.2 1  we  conclude 
that  the  laminar-turbulent  transition  is  marked  by 


This  theoretical  estimate  of  the  transition  Reynolds  num¬ 
ber  agrees  very  well  with  experimental  observations  and 
with  estimates  based  on  hydrody  namic  stability  experiments. 
Experimentally,  it  is  well  known  that  the  transition  in  free 
jet  and  wake  flow  occurs  in  the  vicinity  of  Re  ~  30  (see. 
for  example,  Schlichting  [16]).  The  hydrodynamic  stability 
of  buoyant  plumes  and  wall  layers  has  been  studied  exten¬ 
sively.  as  summarized  by  Gebhart  [17],  For  e.\umplo.  the 
lowest  Reynolds  number  where  instability  has  been  de¬ 
tected  in  vertical  laminar  boundary  layer  flow  over  a 
constant-flux  wall  is  of  order  67  (Fig.  8-21.  [17]).  The 
difference  between  this  estimate  and  Eq.  (5.3)  is  explained 
by  the  fact  that  in  [17]  the  local  Reynolds  number  (<J*)  is 
based  on  the  thickness  of  the  laminar  boundary  layer.  If. 
instead,  the  Reynolds  number  is  expressed  based  on  the 
displacement  or  momentum  thickness  (more  consistent 
with  the  slug  How  model  employed  in  the  present  theory  I 
the  the  transition  Re  is  a  number  of  the  same  order  as  in 
Eq.  (5.3). 


6  Concluding  Remarks 

The  demonstrated  ability  to  correctly  predict  the  transi¬ 
tion  to  turbulence  is  additional  evidence  that  the  sinuous 
structure  (Fig.  5)  and  the  characteristic  time  scale  (Eq.  4.1 ) 
are  real  properties  of  all  turbulent  buoyant  plumes.  It  is 
important  to  recognize  also  that  (lie  /</n.v-v<i/c  hmklniv 
pht-iumwiwn  described  in  this  paper  is  not  observed  only 
in  buoyant  plumes  as  in  Fig.  I.  The  same  phenomenon  has 
been  observed  and  described  as  "orderly  structure",  "cork¬ 
screw  shape",  and  “whiplash  motion"  in  the  starting 
section  of  low  speed  turbulent  jets  (Crow  and  Champagne 
[5],  Reynolds  [18]).  Others  have  recognized  this  snake-like 
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sliape  as  the  cause  of  the  intermittency  phenomenon  (Yih 
[19].  Fig.  Sb.  and  pp.  545.  546).  In  the  field  of  boiling  heat 
transfer  we  are  familar  with  the  occurrence  of  tall  S- 
sliaped  vapor  bubbles  (buckled  vapor  plumes)  on  intense¬ 
ly  heated  horizontal  surfaces:  this  phenomenon  is  recog¬ 
nized  as  the  "continuous  vapor  column  regime"  (Figs.  2d 
and  2j  in  Moissis  and  Berenson  [20]). 

The  reader  can  easily  reproduce  the  natural  buckled 
shape  by  experimenting  with  the  continuous  water  column 
falling  from  the  kitchen  faucet.  This  water  column  is  a 
"sinking"  plume  defined  not  by  a  flexible  plume-ambient 
interface  as  in  Fig.  I.  but  by  the  flexible  hose  provided  by 
capillary  forces.  Placing  his  finger  about  1-2  cm  under 
the  faucet,  the  reader  can  buckle  this  sinking  plume  into  a 
shape  which  resembles  very  closely  the  shape  shown  in 
Fig.  5.  rotated  by  180°.  In  the  kitchen  faucet  experiment 
the  buckled  column  does  not  break  up.  since  the  lateral 
growth  of  its  elbows  is  suppressed  by  the  effect  of  surface 
tension.  The  same  stabilizing  effect  is  present  in  the  S- 
shaped  tall  bubbles  photographed  during  intense  boiling 
by  Moissis  and  Berenson  [20]. 

From  the  point  of  view  of  theoretical  research  in  fluid 
mechanics,  the  Buckling  Theory  reported  in  this  paper 
represents  a  dramatic  departure  from  accepted  methods 
such  as  Hydrodynamic  Stability  analyses.  For  the  first 
time,  the  Buckling  Theory  focuses  on  the  equilibrium  of 
fluid  flow  as  a  finite  size  system:  the  new  result  is  that  the 
inviscid  plume  possesses  a  new  property  -  buckling  -  and 
that  the  buckled  shape  is  behind  the  observed  "large- 
scale  structure"  of  turbulent  plumes.  The  Buckling  Theory 
predicts  for  the  first  time  the  transition  to  turbulent 
plume  How:  the  transition  is  an  internal  property  of  the 
stream,  the  competition  between  two  stream  properties 
Uriui  tuation  and  /viscous)  Thus,  the  Buckling  Theory  explains 
the  origin  of  turbulent  fluid  motion  [21]. 
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Theoretical  considerations  of  transition 
to  turbulence  in  natural  convection 
near  a  vertical  wall 

A.  Bejan  and  G.  R.  Cunnington* 


Hydrodynamic  stability  analysis  of  an  mviscid  wall  jet  shows  that  instability  is 
possible  above  a  characteristic  disturbance  wavelength  which  is  proportional  to  the 
jet  thickness  This  scaling  is  the  basis  for  an  argument  that  transition  occurs  when 
the  fluctuating  time  period  of  the  unstable  (mviscid)  wall  jet  is  of  the  same  order 
as  the  viscous  diffusion  time  normal  to  the  jet.  The  transition  must  occur  when  the 
jet  Reynolds  number  is  of  the  order  of  102  Published  observations  of  transition 
along  a  heated  vertical  wall  are  reviewed  in  order  to  test  the  validity  of  the  proposed 
scaling  argument  Specifically,  numerous  observations  on  buoyant  jets  near  isother¬ 
mal  walls,  near  constant-heat-flux  walls,  and  in  enclosures  with  vertical  isothermal 
walls  are  shown  to  support  the  validity  of  the  transition  mechanism  proposed 

Key  words:  convection,  turbulence,  fluid  flow 


Transition  phenomena  have  been  studied  extensively 
during  the  past  half-century.  This  work  has  been 
brought  in  perspective  in  a  number  of  review  papers, 
for  example  by  Tani1  and  Reshotko2  for  boundary 
layer  flow  and  by  Gebhart3  for  natural  convection 
flow. 

The  object  of  this  study  is  to  propose  a  scaling 
argument  as  basis  for  transition  to  turbulence  in  the 
wall  jet  flows  encountered  in  natural  convection  along 
heated  vertical  walls.  First,  linear  stability  analysis  of 
a  wall  jet  indicates  that  the  flow  is  unstable  to  disturb¬ 
ance  wavelengths  greater  than  a  certain  multiple  of 
the  jet  thickness  D.  Based  on  this  proportionality,  it 
will  be  argued  that  transition  is  ruled  by  the  internal 
competition  between  two  time  scales,  the  jet  fluctua¬ 
tion  period  and  the  time  of  viscous  penetration  normal 
to  the  flow  direction.  A  comprehensive  review  of  the 
published  experimental  observations  on  transition 
supports  the  validity  of  this  scaling  argument.  Finally, 
it  is  shown  that  photographs  and  numerical  simula¬ 
tions  of  transition  in  vertical  enclosures  heated  from 
the  side  visualise  the  transition  mechanism  described. 

Hydrodynamic  instability  of  an 
inviscid  wall  jet 

Consider  the  flow  of  an  inviscid  jet  U(y)  next  to  a 
vertical  wall,  as  shown  in  Fig  1.  We  want  to  know 
the  waviness  of  the  jet  as  it  becomes  unstable,  in  other 
words,  the  frequency  of  its  fluctuation  relative  to  a 
fixed  spot  on  the  wall,  which  can  be  obtained  by 
performing  a  linearised  stability  analysis  of  the  flow. 
This  analytical  approach  is  outlined  by  Lamb4.  By 
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superimposing  on  the  base  flow  ( u=U ,  v=0) 
unspecified  disturbances: 

u=U  +  u’  (1) 

v  =  v  (2) 


the  inviscid  flow  (Euler)  equations  yield  the  vorticity 
transport  equation: 


dr  ,  Br  ,Bt 

—  +  (U  +  u')—+v—  =  0 
Bt  Bx  By 


(3) 


u  =  u0 
u  =  0 


Fig  I  Wall  jet  velocity  profile 


U(y) 
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J 


with  the  vorticitv  defined  us: 
ilv'  rtu'  d  V 

T  =  - - - - —  (4) 

r).t  fltj  d  IJ 

Linearising  Eq  (3)  yields: 

(it  ..  i>\ (iiv'  itu'\  d- U  , 

(3) 

Next  assume  that  the  disturbance  velocity  components 
are  both  periodic  in  x: 


where  k  is  the  wavenumber  (k  =  2tt/\).  Substitution 
of  Eqs  (6,7)  into  the  vorticitv  transport  equation  and 
the  mass  continuity  equation,  and  the  elimination  of 
u :|[  yields  a  single  equation  for  the  transversal  disturb¬ 
ance  amplitude  e*. 


(<r+kui^-kic*)-j?kv*=° 


Our  interest  is  in  the  wavelengths  A  or 
wavenumbers  k  for  which  cr  is  imaginary,  ie  for  which 
the  assumed  disturbance  is  likely  to  be  amplified  (Eqs 
6,7).  As  shown  in  Fig  1,  we  assume  a  triangular 
velocity  profile,  so  that  in  all  regions  of  the  flow 
d2  U  /  d  i/”  =  0: 


Solving  Eii  (8),  the  corresponding  expressions  for  r 
in  the  three  regions  are: 

I)  c *  =  Ae  + 


III)  u,=  Ee  k"  +  Fek" 

The  solid  wall  condition,  f,  =0  at  (/  =  (),  and  the 
condition  that  rf  must  be  finite  infinitely  far  from  the 
wall,  mean  that: 

A  +  B  =  0  (111 

F  =  0  (12) 

Two  more  equations  follow  the  condition  that  v  <  must 
varv  continuouslv  from  region  I  to  II,  and  from  II  to 
III: 

Ae-kDr-+BekD'2=Ce-ku 2  +  D„eko'z  (13) 

Ce  kD  +  D„  ekn  =  £ e~kn  (14) 

Finally,  the  condition  that  the  pressure  must  vary 
continuously  from  one  region  to  the  next  amounts  to 
integrating  Eq  (8)  across  one  region-to-region  inter¬ 
face  (eg  from  y  =  D ~  to  (/  =  D*)4: 


_  uo 
-  D/2  y’ 


0<  y<  D/2 


II)  V=-gj-(D-y),  D/2<  y<  D 

III)  t/  =  0,  D<  y 


Notation 

A,  B,  C, 
D0,  E,  F 

Coefficients  (Eq  10) 

u*e* 

or,. 

a 

Disturbance  amplitudes 

Critical  wavenumber  (2nL/\,.) 
Thermal  diffusivity 

D 

Jet  thickness 

R 

Coefficient  of  thermal  expansion 

g 

Gravitational  acceleration 

y 

e~kD'-  (Eq  18) 

(in  negative  x  direction) 

y 

Dimensionless  stratification  parameter 

C* 

5  (gjP)  49) 

Thermal  boundary  layer  thickness 

H 

k 

Vertical  length  scale 

Wavenumber  (2it/\) 

5t 

L 

Horizontal  dimension  of  enclosure 

*v 

Boundary  layer  where  vertical 

m 

o-D/uo+1  (Eq  18) 

velocity  obeys  a  no-slip  condition 

N* 

Utmin  (Eq  29) 

A 

Wavelength 

Pr 

Prandtl  number 

r 

Vorticitv 

Ra 

Rayleigh  number:  RaH — based  on 

<7 

Disturbance  growth  rate 

vertical  length  scale  H ;  Rau — based 

V 

Kinematic  viscosity 

Re 

S 

on  cavity  width  L 

Reynolds  number 
vertical  temperature  gradient 

Subscripts  and  superscripts 

t 

T 

Time 

Temperature 

* 

Uniform  heat  flux 

u 

Velocity  parallel  to  the  wall 

• 

Disturbances 

U -direction) 

min 

Minimum 

t 

Velocity  normal  to  the  wall 

V 

Viscous 

( (/-direction) 

B 

Buckling  h  M  ;1- 
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Applying  this  condition  twice,  at  y=D/2  and  at 
tj=D,  yields  two  more  equations: 

A^+ljf^  +  Ce  iH  J  =  0  (16) 

[ 1  ■  £?* + /cu<*)] c + [ 1 + 5^(<r + ku»}]D " 

+  -p~(rr  +  ku„)(A-  B)  =  0  (17) 

4u„ 

Eqs  (11,  13,  14,  16,  17)  are  all  homogeneous  and, 
together,  they  constitute  a  system  for  determining  the 
five  unknown  coefficients  A,  B,  C,  D„,  E.  A  non-trivial 
solution  is  possible  if  the  determinant  of  this  system 
is  zero: 

m~+  m(2y2  +  kD  -  3  -  y 4)  -  y'(  1  +  kD) +  2y2  =  0 

(18) 

where  y  =  e~ka  l  and  m=crD/u0+l,  Eq  (18)  is 
obtained  by  systematically  eliminating  the  five 
coefficients  among  Eqs  (11,  13,  14,  16,  17).  Whether 
or  not  cr  is  imaginary  depends  on  the  character  of  m: 
since  Eq  (18)  is  a  quadratic  in  m,  afflJ  +  bm+c, 
imaginary  roots  are  possible  if  the  discriminant  A  = 
b2-4ac  is  negative,  ie  when: 

A  =  (2  y2  +  k£>  -  3  -  y4)2 + 4  y4(  1  +  kD)  -  8  y2  <  0  (19) 

Solving  A  =  0  we  find  that  the  wall  jet  of  Fig  1  is 
unstable  (A<0)  if 

1.337  <kD<  3.427  (20) 

In  other  words,  the  inviscid  wall  jet  is  likely  to  acquire 
a  waviness  described  by  wavelengths  in  the  range 

1.833  <-^<4.701  (21) 

The  conclusion  that  an  inviscid  wall  jet  of  thickness 
D  is  unstable  only  in  a  certain  range  of  disturbance 
wavelengths,  agrees  with  earlier  results  concerning 
other  inviscid  flows.  For  example,  Rayleigh5  found 
that  a  free  two-dimensional  jet  of  triangular  profile  is 
unstable  for  wavelengths  A  >  1.714D.  Rayleigh  also 
found5  that  a  free  shear  layer  of  thickness  D  is  unstable 
if  A  >  4.914D.  It  is  significant  that  for  several  different 
base  flows,  from  the  wall  jet  of  Fig  1  to  Rayleigh's 
free  shear  layer,  the  edges  of  the  wavelength  domains 
for  inviscid  instability  are  marked  by  wavelengths 
which  scale  with  the  flow  thickness  D.  We  feel  that 
this  scaling  is  important  and,  ultimately,  responsible 
for  the  phenomenon  of  transition  to  turbulence6.  This 
scaling  is  the  basis  for  the  transition  criterion  outlined 
next. 

Time  scale  criterion  for  transition 

Each  longitudinal  length  scale  A  and  the  jet  velocity 
U  define  a  time  scale: 


This  time  scale  is  the  period  of  the  jet  fluctuation  as 
seen  by  an  observer  positioned  on  the  vertical  wall. 
Note  that  V / 2  is,  in  an  order  of  magnitude  sense,  the 
wave  velocity  relative  to  a  fixed  spot  on  the  wall. 


Thus,  the  fluctuation  time  of  the  wall  jet  as  an  inviscid 
stream  is: 

(3(  =  ^"*'n  (031 

*  ‘tnm  (J  fO 

where,  according  to  Eq  (21): 

A =  1.833D  (24) 

In  conclusion,  tmlll~£);  this  proportionality  appears 
as  a  straight  line  in  Fig  2,  showing  that  an  inviscid 
wall  jet  of  thickness  D  can  become  unstable  within 
a  time  interval  t> 

The  issue  of  whether  or  not  the  wall  jet  will 
become  unstable  is  decided  by  examining  the  'invis- 
cidity’  of  the  flow.  In  viscidity  or  viscidity  is  a  flow 
property,  not  a  fluid  property.  If  the  wall  jet  tends  to 
fluctuate  (wave),  then  jet  fluid  will  tend  to  make 
contact  with  the  solid  wall  and  the  adjacent  semi¬ 
infinite  fluid  reservoir  intermittently,  at  time  intervals 
t  >  tmin.  The  wall  jet,  as  a  flow,  remains  inviscid  if 
during  each  interval  t  it  cannot  learn  by  viscous 
diffusion  of  the  presence  of  a  restraining  ambient. 
The  characteristic  time  of  viscous  penetration  (tv) 
from  the  wall  and  from  the  outer  edge  of  the  jet  to 
the  jet  centreline  (over  a  distance  D/2)  is  given  by 
the  solution  to  Stokes’  first  problem': 

D/2 

T7=~  1  (25) 

2v  vt„ 


Locally,  no  jet  will  remain  inviscid  forever.  Fig  2 
shows  that  if  the  fluctuation  time  exceeds  the  viscous 


STABLE  (inviscid)  •» — .  UNSTABLE  ( inviscid ) 


/INVISCID- 


VISCID 


Transition 


VISCID 

v(stabie)N> 


Fluctuation  time 


Fig  2  Internal  competition  between  two  characteris¬ 
tic  times,  the  minimum  period  for  inviscid  instability 
(tm,„)and  the  viscous  communication  time 6  (t  J 
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communication  time  f„,  the  jet  will  remain  laminar 
because  it  will  continue  to  be  restrained  viscously  by 
its  ambient. 

The  instability  predicted  by  so  many  hydro- 
dynamic  stability  studies  is  therefore  possible  only  if 
fv  exceeds  tmin.  In  Fig  2,  this  condition  corresponds 
to  the  intersection  of  the  tmin  ~  D  line  with  the  tv  —  D 1 
line: 

=  K  (27) 

The  phenomenon  of  transition  to  nun-laminar  flow 
appears  to  be  governed  by  the  time  criterion: 


0(.V„)=1 

where: 


.Y 


a  — 


DU 
58  v 


(28) 


(29) 


The  object  of  the  remaining  presentation  is  to  test  the 
.Yu  —  1  criterion  against  the  voluminous  experimental 
record  available  on  transition  in  natural  convection 
along  a  vertical  heated  wall.  Note  at  this  point  that 
the  A('b~I  criterion  is  equivalent  to  Re  —  58 »  1 , 
where  the  Reynolds  number  Re  is  based  on  the  local 
U,  D  scales  of  the  buoyant  wall  jet.  Experimentally 
(as  described  later)  it  is  found  that  the  transition  cor¬ 
responds  to  a  Re  constant  considerably  greater  than 
unity  {Re  —  10").  The  contribution  of  the  time  scale 
argument  that  led  to  Eq  (29)  is  to  predict  a  transition 
Re  much  greater  than  unity;  unlike  the  present  argu¬ 
ments,  classical  scaling  arguments  regarding  the  rela¬ 
tive  size  of  viscous  and  inertial  terms  in  the  Navier- 
Stokes  equations  or  in  the  Orr-Sommerfeld  equation 
reveal  Re  —  1  as  a  critical  dimensionless  parameter. 

It  is  also  worth  noting  that,  theoretically,  the  same 
transition  criterion  ( iVB  ~  1,  or  Re  —  10")  is  recommen¬ 
ded  by  the  buckling  (meandering)  property  of  inviscid 
streams8.  This  coincidence  arises  because  the  buck¬ 
ling  wavelength  of  a  two-dimensional  inviscid  jet8, 
AB=  W3£)=  1.81D,  is  practically  the  same  as  the 
minimum  wavelength  for  inviscid  instability, 
Amin=1.83D.  It  has  been  shown8  that  the  NB~l 
criterion  anticipates  correctly  the  transition  to  tur¬ 
bulence  in  free  jet  and  wake  flow. 


Scale  analysis  of  natural  convection  along  a 
vertical  heated  wall 


To  be  able  to  apply  the  time-scale  criterion,  a  theoreti¬ 
cal  understanding  of  the  two  wall  jet  scales  ( U,  D)  is 
essential.  Consider  the  flow  near  a  vertical  wall,  driven 
by  the  temperature  difference  AT  between  wall  and 
fluid  reservoir.  In  general,  the  flow  thickness  ( D )  will 
differ  from  the  thickness  of  the  fluid  layer  heated  by 
the  wall9  (5t).  The  boundary  layer-approximated 
equations  governing  the  conservation  of  mass, 
momentum  and  energy  in  the  system  shown  in  Fig  3 
are1”: 


bu  b  v  „ 
— +  — =  0 
,ix  by 


bu  bu  b“ u  _  , 

u  — +  v—-=v  —  +  gP(T-Tx) 
bx  by  by 


(30) 

(31) 


Wiiere  x,  y,  u,  c,  T,  v,  g,  f3  and  a  are  the  co-ordinates, 
velocity  components,  temperature,  kinematic 
viscosity,  gravitational  acceleration  in  the  negative 
x  direction,  coefficient  of  thermal  expansion  and 
thermal  diffusivity,  respectively. 

Let  8t  be  the  thermal  boundary  layer  thickness, 
ie  the  slender  fluid  region  in  which  the  wall  heating 
effect  is  felt.  In  this  region,  the  heat  conducted 
horizontally  from  the  wall  into  the  fluid,  represented 
by  the  scale  a\T/8j  from  the  energy  equation  (Eq 
(32)),  is  converted  into  enthalpy  flow  in  the  vertical 
direction,  uAT/H.  Thus,  the  balance  between  con¬ 
duction  and  convection  in  the  layer  of  thickness  ST 
requires  the  following  equivalence  between  the 
corresponding  scales: 


uAT  AT 

H  5^ 


(33) 


or: 

u  —  aH/&\  (34) 

The  momentum  equation  (31)  accounts  for  the  com¬ 
petition  between  three  forces:  inertia,  friction  and 
buoyancy.  The  scales  of  these  forces  are,  in  order: 

u 1  u 

—  i'tt ,  g/3  AT  (35) 

n  o  j 

Assuming  first  that  the  effect  of  inertia  is  negligible, 
and  that  Eq  (31)  is  a  balance  between  buoyancy  and 
friction,  we  write: 

p^r-g(i±T  (36) 

5t 


r 

u 


Pr>  1 
Scaling 


Scaling 


u 


bT  bT 

- f  c  — 

bx  by 


Fig  3  Relative  sizes  of  thermal  and  velocity  boun¬ 
dary  layers  for  high  and  low  Pr  fluids 
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Using  Eq  (34)  we  find  that,  in  a  ST  layer  dominated 
by  a  balance  between  buoyancv  and  friction,  the 
scales  are: 


8r~  HRa,,'  4 
u-aRa)r/H 


The  Rayleigh  number  based  on  vertical  length  scale 
H,  RuH,  is  defined  as 

g fiSTH' 


(38) 


With  the  same  assumptions  that  led  to  the  scales  above 
( Eq  (37)),  the  negligibly  small  ratio  inertia/friction, 
or  inertia/buoyancy,  requires 


u2/H 

u/Sj 


(39) 


In  conclusion,  the  scales  are  valid  in  high  Prandtl 
number  fluids  (v/a»  1):  as  shown  on  the  left  side  of 
Fig  3,  Pr  »  1  fluids  develop  an  additional  length  scale 
D.  which  is  the  thickness  of  the  wall  jet  referred  to 
in  the  time  criterion  for  transition  (Eqs  (28,  29)).  The 
D  scale  follows  from  the  momentum  equation  scales 
(Eq  (35)).  Outside  the  thermal  boundary  layer  the 
fluid  is  isothermal,  hence,  the  buoyancy  force  is  negli¬ 
gible  relative  to  both  inertia  and  friction.  The 
equivalence  of  inertia  and  friction  scales  in  the  layer 
of  thickness  D: 


u“  u 
—  —  v — ; 

H  £>- 


(40) 


in  conjunction  with  scales  (37)  yields: 

D  —  Pr1  ~8r  —  HPr1'2  Raj{1  4  (41) 

for  Pr»l.  The  size  of  the  wall  jet  relative  to  the 
thermal  boundary  layer  is  shown  schematically  in  Fig 
3.  The  D  scale  (Eq  (41)),  will  be  used  later  in  the 
application  of  the  transition  criterion. 

It  remains  to  establish  the  u,  5r  and  D  scales 
prevailing  in  the  case  where  the  momentum  equation 
represents  a  balance  between  buoyancy  and  inertia 
in  the  layer  of  thickness  5T  (note  that  the  Pr  »  1  scales 
(37)  and  (41)  are  based  on  a  friction  ~  buoyancy 
balance  in  the  5T  layer).  Writing: 

jf-gP^T  (42) 

and  using  Eq  (34)  yields: 

8T~H(Ra„Prr1'* 

u~a(Ra„Pr)w2/H 

The  inertia  ~  buoyancy  balance  governs  the  ST  layer, 
except  in  a  layer  immediately  adjacent  to  the  wall 
where  the  vertical  velocity  u  obeys  the  no-slip  condi¬ 
tion.  In  the  5V  layer  friction  is  an  important  effect, 
hence: 


(44) 

o~ 

Combining  this  result  with  the  u  scale  given  by  Eq 
(43)  yields: 


<$,  —  ///fa H 1  4  Pr'  4  —  5T  Pr 1  2 


(45) 


The  relative  magnitude  of  the  6X  and  5,  scales  is 
shown  in  the  right  half  of  Fig  3.  Thus,  we  draw  two 
important  conclusions:  that  the  flow  scales  (43,  45) 
which  follow  from  a  balance  between  inertia  and 
buoyancy  in  the  8r  layer  correspond  to  low  Prandtl 
number  fluids;  and  that  the  thickness  of  the  wall  jet 
in  this  case  is  the  thickness  of  the  entire  layer  heated 
by  the  wall. 

D~-8r~IHRa„  Pr)  1  4  (46) 

for  Pr«  1.  In  what  follows  we  shall  rely  on  the  above 
dimensions  to  translate  the  time-scale  criterion  (28. 
29)  into  the  terminology  in  which  the  phenomenon 
of  transition  has  been  recorded  by  previous  studies  of 
natural  convection. 


Transition  along  a  vertical  heated  wall 

A  relatively  wide  selection  of  experimental  observa¬ 
tions  on  the  beginning  of  transition  along  an  iso¬ 
thermal  wall  was  compiled  already  by  Mahajan  and 
Cebhart11.  Table  1  reproduces  this  compilation  and 
shows  the  number  G  above  which  the  burn  ant  wall 
jet  was  noticed  to  become  nonlaminar.  The  number 
G  is  defined  as: 

in  other  words: 

G  =414  Raj, 4  PG'  4  (47) 

The  time  criterion  (28.  29)  can  be  rewritten  in  terms 
of  G,  by  using  scales  (37,  41)  for  Pr  >  1  fluids;  taking 
U  ~u  we  have: 


1,4 

H 


Pr'1'2 


C 

164  Pr1  4 


(48) 


Table  1  shows  the  Xa  value  corresponding  to 
each  experimental  report:  in  all  cases  0(.V„)=1, 
which  is  the  same  as  the  time  criterion  for  transition. 
In  conclusion,  the  experimental  data  on  transition 
along  isothermal  walls  supports  the  theoretical  argu¬ 
ment  that  the  transition  phenomenon  is  marked  by 
the  equivalence  of  time  scales  tmin—  f,. 

Table  2  is  a  compilation  of  transition  observa¬ 
tions  made  using  a  vertical  wall  with  constant  heat 


Table  1  Experimental  observations  on  beginning  of 
transition  along  a  vertical  isothermal  wall 


Pr 

G 

"b 

Warner  and  Arpaci’2 

0.72 

466 

308 

Colak-Antic'3 

0.72 

572 

3.79 

Cheesewright'4 

0  72 

600 

3.97 

Regnier  and  Kaplan15 

0.72 

622 

4.12 

0.77 

460-547 

299-3.56 

0.77 

645-702 

4  20-4  57 

0.77 

541 

352 

0.77 

605 

394 

077 

378 

2  46 

Eckert  and  Soehngen16 

0.72 

400 

2  65 

Hugot  ef  a/17 

07 

724 

483 

Szewczyk1* 

07 

665 

4  43 

6.7 

534 

2  02 
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flux  (/".  Tilt*  observations  have  been  catalogued  in 
terms  of  Pr  and  the  number  G*  defined11  as: 


In  order  to  obtain  the  relationship  between  jV„  and 
C*.  we  make  the  observation  that  for  Pr>  1  fluids 
the  thermal  boundary  layer  thickness  and  vertical 
velocity  scale'1  as: 


Substituting  D  -  ST  Pr'  2  and  V  -  u  into  Eq  (29) 
yields: 


x  G 1 

‘  b"'-210 Pr>  10  {o2) 

The  NB  values  corresponding  to  the  experimental 
observations  are  listed  in  the  last  column  of  Table  2. 
Once  again,  the  measured  -VB  number  is  of  order  one, 
in  agreement  with  the  theoretical  criterion. 

Fig  4  summarizes  Tables  1  and  2.  The  experi¬ 
mental  observations  on  transition  fall  consistently  in 
the  0(.VB)=  1  domain  in  the  Prandtl  number  range 
0.7-11.4. 


Table  2  Experimental  observations  on  beginning  of 
transition  along  a  wall  with  constant  heat  flux 


Pr 

G 

% 

Mahajan  and  Gebhart" 

0.71 

388-620 

2.05-3.27 

0.71 

400-650 

2.11-3.43 

Jaluna  and  Gebhart'9 

6.7 

504-802 

1.36-2.16 

Godaux  and  Gebhart20 

6.7 

563-802 

1.52-2.16 

6.7 

528-979 

1.42-2.63 

Vliet  and  Liu2' 

6.2 

855 

2.36 

6.4 

955 

2.61 

6.4 

900 

246 

5.05 

960 

2.81 

Lock  and  Trotter22 

11.0 

293 

066 

11.4 

368 

0.84 

Fig  4  Variation  of  time  ratio  NH  with  Pr 


Transition  in  enclosures  heated  from  the  side 

In  the  preceding  section  we  saw  that  the  experimental 
observations  on  transition  accumulated  over  the  past 
three  decades  support  the  idea  that  the  fluctuation 
(waving)  time  flnm  matches  the  time  of  viscous 
diffusion  across  the  wall  jet.  Now,  by  focusing  on 
enclosures  heated  from  the  side,  we  have  the  oppor¬ 
tunity  to  actually  see  the  incipient  waviness  of  the 
two  wall  jets  during  transition.  The  visualisation  of 
the  wall  jet  wavelength  A  is  made  possible  by  the 
slow  motion  of  the  ‘core’"',  ie  the  motion  of  the  cavity 
fluid  sandwiched  between  the  two  fast-moving  wall 
jets.  It  has  been  discovered  experimentally2'  that  in 
a  characteristic  Rayleigh  number  range  the  core  fluid 
engages  in  a  cellular  motion  of  the  type  shown 
schematically  in  Fig  5.  A  large  number  of  experi¬ 
mental  and  numerical  studies  have  confirmed  this 
phenomenon,  especially  the  fact  that  the  number  of 
core  cells  increases  as  the  Rayleigh  number  increases. 
Table  3  shows  a  representative  sample  of  experimental 
and  numerical  observations. 

In  view  of  the  theoretical  discussion  presented 
earlier,  it  is  reasonable  to  regard  the  cellular  structure 
of  the  slow-moving  core  as  the  reflection  (the  finger¬ 
print)  of  the  waviness  acquired  by  the  wall  jets  during 
transition.  Thus,  the  eell-to-cell  distance  visible  in 


Fig  5  Representation  of  buoyancy  induced  flow  in 
a  cavity  of  large  aspect  ratio  with  one  wall  heated, 
the  other  cooled,  and  the  top  and  bottom  insulated 
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scale  predicted  by  Eqs  (54  and  55).  The  best  agree¬ 
ment  between  (A,,,,,,),,,,  and  occurs  at 

extreme  Prandtl  numbers  (Pr  =  0.035  and  Prs=  480): 
This  is  a  direct  consequence  of  the  fact  that  scales  (54 
and  55)  are  valid  strictly  in  the  limits  Pr  -*oc  and 
Pr-*0,  respectively.  The  in-between  experimental 
results  (Pr  =  0.7)  are  least  accurately  represented  by 
either  Eq  (54)  or  Eq  (55);  nevertheless,  the  ratio 
(Amll|),.Xp/(AMllll)lhi.m%  for  Pr  =  0.7  is  practically  indepen¬ 
dent  of  Rau,  stressing  the  earlier  conclusion  that  A,,,,,, 
alwavs  scales  with  the  thickness  of  the  wall  jet. 

It  is  worth  noting  that  the  (A miJ TI)~  Ra,,'  4 
scaling  law  recommended  by  the  time  scale  criterion 
(28)  is  consistent  not  only  with  the  27  experiments 
reviewed  in  Table  3  and  Figs  6  and  7,  but  also  with 
theoretical  results  known  already  from  the  hydrody¬ 
namic  stability  analysis  of  the  vertical  enclosure  flow. 
Attention  is  drawn  to  Bergholz’s  comprehensive  study 
of  the  flow  stability  in  a  vertical  slot’11.  As  ‘base  flow’ 
for  the  stability  analysis,  Bergholz  considered  a 
counterflow  velocity  profile  independent  of  altitude, 
as  would  be  the  case  only  in  an  infinitely  tall  cavity31. 
He  then  accounted  for  the  finiteness  of  the  cavity 
aspect  ratio  H/L  by  postulating  the  existence  of  a 
constant  vertical  temperature  gradient  through  the 
slot,  S[K/m]:  note  that  in  the  vertical  boundary  layer 
regime,  the  thermal  stratification  S  is  of  order  XT / H, 
where  AT  is  the  temperature  difference  in  the 
horizontal  direction  (Fig  5).  By  increasing  the 
dimensionless  stratification  parameter  y  = 
(\S(L/XT'Rai)1’*,  Bergholz  was  able  to  make  the 
base  flow  more  jet-like,  that  is  more  like  the  vertical 
wall  jets  of  Figs  1  and  3,  which  are  known  to  prevail 
in  vertical  enclosures  in  the  boundary  layer 
regime  . 

In  Fig  7  of  his  study,  Bergholz  reported  the 
critical  wavenumber  aL =  2itL/\l.  versus  the  stability 
parameter  y  and  the  Prandtl  number.  One  very  inter¬ 
esting  aspect  of  Fig  7  is  that  in  the  boundary  layer 
limit  (■y»l)  and  in  the  high  Prandtl  number  limit 
Pr-* oc,  the  wavenumber  of  travelling  modes  is  pro¬ 
portional  to  y: 

«<*-%  (56) 

the  proportionality  constant  being  a  number  of  order 
0(1)  (in  the  same  limits,  the  wavenumbers  of  the 
stationary  modes  are,  numerically,  not  much  different 
than  those  of  the  travelling  modes,  however,  they  do 
not  appear  to  follow  the  line  represented  by  Eq  (56)). 


0  _ . _ i _ ; _ i _ i _ I 

tO*  ‘O'  *0‘  "0"  10'  'O’"  10'' 


Fig  7  Experimental/theoretical  inviscid  wavelength 
ratio  as  a  function  of  Ra  and  Pr  t  Pr  notation  as  for 
Fig  6) 


Noting  the  L-based  definitions  of  both  cr,  and  y 
employed  by  Bergholz’11,  it  is  easy  to  prove  that  the 
a,  —  y  scaling  illustrated  by  Bergholz  is  actually: 

“~Ra"'  (57) 

This  scaling,  predicted  by  the  conventional  stability 
theory 10  is  the  same  as  the  scaling  law  (53)  produced 
by  a  much  more  direct  argument  above. 

The  agreement  between  the  present  approach 
and  conventional  stability  theory  is  illustrated  further 
by  the  fact  that  aL  reaches  asymptotically  a  0(1)  con¬ 
stant  as  y  approaches  zero.  This  is  shown  by  Fig  7  in 
Bergholz1'1,  but  is  also  predicted  by  the  A ,,,,,,  —  D  scal¬ 
ing  on  which  the  present  argument  is  based.  Note  that 
as  y  decreases  towards  zero,  the  vertical  jet  thickness 
increases  and  so  does  the  wavelength  (Eq  (56);;  but 
this  process  cannot  continue  beyond  the  point  where 
the  jet  thickness  becomes  of  order  L  because,  regard¬ 
less  of  how  small  y  is,  the  slot  of  width  L  must  house 
two  jets  in  counterflow.  According  to  the  A D 
scaling,  both  D  and  Alni„  must  be  of  order  L  when 
y  <  0(1):  this  also  means  ’at.  -»  0(  1)  constant  as  y  -»  O', 
which  is  precisely  the  behaviour  unveiled  by  classical 
stability  analysis  (travelling  as  well  as  stationary 
modes.  Fig  7  of  Ref  30). 

Although  the  compatibility  between  the  results 
of  conventional  stability  theory  and  the  scaling  results 
based  on  the  argument  given  earlier  is  relevant  and 
interesting,  one  should  not  expect  one  theory  to  ‘repro¬ 
duce’  the  results  of  a  different  theory.  Unlike  stability 
theory,  which  is  a  mathematically  precise  approach, 
the  present  theoretical  argument  is  approximate  and 
based  on  the  comparison  of  scales.  This  is  why  the 
success  of  the  present  argument  should  be  measured 
in  terms  of  its  ability  to  predict  the  trends  and  orders 
of  magnitude  revealed  by  the  many  independent 
experiments  collected  in  Tables  1-3. 


Concluding  remarks 

The  published  observations  on  transition  to  tur¬ 
bulence  in  natural  convection  were  rev  iewed  in  order 
to  test  the  validity  of  the  time  scale  criterion  (Fig  2) 
for  transition  0(.VB)  =  1  or  O (Re)  =  10‘  formulated.  It 
was  found  that  laboratory  observations  and  numerical 
simulations  support  the  theoretical  viewpoint  that: 

1.  At  transition,  the  wall  jet  exhibits  a  unique 
wavelength  which  always  scales  with  the  thickness 
of  the  jet. 

2.  Transition  is  marked  by  the  equivalence  of  two 
time  scales,  both  properties  of  the  jet  region  of  the 
flow,  the  minimum  fluctuation  period  (fmin)  and 
the  time  of  viscous  diffusion  normal  to  the  jet  (fj. 

The  empirical  evidence  on  natural  convection  along 
vertical  heated  walls  suggests  that: 

3.  Transition  along  an  isothermal  wall  is  correctly 
anticipated  by  the  criterion  0(.VB)=  1.  where  .V„ 
is  given  by  Eq  (48). 

4.  Transition  along  a  vertical  wall  with  uniform  heat 
flux  is  correctly  anticipated  by  the  criterion 
0(\B)  =  1,  where  .V,,  is  given  by  F.q  (52). 
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5.  The  cellular  flow  exhibited  by  the  core  fluid  in  an 
enclosure  is  a  reflection  of  the  waviness  of  the  two 
wall  jets  during  transition.  The  cell  size,  or  the 
wall  jet  wavelength,  are  correctly  described  by  the 
scales  shown  in  Ecjs  (54)  and  (55). 

Finally,  the  theoretical  argument  leading  to  the 
0(.VH)  =  1  criterion  provides  a  theoretical  basis  for  the 
empirical  notion  that  the  transition  Reynolds  number 
is  a  flow  constant  considerably  greater  than  unity. 
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MECHANISM  FOR  TRANSITION  TO  Tl  RBLLENCE  IN 
BLO'  ANT  PLUME  FLOW 

Sniiii  o  K  i \n  k  v  jiiJ  -\ i) k i  v\  Hi  j  \n 

Department  of  Mechanical  I  neiiicei'iiic.  I'.unpu'  llov  OT 
I'nivorsiiv  of  Colorado.  UouiJer.  I  O  Mi.luu.  I  S  \ 

I  /\tt  i'ii  til  21  Jul\  1 9h2  iitul  in  h-rm  25  Jun:uir\  I'ts.u 

Abstract  Th  t?.  paper  reports  j  theoretical  ami  experiment.:!  Mtult  ofthe  lundament.il  me.  nanism  responsible 
for  transition  in  nalural  convection  plume  lime  Theoretically,  u  is  argued  that  the  Iran' toon  occurs  a  lien  ihe 
time  of  \  t'Cous  penelralion  normal  to  the  plume  becomes  comparable  "till  the  minimum  nine  period  with 
which  the  plume  ean  llucluate  as  an  unstable  inviscid  stream,  it  is  also  argued  dial  at  transition  the  plume 
wavelength  must  always  scale  with  the  local  plume  diameter  The  experimental  part  of  the  study  focused  on 
transition  in  the  axisymmetric  air  plume  above  a  point  heal  source.  Smoke  v  tsuali/ation  of  the  plume  shape  ai 
transition  led  to  extensive  observations  that  support  strongly  ihe  transition  mechanism  proposed 
theoretically.  The  transitional  plume  is  seen  to  meander  in  a  plane liwo-dimensionally  i  and  with  a  wavelength 
which  scales  with  the  plume  diameter.  If  evened  externally  by  many  such  wavelengths,  the  plume  has  the 
properly  to  select  the  natural  wavelength  proposed  theoretically.  The  equivalence  beivveen  the  present 
transition  mechanism  and  the  transition  predicted  by  the  buckling  theory  is  discussed. 


NOMUNCI.  VTL  RC 

D  local  plume  diameter  [m] 

f  disturbance  frequency  [s  '] 

(/  gravitational  acceleration  [  m  s  " :] 

H  loudspeaker  height  [ml 

k  thermal  conductivity  [\V  m  1  K  1  ] 

V  ratio  between  t.  and  same  as  buckling 
number  V „  [16.  21] 

Q  heat  input  [W] 

t  time  [s] 

rmin  minimum  plume  fluctuation  time  [s] 

(,  viscous  communication  time  [s] 

U  plume  velocity  [ms  '] 

V  transition  height  [m] 

Greek  symbols 

:i  thermal  dilTusivity  [m:  s  '1 

l>  thermal  expansion  coefficient  [K  ;] 

/.  wavelength  [m] 

/.8  buckling  wavelength  [m] 

7.min  minimum  plume  lluetuation  wavelength  [m] 
v  kinematic  viscosity  [m- s  '] 

Subscripts 

B  buckling  property 

0  reference  state 

I.  INTROlHt  TION 

Tuts  is  a  study  of  the  fundamem.il  mechanism  which 
causes  the  transition  to  turbulence  tn  buoyant  plumes 
rising  from  a  point  heat  source  The  transition  to 
turbulence  is  one ot  the  most  basic  phenomena  w  hiv  It  is 
not  vet  fully  understood  I  he  importance  of  under¬ 
standing  this  phenomenon  is  self-evident,  coin,  Icr- 
tng  the  importance  of  predicting  the  ensuing  turbulent 
motion  ot  lluids.  I  lie  transition  phenomenon  s  par¬ 
ticularly  important  m  the  held  of  heal  transfer,  env  iron- 

°Kr  ?*. 


menial  engineering,  atmospheric  research,  because 
turbulence  is  the  most  effective  transport  mechanism 
known  to  man. 

The  buoyant  plume  is  one  frequent  type  among 
the  many  occurrences  of  free-eenvection  flows  in 
engineering  and  other  applications.  In  this  paper  we  are 
specifically  interested  in  the  axisymmetric  buoyant 
plume  rising  from  a  point  heat  source  in  a  quiescent 
environment.  Existing  studies  on  such  plumes  and  the 
plume  transition  phenomenon  have  been  summarized 
by  Gebhart  [I],  who  showed  that  considerable  effort 
has  been  devoted  to  this  problem  over  the  past  few 
decades.  For  example,  the  laminar  2-dim.  and 
axisymmetric  plume  has  been  studied  by  Yih  [2.  ?].  and 
Brand  and  Luhey  [4],  Probably  the  most  thorough 
treatment  of  this  problem  is  the  numerical  analysts  bv 
r  um  [5], 

Experimental  studies  on  the  plane  plume  above  the 
line  source  were  reported  by  Brodowtc/  and  kierkus 
[ti],  Forstrom  and  Sparrow  [7],  and  Schorr  and 
Gebhart  [X],  The  latter  two  works  are  concerned  with 
the  plume  behavior  in  the  transition  regime  as  well  as  in 
the  laminar  regime.  Schorr  and  Gebhart  observed  by 
means  of  interferometric  flow  visualization  a  regular 
laminar  'swaying'  motion  at  a  large  distance  above  a 
line  heat  source.  This  type  of  boundary  layer  swaving 
motion  is  amplified  and  eventually  the  llovv  becomes 
turbulent.  Forstrom  and  Sparrow  also  observed  the 
existence  of  swaying  motion  at  a  fixed  point  in  space 
near  transition,  by  means  of  a  thermocouple  placed 
midway  between  the  mid-plane  and  ihe  edge  of  the 
thermal  houndarv  layer 

The  theoretical  research  >>n  transition  in  plume  llovv 
proceeded  along  the  lines  of  hydrodynamic  stability 
theory  for  example.  Peru  and  Gebhart  ( v>  ]  showed  by 
integrating  mviscid  cases  of  the  On  Somnierlelt  type 
equation  that  the  assumed  2-dim  plume  base  llovv  is 
less  stable  for  the  asymmetric  mode  ill. in  for  the 
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symmetric  one.  A  thorough  study  of  buoyancy  elTect  on 
hydrodynamic  stability  in  the  vertical  round  jet  has 
been  conducted  by  Vlollendorf  [10]  who  found  that 
buoyancy  strongly  alTects  the  stability  of  jets. 

In  the  discussion  of  existing  information  on 
transition,  the  adopted  theoretical  view  is  rooted  in  the 
theory  of  hydrodynamic  stability.  This  point  of  \:evv 
consists  of  recognizing  the  existence  of  external 
disturbances  of  many  wavelengths  which,  when 
superimposed  on  the  laminar  (low  of  interest,  might 
render  the  flow  unstable,  thus  triggering  turbulence. 
The  same  comment  applies  to  the  existing  experimental 
work  on  transition:  this  time  the  disturbances  are 
introduced  into  the  flow  externally,  for  example,  by 
using  a  loudspeaker.  Thus,  as  we  look  back  at  the 
important  advances  made  in  this  field  of  transition 
research,  it  is  important  to  keep  in  mind  that  this 
research  is  not  about  tlows  alone,  but  about  the 
response  of  certain  tlows  to  certain  disturbances. 

From  the  outset,  it  must  be  said  that  the  point  of  view 
which  stimulated  the  present  study  differs  somewhat 
from  the  classical  hydrodynamic  stability  approach. 
The  difference  lies  in  the  fact  that  in  this  study 
’transition-  is  viewed  as  an  intrinsic  property  of  the  flow 
alone,  i.e.  a  property  which  is  not  related  to  the  nature 
or  questionable  presence  of  an  external  disturbance.  It 
is  shown  in  the  next  section  that  this  intrinsic  property 
stems  from  hydrodynamic  stability  results  which  have 
been  known  for  one  hundred  years.  However,  it  is 
apparent  that  the  significance  of  these  results  vis-a-vis 
transition  has  not  been  emphasized  untii  now. 

1  THE  MINIMUM  WAVELENGTH  FOR  INVISCID  FLOW 

INSTABILITY:  TIME  CRITERION  FOR  TRANSITION 

The  type  of  classical  stability  results  which 
stimulated  the  present  study  is  exhibited  in  Table  1.  In  a 
1 880  paper.  Rayleigh  [11]  showed  that  an  inviscid  jet  of 
triangular  profile  is  unstable  to  disturbances  whose 
wavelengths  exceed  a  certain  multiple  of  the  jet 
thickness  D  Rayleigh  did  not  calculate  explicitly  the 
minimum  wavelength  of  ’neutral’  stability:  his 
discussion  focused  primarily  on  another,  longer 
wavelength  ( 2 \D)  for  which  an  assumed  disturbance 
exhibits  the  highest  amplification  rate  ([1 1  ],  p.  65).  The 
minimum  wavelength  for  instability.  /.min,  which  results 
from  Rayleigh’s  jet  analysis  is  listed  in  rowfalofTable  1. 

Similar  results  have  been  known  from  stability 
studies  involving  other  basic  flows.  For  example. 
Rayleigh  considered  also  the  free  shear  flow  profile  (b) 
and  found  instability  for  wavelengths  greater  than  ’5 D' 
([11],  p.  65).  Thus,  for  several  velocity  profiles  of  the 
base  flow,  the  minimum  wavelength  for  inviscid 
instability  always  scales  with  the  transversal  dimension 
of  the  flow.  As  shown  in  row  (cl  of  Table  I.  this  scaling  is 
consistent  with  another  classical  result,  namely,  the 
instability  of  a  plane  of  velocity  discontinuity  to  any 
wavelength  [  I  2]  In  thisea.se.  the  minimum  wavelength 
is  zero.  i.e.  oft  he  same  order  as  the  shear  layer  thickness 
The  proportionality  /„„„  -  /)  is  also  encountered  in  the 
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stability  analysis  of  radially  symmetric  flows,  for 
example,  round  jets  [  1 3]  and  annual  shear  layers  [  1 4], 

The  object  of  this  study  is  the  transition  to  turbulence 
in  a  buoyant  plume.  This  flow  is  represented 
approximately  by  profile  (ai  in  Table  1.  The  theoretical 
basis  of  the  present  research  is  the  idea  that  the  zmin  -  D 
scaling  discussed  in  the  preceding  paragraphs  is  an 
intrinsic  property  of  the  inviscid  flow,  and  that  this 
property  is  responsible  for  transition.  The  mission  of 
the  experimental  work  outlined  later  in  this  paper  is  to 
verify  the  validity  of  this  theoretical  viewpoint.  Below,  a 
simple  scaling  argument  is  offered  as  a  basis  for  the 
transition  phenomenon,  and  as  an  analytical  result  to 
be  verified  by  experiment. 

Each  longitudinal  length  scale  z.(>z.mln)  and  the 
plume  velocity  U  define  a  new  time  scale. 


This  is  the  period  in  which  the  stream  will  fluctuate 
relative  to  the  still  environment.  Note  that  C/  '2  is  the 
plume  mean  velocity  which,  from  symmetry  consider¬ 
ations.  represents  the  order  of  magnitude  of  the  velocity 
with  which  the  /.  wave  rises.  The  same  flow  is  unstable  to 
an  infinity  of  wavelengths  /.  >  zm,n  [11],  hence,  the 
same  flow  can  fluctuate  with  an  infinity  of  periods 


However,  since  /.m,„  is  proportional  to  D  (Table  1),  the 
minimum  fluctuation  period  rmm  is  proportional  to  D 
also.  The  proportionality  imw  -  D  is  shown  as  a 
straight  line  on  Fig.  1 .  where  D  is  plotted  on  the  ordinate 
because  in  natural  convection  the  plume  becomes 
thicker  with  altitude.  For  any  inviscid  stream  of 
thickness  D.  fluctuations  with  a  period  shorter  than  tm,„ 
are  stable. 

The  issue  of  whether  the  stream  (l  .  Dl  will  become 
unstable  is  decided  by  examining  the  Ynviscidity’of  the 
flow  Inustulil  \  sir  rtsvuhl  visa  ’.'<•«  properly,  nut  a  ■ luul 
proper!  v.  If  the  stream  tends  to  fluctuate  (wave),  then 
phime  fluid  will  tend  to  collide  with  the  stagnant 
ambient  intermittently,  at  time  intervals  i  >  tmin.  The 
plume-stream  remains  inviscid  if  during  each  interval  t  it 
is  not  overcome  by  viscous  effect-,  i  c  it  does  not  learn 
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Fio.  I.  The  internal  competmon  between  the  minimum  period 
for  inviscid  instability  |rmlll)  and  the  viscous  communication 
time  (rj. 


by  viscous  diffusion  of  the  presence  of  a  restraining 
ambient.  The  characteristic  time  of  viscous  penetration 
Uv)  from  the  plume-ambient  interface  to  the  plume 
centerline  (over  a  distance  of  order  D  2)  is  given  by  the 
classical  solution  to  Stokes'  first  problem  [15] 


in  other  words. 


Thus,  at  any  level  (height)  in  its  development,  no  plume 
will  remain  inviscid  forever.  Figure  1  shows  that  if  the 
fluctuation  time  exceeds  the  viscous  penetration  time  tv, 
the  plume  will  remain  laminar  because  its  ambient  will 
successfully  continue  to  restrain  it  viscously. 

The  inviscid  instability  predicted  by  so  many 
hydrodynamic  stability  studies  (Table  I)  is  therefore 
possible  only  if  /„  exceeds  fmin.  In  Fig.  1.  this  condition 
corresponds  to  the  intersection  of  the  tmin  -  D  line  with 
the  r,  ~  Dz  parabola, 

U*  ~  fv  <51 


The  phenomenon  of  transition  to  non-laminar  plume 
flow  appears  to  be  governed  by  the  time  criterion 

OIS)  =  1  (6) 


where 


\ 
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Noting  that  D  is  a  constant  (Table  I  (.the  Ol.Vl  =  1 
criterion  is  equivalent  to  stating  that  at  transition  the 
stream  t local)  Reynolds  number  D C  >•  is  a  certain 
(critical)  constant  considerably  greater  than  unity. 
Thus,  approximating  -  I  7140  from  fable  1.  the 


transition  criterion  «>l  becomes 
DC 

>  55.  (XI 

r 

for  transition  to  non-laminar  flow 

The  obteet  of  the  following  experimental  study  is  to 
test  the  validity  of  the  r„lm  -  /,  sealing,  as  mechanism 
for  transition  in  plume  flow.  It  is  worth  noting  from  the 
outset  that  the  theoretical  time  criterion  is  already 
compatible  with  two  earlier  conclusions  regarding 
transition  : 

t!)  Lxpcrimcntally.  it  is  a  universal  conclusion  that 
transition  is  associated  with  a  characteristic  Reynolds 
constant  considerably  greater  than  unity. 

(2)  Theoretically,  the  same  transition  criterion  is 
recommended  by  the  buckling  theory  of  inviscid  jets 
[16],  where  /.m,„  is  replaced  by  the  bucklii u/  wavelength 
of  the  stream.  /.H  =  (n  2 ID  =  1.570 

An  important  distinction  must  be  made,  however, 
between  the  above  theoretical  criterion  [equations  (6)- 
sxt]  and  the  universally  accepted  fact  that  the  transition 
is  characterized  by  Di  V  =  constant.  To  begin  with, 
the  notion  of  a  ’critical'  Reynolds  number  of  transition 
is  of  purely  empirical  origin.  On  the  theoretical  side,  the 
linearized  (small-disturbance)  equations  of  hydro- 
dynamic  stability  theory  can  easily  be  subjected  to  scale 
analysis  to  show  that  the  Rey  nolds  number  is  indeed  an 
important  dimensionless  parameter:  comparing  the 
Reynolds  number  with  unity  (one),  we  can  say  w  hether 
or  not  the  v iscous  terms  can  be  neglected  in  the  stability 
analysis.  Note,  however,  that  this  scaling  argument  is 
not  about  'transition',  rather,  it  is  about  the  simpli¬ 
fication  of  stability  analysis.  Also  on  the  theoretical 
side,  the  Reynolds  number  appears  in  the  solution  to 
the  complete  Orr-Sommerfeld  equation :  however, 
unlike  in  criterion  (8),  the  stability-derived  transition 
Reynolds  number  is  not  a  constant,  lit  is  a  function  of 
the  wavelength  of  each  postulated  disturbance.) 

In  view  of  this  discussion,  the  time  criterion  (6HS) 
seems  to  provide  for  the  first  time  a  hydrodynamic 
stability  scaling  basis  for  "the  Reynolds  number  =  a 
constant  considerably  greater  than  unity”  as  transition 
criterion.  The  experimental  observations  summarized 

later  in  Section  6  show  that  at  transition  the  order  of 
«  , 

CD  v  is  10'.  in  agreement  with  the  time  criterion  (6H 81. 

3.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  experiments  focused  on  a  controlled  version  of 
the  cigarette  smoke  phenomenon  with  which  we  are  all 
familiar.  As  shown  in  Fig.  2.  the  apparatus  consisted  of  a 
man-size,  airtight.  Plexiglas  enclosure  which  was 
needed  to  isolate  the  experiments  from  ambient  air 
currents  present  in  the  laboratory  Two  adjacent  side- 
walls  and  the  top  and  bottom  walls  oft  he  chamber  were 
constructed  of  wood  The  remaining  two  adjacent  side 
walls  were  made  of  Plexiglas  in  order  to  permit  the 
lighting  and  viewing  ol  the  smoke  plume 

\n  axisymmelriv  ur  plume  was  generated  above  a 
small  heat  source  placed  in  the  center  ol  the  box. 
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Fig.  2.  Schematic  of  the  experimental  apparatus. 


approximately  45  cm  from  the  bottom.  The  heat  source 
was  constructed  in  the  form  of  a  nichrome  resistance 
coil,  as  shown  in  Fig.  2.  Electric  power  in  the  0-50  W 
range  was  measured  and  dissipated  in  the  coil :  most  of 
this  energy  was  released  into  the  buoyant  air  plume. 
(The  estimated  radiation  heat  loss  was  less  than  6°;,  at 
22  W  and  less  than  20°o  at  55  W.) 

Cigarette  smoke  was  generated  in  an  external  device 
constructed  after  a  design  described  by  Bradshaw  [17], 
The  smoke  was  released  directly  beneath  the  nichrome 
coil  very  slowly  so  that  it  did  not  affect  the  air  plume. 
This  simple  flow  visualization  technique  worked  very 
well,  and  the  plume  shape  visualized  by  the  smoke  trace 
was  photographed. 

The  experiments  were  designed  to  test  the  validity  of 
the  fmin  ~  rv  scaling  during  transition.  For  a  certain 
(reference)  power  dissipated  in  the  coil.  Q0,  the  plume 
shape  was  photographed  3-6  times.  The  photographs 
showed  statistically  the  existence  of  a  characteristic 
height  x0  and  wavelength  a0  for  the  beginning  of 
transition  (Fig.  2).  The  power  setting  Q  was  changed 
during  the  course  of  experiments  and  these  changes 
reflected  in  the  measured  x  and  L  The  object  of  the 
experiment  was  todiscover  thedimcnsionless  functions 

—  =  function  (  — ] .  (9) 

*»  \Qo) 

~  =  function  (  —  ).  (HU 

/-n  ' QoJ 

The  experimental  findings  were  then  compared  w  ith 
the  theoretical  functions  recommended  by  the  fm,„  -  t, 
scaling.  The  theoretical  functions  x  v„  and  /.  can  be 
obtained  by  recalling  that  the  diameter  and  velocity  of  a 


laminar  plume  above  a  point  heat  source  scale  as  [18] 


At  transition  we  expect  /.  -  D  and  tm,„  -  r,  [or  DU  v 
^constant.  equation  (8)].  hence,  the  theoretical 
functions  to  be  tested  are 


X 

fQ\-'2 

(13) 
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Note  that  the  scaling  represented  by  equations  (11)  and 
( 1 2)  is  valid  for  Prandtl  numbers  of  0(1)  or  greater. 

4.  RESU  LTS 

In  the  first  series  of  observations  the  plume  was 
photographed  in  the  absence  of  any  external  noise 
which  might  act  as  a  trigger  for  transition.  The  measured 
transition  heights  and  wavelengths  are  shown  in  Fig.  3  : 
both  v  and  /.  decrease  as  the  energy  content  of  the  plume 
Q  increases.  The  variation  of  /.  /.0  vs  Q  Q„  parallels  the 
theoretical  curve,  equation  1 14)  but  the  measured 
transition  heights  are  consistently  greater  than  the 
theoretical  levels.  Comparing  this  first  series  of 
observations  with  the  theoretical  expectations  [equa¬ 
tions!  13)  and  i  I4i],  we  conclude  that  in  the  absence  of 
external  triggers  the  transition  wavelength  scales  with 
the  plume  diameter,  however,  the  transition  is  delayed 
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Fig.  3.  The  transition  height  and  the  wavelength  as  a  function 
of  heal  input,  in  the  absence  of  external  noise.  O.  xx0, 
transition  height.  A.  >■  7.0,  wavelength. 


Fig.  4.  The  transition  height  and  wavelength  as  a  function  of 
heat  input.  The  enclosure  wall  was  lapped  by  a  linger.  Q,x  x0. 
transition  height.  A. '  A-  wavelength.  Q0  =  55  W,  x0  -  0.25  m. 
/.0  =  0.04  m. 


and  occurs  further  downstream  from  the  theoretical 
level. 

Considerably  more  conclusive  results  were  obtained 
by  photographing  the  plume  shapes  immediately  after 
tapping  the  enclosure  once,  with  a  finger.  This  sort  of 
noise  served  to  introduce  disturbances  of  many 
(unspecified)  wavelengths  and  amplitudes  into  the  air 
plume  flow.  Figure  4  shows  the  measured  variation  of 
and  X'.v0  with  Q/Q0 :  the  agreement  with  theory  is 
very  good  in  a  relatively  wide  range  of  power  settings 
Q'Qo ■  It  appears  that  the  stream  has  the  natural  ability 
to  filter  [19]  out  of  the  disturbance  spectrum  the 
natural  wavelength  of  transition.  The  measurements 
indicate  that  the  natural  wavelength  scales  with  the 
plume  diameter  (because  7.  These  results 

validate  the  theoretical  basis  for  adopting  tm(n  ~  f,  as 
transition  criterion  (Section  2). 

The  repeatability  of  the  above  observations  is 
demonstrated  by  the  sequence  of  photographs 
presented  as  Figs.  5(a)— (c).  These  three  photographs 
belong  to  the  same  plume,  as  the  plume  strength  Q  was 
held  constant  (Q  =  31.1  W).  The  transition  wavelength 
and  height  are  recorded  instantly  by  means  of  the 
vertical  scale  mounted  next  to  the  plume,  at  the  same 
distance  from  the  camera  [note  that  Figs.  5(aHcl  were 
taken  at  dilTerent  times,  using  different  focusing 
lengths  |.  The  photographs  show  clearly  that,  gtven  a 
plume,  the  transition  to  non-lamtnar  flow  is 
characterized  by  a  characteristic  wavelength  /.  and  a 
characteristic  height  v. 

The  relationship  between  x  and  Q  at  transition 
I  Fig  4)  is  illustrated  in  Figs.  (Sal  and  lb).  From  hi”  6(a) 
to  Fig.  6<bl  the  source  strength  Q  increases  by  almost  a 
factor  of  2 ;  correspondingly  both  /.  and  v  decrease  by  a 


factor  of  the  order  of  1/^/2  predicted  theoretically. 
Again,  transition  is  characterized  by  a  well-defined 
‘meandering’  shape  with  a  unique  wavelength  and  at  a 
unique  height. 

An  important  aspect  of  the  plume  shape  during 
transition  is  its  two-dimensionality.  We  investigated 
this  aspect  by  conducting  a  separate  series  of 
experiments  in  which  the  plume  was  photographed 
simultaneously  from  two  angles,  from  the  front  and 
from  the  side.  The  side-view  was  visible  in  a  tall  mirror 
placed  vertically  near  the  plume,  at  a  45°  angle  with 
respect  to  the  camera-plume  direction.  The  mirror  view 
appears  on  the  LHS  of  each  of  the  photographs  show  n 
in  Figs.  7(aHc). 

By  tapping  the  side  of  the  box  once,  we  had 
absolutely  no  control  on  the  plane  in  which  the  plume 
would  choose  to  meander  during  transition.  Thus,  we 
had  to  take  many  photographs  in  order  to  come  across 
a  few  cases  where  the  plane  of  deformation  happened  to 
be  nearly  perpendicular  to  the  camera-plume 
direction.  Two  such  cases  are  exhibited  in  Figs.  7(a)  and 
(b):  plume  deformation  during  transition  is  clearly  in 
one  plane.  Figure  7(c)  shows  one  of  the  many  cases  in 
which  the  plane  of  deformation  did  not  coincide  with 
either  the  camera-plume  direction  or  with  the  direction 
perpendicular  to  the  camera-plume  line :  regardless  of 
the  misalignment.  Fig.  7(c)  shows  that  the  elbows  of  the 
sinusoidal  shape  are  all  in  the  same  plane. 

The  2-dim.,  planar,  character  of  the  plume 
deformation  during  transition  [Figs.  7(a) -(c)]  is  an 
important  conclusion,  it  contradicts  the  hydro- 
dynamic  stability  iis.sumption  [13]  that  the  initial 
deformation  (disturbance)  in  free  jet  flow  is  helical 
1 3-dim.  I 
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5.  TRANSITION  IN  V  II  \R\H)M<  \l  I  Y-l  OK<  H)  1*1.1  Ml 

The  scries  of  experiments  described  in  the  preceding 
section  yielded  cvvnvincing  evidence  that  the  transition 
phenomenon  is  the  result  of  the  scale  properties 
discussed  in  Section  1  I  hcse  properties  were  studied 
further  in  another  senes  ol'experiments  where,  unlike  in 
Section  4.  the  location,  intensity  and  frequency  of  the 
'trigger'  disturbance  were  controlled.  The  disturbance 
was  generated  by  a  loudspeaker  suspended  at  a  certain 
height.  //.  next  to  the  plume  d  ig.  2).  The  face  of  the 
loudspeaker  was  covered  by  a  wooden  panel  with  a 
3  mm  <  2nd  mm  horizontal  slit  cut  into  it.  Thus,  we  were 
able  to  harmonically  force  or.iy  a  narrow  section  of  the 
rising  plume,  at  a  height  determined  by  the  position  of 
the  loudspeaker.  The  intensity  of  the  harmonic  forcing 
l relative  to  a  reference  intensity)  was  monitored  by 
measuring  the  power  needed  to  drive  the  loudspeaker. 

The  experiments  were  conducted  in  a  manner  similar 
to  what  led  to  the  observations  summarized  in  Fig.  4. 
For  a  fixed  heat  source  strength  0.  the  plume  was 
disturbed  isltakenl  at  certain  frequencies  /  of  the  same 
amphtudelthe  loudspeaker  frequency  varied,  however, 
the  maximum  travel  of  us  cone  was  held  constant).  It 
was  observed  that  the  transition  height  v  depended 
strongly  on  the  frequency  of  harmonic  excitation  i.  It 
was  found  that  there  exists  a  characteristic  frequency  / 
such  that  the  transition  takes  place  at  a  minimum 
height:  frequencies  higher  and  lower  than  this 
characteristic  value  triggered  transition  at  higher 
altitudes.  These  'resonance'  characteristics  have  been 
studied  extensively  (visually  and  photographically 
[20] )  and  are  amply  documented  in  Figs.  S(ai-(fi  In 
Figs.  8(a)  -tc)  the  loudspeaker  was  held  at  a  level  4.5  cm 
above  the  heat  source,  while  in  Figs.  8ldl-(fi  the 
loudspeaker  height  was  13  cm. 

The  minimum  transition  height  and  the  correspond¬ 
ing  wavelength  were  found  to  decrease  with  the 
increasing  heat  source  strength  Q.  These  observations 
are  summari/ed  in  Fig.  0:  they  are  nearly  identical  to 
the  results  of  Fig.  4  obtained  by  tapping  the  enclosure 
once  with  a  linger.  Thus,  the  plume  resonates,  hence,  is 
deformed  most  effectively  when  it  is  harmonically 
forced  at  us  natural  frequency,  with  a  svavelength  that 
scales  with  the  plume  thickness  at  the  transition  height. 

The  acoustic  excitation  provided  by  the  loudspeaker 
introduces  two  more  variables  in  the  experiment,  the 
loudspeaker  height  and  the  excitation  amplitude.  In 
Figs.  Xial-lf)  we  show  the  effect  of  disturbance 
amplitude  Increasing  the  disturbance  has  the  effect  of 
precipitating  the  transition,  i.e.  the  effect  of  decreasing 
the  transition  height  Raising  the  loudspeaker  from 
5.4  cm  to  I  3  cm  above  the  plume  origin,  has  the  effect 
of  placing  the  transition  further  downstream 
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[lie characteristic  wavelength  and  transition  height  (the 
same  measurements  are  listed  in  Table'  2  and  ’i  It  is 
ev  idem  that  at  transition  the  height  v  is  proportional  to 
the  characteristic  wavelength  /.  in  fact. 

v  -  tttz  tf  5» 

is  an  order-of-magnitude  lit  for  the  data  plotted  in  I  ig 
Id  rite  measured  proportionality  between  v  and  at 
transition  is  anticipated  correctly  by  the  sealing 
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argument  presented  in  the  beginning  of  this  paper  [see 
equations  (13)  and  (14)] 

In  addition,  note  (hat  the  coefficient  in  equation  1 1  s| 
is  a  number  greater  than  unin.  This  finding  validates 
the  time  criterion  of  transitions,  equation  (X).  \v hicji 
translated  into  a  transition  Reynolds  number  con¬ 
siderably  greater  than  units  f  or  example,  by  using  the 
D  and  (  scales  of  the  plume  [equations  1 1  1 1  and  ( 1 2 )]. 
the  transition  height  v  is  eliminated  Irom  equation  1 1 5i 
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Fig.  8  Effect  of  disturbance  frequency  and  amplitude  on  the 
resonance  characteristics  of  the  plume  The  voltage  1'  repre¬ 
sents  the  amplitude  of  the  input  signal  to  the  loudspeaker 
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According  to  Table  I./.  Ocanonls  he  of 'the  order  of  3 
or  greater,  hence,  the  transition  Rcsnolds  number 
IH  i  must  he  a  constant  considerably  greater  than 
units  i in  the  range  ID  ln-i  Thus,  the  present 
experimental  observations  support  qualitatively  and 
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Fig.  9.  The  transition  height  and  wavelength  as  a  function  of 
heat  input,  in  air  plumes  excited  harmonically  by  a 
loudspeaker.  O  .x  v0  transition  height:  A  '  /.,>  wavelength: 
Qo  =  21  W;  x0  -  0.35  m:  =  0.06  m. 


quantitatively  the  transition  mechanism  envisioned  in 
Sections  2  and  3  and  in  Fig.  1. 

It  is  worth  noting  that  the  local  Reynolds  number 
transition  criterion  Ol.\)  =  1  or  DV  v  =  OUO1'  stated 
in  equations  (6)  and  (8)  and  verified  above,  can  be  easily 
translated  into  the  Rayleigh  or  Grashof  number 
language  used  in  natural  convection.  In  fluids  with 
Pr  >  Ot  1 )  the  velocity  boundary  layer  thickness  scales 
as \/V  2Ra~‘  4  and  the  vertical  velocity  astx.vlRu1 
where  Ru  is  the  Rayleigh  number  (xvA). 

Therefore  DUv  =  0(10*)  means  that  at  transition  the 
Rayleigh  number  is  Ru  =  Ol  l()HPr:i.  Fortheairpluntes 
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considered  in  our  study  this  prediction  reduces  to 
Gr  =  U i  It)3),  which  the  commonly  observed  range  of 
transition  Grashof  numbers  [note  that  Gr  =  Ra  Pr). 

7  CONCLl  stoss 

This  paper  described  a  fundamental  study  of  the 
phenomenon  of  transition  to  turbulence  in  natural 
convection  plume  flow.  The  experimental  part  of  the 
study  focused  on  a  controlled  version  of  the  'cigarette 
smoke'  plume  (low.  The  experiment  demonstrated 
that : 

III  at  transition,  the  plume  assumes  a  sinusoidal 
(meandering)  shape  of  characteristic  wavelength,  and 
at  a  characteristic  height  above  the  plume  origin  [Figs. 
5lal-(c)] : 

(2)  the  transition  wavelength  scales  with  the  local 
plume  diameter  (Figs.  3.  4  and  9): 

(3)  the  transitional  meandering  shape  is  2-dim.,  i.e.  in 
one  plane  [Figs.  7(al-(cl] ; 

|4|  if  the  transition  is  triggered  by  'noise',  then  it  is  a 
plume  property  to  'filler'  out  of  the  noise  the 
characteristic  transition  wavelength  which  is  propor¬ 
tional  to  the  plume  diameter  at  the  transition  height 
(Figs.  3  and  4); 

|5)  if  ihe  transition  is  triggered  by  single-frequency 
forcing,  then  ii  is  a  property  of  the  plume  to  'resonate', 
i.e.  to  deform  most  when  disturbed  with  a  wavelength 
which  scales  with  the  plume  diameter  at  transition 
height  (Fig.  9l. 

These  experimental  conclusions  support  strongly 
the  theoretical  argument  presented  in  Section  2. 
According  to  this  argument,  the  phenomenon  of 

Table  3.  Transition  height  and  wavelength  measurements 
plotted  in  Fig.  9 
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Fig.  10  The  measured  characteristic  wavelength  and  tran¬ 
sition  height  ~j  data  from  Table  2.  5  data  from  Table  ' 
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transition  is: 

(a)  an  intrinsic  property  of  the  flow: 

(hi  the  result  of  the  internal  competition  among  two 
characteristic  time  scales  of  the  stream,  the  minimum 
fluctuation  period  /„„„  and  the  v  iscous  penetration  time 
r, : 

(c)characterized  by  a  unique  meandering  wavelength 
which  is  always  proportional  to  the  local  stream 
thickness. 

The  transition  mechanism  proposed  in  Section  2 
evolved  from  a  review  of  classical  hydrodynamic 
stability  results  (Table  I )  and  led  to  the  flow  properties 
fscaling)  embodied  in  the  time  criterion  tmin  ~  t, 
[equations  |6i  and  (7)].  The  same  scales  and  time 
criterion  have  been  brought  to  light  earlier  by  the 
buckling  theory  of  inviscid  flow  [16].  It  was  shown 
recently  [21]  that  the  buckling  time  criterion  ,VB  =0(1) 
also  predicts  correctly  the  transition  in  free  jet  flow, 
shear  flow  and  wake  flow,  as  well  as  a  series  of  turbulent 
boundary  layer  features  such  as  the  bursting  frequency 
and  the  viscous  sublayer  thickness.  Thus,  the  present 
study  reveals  an  important  relationship  (equivalence) 
between  stability  theory  and  buckling  theory  with 
regard  to  explaining  transition. 
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MECANISME  DC  LA  TRANSITION  VERS  LA  TURBULENCE 
DANS  UN  ECOULEMENT  DE  PANACHE 

Resume-  On  etudie  theoriquement  et  expcrimentalement  le  mecanisme  fondamental  responsable  de  la 
transition  dans  la  convection  naturelledun  panache.  Theoriquement,  la  transition  apparait  quand  lc  temps  de 
penetration  visqueuse  normale  au  panache  devient  comparable  a  la  periode  de  temps  minimale  de  fluctuation 
du  panache  comme  un  ecoulement  instable  non  visqueux.  On  suppose  aussi  qu  a  la  transition,  la  longueur 
d'onde  du  panache  est  en  proportion  du  diametre  local  du  panache.  La  partie  experimenlale  de  letude  est 
focalisee  sur  la  transition  dans  un  panache  axisymetrique  d'air  au-dessus  d'une  source  ponctuelle.  Une 
visualisation  par  fumee  de  la  forme  du  panache  a  la  transition  conduit  a  des  observations  extensives  qui 
soutiennent  fortement  le  mecanisme  propose  dans  la  theorie.  Le  sillage  de  transition  serpente  dans  un  plan 
(bidimensionncllementlct  avecune  longueur  d'ondequi  est  a  Techelle  du  diametre  du  panache.  S‘il  est  excite 
exterieurement  par  des  longueurs  d'onde.  le  panache  a  la  propnete  de  selcctionner  la  longueur  d'onde 
proposee  theoriquement  On  discute  ('equivalence  entre  le  mecanisme  de  transinon  propose  et  celui  de  la 

theorie  du  llambement. 


Buckling  of  a  turbulent  jet  surrounded  by  a  highly  flexible  duct 
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An  experimental  and  theoretical  study  of  the  static  buckling  of  a  downward  flowing  stream 
surrounded  by  a  highly  flexible  duct  is  reported.  It  is  shown  analytically  that  the  stream  buckles 
into  a  meandering  shape  that  depends  strongly  on  the  flow  velocity  in  the  duct.  The  amplitude  of 
the  buckled  shape  is  found  to  be  governed  by  the  magnitude  of  the  transverse  force  created  by 
nonzero  curvature  at  the  free  end  of  the  duct.  Experimental  measurements  confirm  the  buckled 
shape  described  analytically. 


I.  INTRODUCTION 

In  the  course  of  flow  visualization  experiments  with 
turbulent  water  jets  we  discovered  that  vertical  jets  sur¬ 
rounded  by  thin  rubber  ducts  naturally  buckle  into  a  stable 
shape  that  is  a  strong  function  of  the  flow  velocity  in  the 
duct.  At  small  flow  rates  the  buckling  of  the  duct  is  limited  to 
regions  near  its  free  end.  As  the  flow  rate  is  increased  the 
buckled  portion  of  the  duct  also  increases,  appearing  to 
“climb”  up  the  duct  until  its  entire  length  is  filled  with  small- 
amplitude  standing  waves.  The  equilibrium  shape  of  the 
duct  is  roughly  sinusoidal  with  a  wavelength  that  is  inversely 
proportional  to  both  flow  velocity  and  distance  from  the 
fixed  end  of  the  duct.  Three  photographs  of  a  buckled  duct 
containing  flowing  water  are  shown  in  Fig.  1.  The  flow  ve¬ 
locity  in  these  photographs  increases  from  left  to  right. 

Further  increases  in  flow  velocity  beyond  what  is 
shown  in  Fig.  1  eventually  lead  to  loss  of  stability  and  the 
initiation  of  self-excited  vibrations  of  the  duct  and  enclosed 
jet.  These  vibrations  have  been  the  major  focus  of  previous 
studies  using  cantilevered  ducts  containing  flowing  fluid. 
The  phenomenon  of  static  buckling  has  not  received  as  much 
attention  as  the  self-excited  vibrations,  possibly  because  of 
misunderstandings  concerning  the  conditions  under  which 
buckling  develops.  For  example,  Thompson'  reports  that 
cantilevered  ducts  which  are  initially  straight  cannot  buckle 
under  the  influence  of  flowing  fluids.  He  is  correct  in  theory 
but  not  in  practice,  because  real  ducts  are  never  perfectly 
straight.  As  shown  below,  any  slight  initial  curvature  of  the 
duct  causes  a  nonzero  transverse  force  at  its  free  end  and 
leads  to  buckling,  provided  that  the  flow  velocity  is  large 
enough. 

Interest  in  the  flow-induced  instabilities  of  ducts  con¬ 
veying  fluids  has  been  largely  motivated  by  the  practical 
problems  of  pipeline  and  heat  exchanger  design.  Niordsen2 
was  one  of  the  first  to  derive  the  correct  form  of  the  equations 
of  motion.  In  the  course  of  a  study  of  flow  through  segment¬ 
ed  pipes  with  flexible  joints  Benjamin3  4  confirmed  Niord- 
sen’s  derivation  by  using  Hamilton's  principle.  Gregory  and 
Paidoussis5  6  and  Paidoussis7  extended  Benjamin’s  analysis 
to  the  case  of  continuous  elastic  tubing.  Benjamin's  study, 
like  those  of  Gregory  and  Paidoussis,  is  devoted  primarily  to 
stability  analyses  aimed  at  determining  the  onset  of  self-ex¬ 
cited  vibrations.  Bishop  and  Fawzy*  observed  static  buck¬ 


ling  during  their  experiments  with  forced  oscillations  of  ver¬ 
tical  ducts  filled  with  flowing  fluid,  but  reported  their 
observations  without  attempting  to  explain  them.  More  re¬ 
cently  Lundgren,  Sethna,  and  Bajaj9  observed  large-ampli¬ 
tude  buckling  of  a  cantilevered  duct  when  they  attached  an 
inclined  nozzle  to  the  free  end  of  the  duct.  They  calculated 
the  static  solution  for  the  buckled  shape  which  resulted  from 
the  attachment  of  the  nozzle,  and  examined  the  stability  of 
the  static  solution.  They  found  that,  with  the  nozzle  at¬ 
tached,  the  orientation  of  self-excited  oscillations  was  no 
longer  random  but  occurred  either  in  the  plane  of  the  nozzle 
or  perpendicular  to  that  plane. 


(°>  (b)  (c) 


FIG  I  Static  buckling  of  a  vertical  hose  with  a  diameter  of  12.7  mm  lai 
y=  1  9m/sec.  a0  =  2.78  Ibi  r  =  3  15  m/sec.a#  =  I  01  ici  1=  4  *3  m, 
sec.  a0  =  0  45 
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Our  objective  in  the  present  study  is  to  focus  on  the 
small-amplitude  static  buckling  which  is  exhibited  by  a  tur¬ 
bulent  jet  which  is  bounded  by  a  highly  flexible  duct.  This  is 
the  “flow-dominated”  regime  of  the  pipe  vibration  pheno¬ 
menon  studied  by  earlier  investigators.  In  the  analysis  pre¬ 
sented  below  we  will  include  a  gravitational  tension  term  not 
considered  by  Lundgren  et  al?  in  their  study  of  horizontal 
ducts.  This  gravitational  term  causes  the  wavelength  of  the 
buckled  duct  to  decrease  with  increasing  distance  from  the 
fixed  end  of  the  duct.  In  addition  we  will  show  that,  if  the 
flow  through  the  duct  is  considered  to  be  invisctd,  then  the 
effective  bending  stiffness  of  the  duct/jet  system  is  increased 
by  nonuniform  pressure  and  tension  distributions  created  by 
the  local  curvature  of  the  duct.  After  deriving  the  equilibri¬ 
um  equation  describing  the  system,  we  will  find  perturbation 
solutions  validforshort,flexibleducts  in  the  limits  V2  =  0{e) 
and  V2  —  O  ( 1/e).  We  will  conclude  the  paper  by  comparing 
the  wavelength  predictions  made  by  a  perturbation  solution 
appropriate  for  long  ducts  to  experimental  measurements 
derived  from  photographs  such  as  those  shown  in  Fig.  1. 

II.  MATHEMATICAL  FORMULATION 

The  basic  geometry  of  the  problem  considered  in  this 
paper  is  shown  in  Fig.  2.  A  flexible  tube  of  radius  r  and  length 
L  houses  a  stream  of  mean  velocity  V  and  density  p.  The 
centerline  of  the  flexible  duct  is  described  by  the  curve  y<.x). 
The  duct  is  assumed  to  undergo  small-amplitude,  beam-like 
deflections  with  no  change  in  the  shape  of  its  cross  section. 
These  deflections  are  assumed  to  be  small  enough  so  that 
secondary  flows  can  be  ignored.  The  fluid  mass  per  unit  of 
duct  length  is  mf ,  while  the  wall  mass  per  unit  length  is  m.  As 
shown  in  Fig.  2,  the  average  pressure,  tension  force,  shear 
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force,  and  bending  moment  are  represented  by  ptl,  F,  S,  and 
M.  Finally,  in  accordance  with  the  notation  employed  in 
membrane  stress  analysis,10  we  write  t,  sm,  and  s,  for  wall 
thickness,  meridional  stress,  and  tangential  stress,  respec¬ 
tively. 

A  force  balance  in  the  vertical  direction  (jc)  yields 
[mf  +  m)g  dx  +  ( PfyA  —  T  )cos  6 

+  [F+  dT -  A  I p„  +  dpui ]cos(0  +  dd )  =  0.  |1) 

For  small  amplitudes,  cos  [dd )  =:  cos(#  +  dd )  *  1 ,  and  Eq.  ( 1 ) 
can  be  simplified  to  the  form 

(mf  +  m]g  =  -^-(p0A-T).  (2) 

dx 

Integrating  Eq.  (2)  we  obtain 

[mf  +  m\g(x  —L)  =PoA  —  F.  (3) 

The  constant  of  integration  in  Eq.  (3)  was  evaluated  by  ap¬ 
plying  the  free-end  condition  F  =  p0  =  0  at  x  =  L. 

The  fluid  velocity  in  the  horizontal  (y)  direction  is 

£l>Wnil-££-r£,  m 

dt  dx  dt  dx 

and  the  acceleration  of  the  fluid  in  the  horizontal  direction  is 

jjrbM'lll  =  v‘^r-  HI 

Applying  Newton’s  second  law  of  motion  and  ignoring  sec¬ 
ond-order  terms,  we  find 


dy  _ 


d  /_  ..  dy  \  dS 

—  (,7--^)  —  )-  —  . 


FtG  1  Forces  and  moments  acting  on  a  duct  element 
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Substituting  for  the  pressure  p0  and  tension  T  from  Eq.  (3), 
reduces  Eq.  (6)  to 

£•  (7) 

It  is  instructive  at  this  point  to  examine  the  physical 
meaning  of  each  of  the  terms  in  Eq.  (7).  As  the  fluid  moves 
through  the  duct  it  undergoes  centripetal  acceleration  by 
virtue  of  the  fact  that  it  is  following  a  slightly  curved  path. 
The  centripetal  term  on  the  left-hand  side  of  the  equation  is 
balanced  on  the  right-hand  side  by  gravitationally  induced 
wall  tension  and  by  a  transversal  shear  force.  In  the  limiting 
case  of  zero  fluid  flow  and  zero  bending  stiffness  (i.e.,  zero 
shear  force),  Eq.  (7)  reduces  to  the  equation  describing  a 
hanging  chain.  A  surprising  implication  of  Eq.  (7)  is  that  the 
static  equilibrium  of  the  duct  does  not  depend  on  the  viscos¬ 
ity  of  the  fluid  flowing  through  the  duct.  This  result  is  due  to 
the  fact  that  the  tension  induced  in  the  duct  walls  by  fric¬ 
tional  fluid  drag  is  exactly  offset  by  the  pressure  force  re¬ 
quired  to  force  the  fluid  through  the  duct. 

The  transversal  shear  force  S  can  be  evaluated  by  re¬ 
quiring  rotational  equilibrium  about  the  bottom  face  of  the 
duct  segment  shown  in  Fig.  2. 

?  dM 

$  =  — .  (81 
dx 

where  M  is  the  total  bending  moment  acting  over  the  duct 
cross  section.  As  shown  next,  this  bending  moment  is  due  to 
three  independent  effects:  (a)  the  asymmetric  distribution  of 
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fluid  velocity  and  pressures;  (b)  the  asymmetric  distribution 
of  tangential  stresses  in  the  wall  caused  by  jet  flow;  and  (c)  the 
El  bending  stiffness  of  the  tube,  predicted  by  classical  beam 
theory.  Since  in  the  limit  of  small-amplitude  deflections 
these  effects  are  not  coupled  with  one  another,  each  will  be 
coqsidered  separately. 

In  an  in  viscid  jet  with  an  infinitely  large  radius  of  curva¬ 
ture  R  s ,  the  velocity  and  pressure  distributions  are" 

V[w\=V[l-w/R„),  (9) 

p{w)=p0+pV2w/Rx.  (10) 

The  “jet”  bending  moment  associated  with  these  asym¬ 
metric  profiles  is11 

M,=pV%^L.  (11) 

dx~ 

In  Eqs.  (9)  and  1 10),  coordinate  w  is  measured  radially  from 
the  jet  centerline,  away  from  the  centerline  of  curvature  (see 
Fig.  2).  The  symbol  I0  denotes  the  area  moment  of  inertia  of 
the  jet  cross  section  (for  a  round  jet,  I0  =  irr4/ 4). 

The  membrane  stresses  induced  in  the  wall  by  the  pres¬ 
sure  variation  across  the  duct  are  shown  in  the  bottom  por¬ 
tion  of  Fig.  2.  The  horizontal  components  of  these  mem¬ 
brane  stresses  must  balance  the  internal  pressure  force 
acting  on  the  duct  wall 

s„  t  ds  d6m  +s,tds  d9,  =  p  dr.  ( 12) 


—  =  dd„ 


Eq.  (12)  becomes 


The  relation  between  R,  and  R  is  shown  in  Fig.  3.  The  hoop 
stress,  s„ ,  can  be  found  by  integrating  the  pressure  across  a 
section  of  the  duct  as  is  shown  in  Fig.  3, 

[R^  +  w)smt  cos<{>  =  J  jpn  +  (  ~  )](*«  -I -w]dw'. 


Performing  this  integral  and  recognizing  that 


-V 


FIG.  3.  Force  balance  for  calculating  the  hoop  stress  sm . 
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w0  =  r  cos  <t>,  (18) 

yields 

sm  ={r/t  )[pQ  +  (r  cos  <t>  /2R »)\pV2  - p0)).  (19) 

The  tangential  stress  distribution  s, ,  which  is  needed  for  cal¬ 
culating  the  bending  moment,  follows  immediately  from 
Eqs.  (15)  and  (10), 

s,  =  [r/2t][pV2  -4- />,,)(  1  4-  rcosd/R^  ).  (20) 

The  “membrane”  contribution  to  bending  moment,  M„,  is 
therefore 


Mm  =2  s,tr  cos  <b  d<t>, 
Jo 


Mm  =^{pV2+p0)^2IQ\pV2+p0)*JL .  (22) 

R  v  dx 

So  far,  the  first  two  contributions  to  the  bending  mo¬ 
ment,  Eqs.  (1 1)  and  (22)  are  due  to  the  jet  flow  through  the 
highly  flexible  envelope.  The  third  and  final  contribution  to 
the  bending  moment  is  a  well-known  result  of  classical  beam 
theory.  The  bending  moment  due  to  the  bending  stiffness  of 
the  annular  duct  is12 

d  2v 

Mb=EI±JL,  (23) 

dx- 

where  I  =  vPt. 

The  total  bending  moment  is  obtained  by  adding  Eqs. 
(11),  (22),  and  (23) 

M  =  [El  +  I0(lpV2  +  2 p0))  *pL .  (24) 

dx4 

Combining  (24),  (8),  and  (7),  we  obtain  the  equilibrium  condi¬ 
tion 

[El  +  /0(3 pV2  +  2 p0)]  +  mfv2  4 \ 


This  result  indicates  that  the  effective  bending  stiffness  of  the 
duct  is  increased  by  nonsymmetric  pressure  and  tension  dis¬ 
tributions  created  by  the  slight  local  curvature  of  the  duct. 

The  boundary  conditions  which  much  be  satisfied  at  the 
upper  (fixed)  and  lower  (free)  ends  of  the  duct  are 

y  =  0  at  x  —  0,  (26) 

-^-  =  0  at  x  =  0,  (27) 


[El  +  3 pV%) 


d'y  = 


m,V2S  at  x  —  L. 


Equations  (26H28)  are  the  standard  end  conditions  applied 
to  the  fixed  and  free  ends  of  a  simple  cantilever  and  require 
no  further  comment.  Equation  i29|  ts  unique  to  the  problem 

R  Anderson  ana  A  Seian  3 '95 


FIG.  4.  Origin  of  the  transverse  force  Fs.  The  magnitude  of  5  has  been 
greatly  exaggerated. 

of  a  highly  flexible  duct  conveying  fluid,  and  incorporates 
the  important  physical  observation  that  no  real  duct  is  ever 
perfectly  straight.  The  weight  of  the  fluid  flowing  through 
the  duct  will  tend  to  straighten  it  except  for  a  short  region  of 
unknown  length  ds  near  its  free  end  (Fig.  4),  This  results  in  a 
small  transverse  force  Fs  due  to  the  reaction  force  of  the 
fluid  as  it  leaves  the  duct: 

Fs  =  mfV2  sin[ds/R0}xm/V2  dsl^oL  ~l-  (30) 

This  reaction  force  is  expressed  in  (29)  in  terms  of  the  un¬ 
known  constant 

S-d*yg)amL.  (31) 

It  is  shown  below  that  the  sole  function  of  5  is  to  determine 
the  amplitude  of  the  buckled  duct. 

Dividing  (25)  by  El  +  3pF2/0and  introducing  the  non- 
dimensional  variables 

y^y/L,  (32) 

x  =  x/L,  (33) 

reduces  the  equilibrium  equation  and  boundary  conditions 
to  the  form 

[i  +  rti-*)]4£ 

dx 

-~(pz[\+a0{x-\)}^=O,  (34) 

y  =  0  at.x  =  0,  (35) 

0  ati  =  0,  (36) 

dx 


q-o, 

=  at*=l.  (38) 

The  centerline  pressure  in  the  duct,  p0,  has  been  ex¬ 
pressed  as  a  linear  function  of  x  by  relating  the  pressure 
gradient  to  the  dynamic  head  via  a  resistance  coefficient  co: 


dp  _  (0  (  L_ 


For  turbulent  flow  the  resistance  coefficient  is  small  and  ap¬ 
proximately  constant 

w  =  D(f).  (41) 

The  nondimensional  parameters  appearing  in  (34)  are 

*• 

t oL/D  , , ,, 

y  — - - - .  44 

EI/pV%  +  3 

Each  of  these  parameters  depends  upon  F2,  reflecting  the 
strong  influence  of  fluid  velocity  shown  in  Fig.  1. 

The  final  character  of  the  solution  to  (32)— (36)  depends 
to  a  large  degree  upon  the  sign  of  the  function 

z—  1  +aP(x-  1),  (45) 

which  appears  in  the  second  term  of  Eq.  (34).  Function  z 
changes  sign  when  x  passes  through  the  “turning  point” 

xTP  =  1  —  1  /a0.  (46) 

The  location  of  this  point  defines  the  transition  between  the 
buckled  and  unbuckled  portion  of  the  duct  [Fig.  1(a)],  and  is 
positive  and  nonzero  provided  that  a/3>  1.  This  can  be  seen 
explicitly  by  considering  the  velocity-dependent  behavior  of 
a  short,  flexible  duct  with 

L  /D  =  0(1),  (47) 

El  „  (  (length/time)2  ^ 

-  j,  (48) 

gD  ...  (length/time)2  j  (49) 

(mf  +  m)/m/  =  0(1).  (50) 

Conditions  (48)  and  (49)  are  intended  to  simulate  experimen¬ 
tal  conditions  when  the  physical  properties  of  the  duct  are 
held  constant  while  the  fluid  velocity  is  varied.  Under  these 
conditions,  a,  P,  and  y  become 

tt.o[ "  (-^+3n  (5" 

l  pF2I0  l  pF%  J  J 

0  =  O[4/lEI/pV2Io  +  S)],  (52) 

y  =  0[e/(EI/pF2I0  +  3)].  (53) 

The  velocity  dependence  of  these  parameters  is  summarized 
in  Table  I. 

To  demonstrate  the  general  behavior  of  the  solution  we 
will  first  consider  the  small  velocity  limit pV2I0/EI  =  O  (e). 
By  neglecting/? 2  and  y,  integrating  once  and  introducing  the 
change  of  variables 


TABLE  I.  Velocity  dependence  of  a,  0,  and  y  for  the  case  of  a  short,  flexible 
duct. 
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$  =  -a'r-0(x-  1),  (54) 

we  can  transform  (34)  into  the  homogeneous  Airy  equation 
^7-0  =  0.  (56) 

with  the  general  solution 

j>  =  C,  Ai,f)  +  C:Bi(£).  (57) 

However,  application  of  conditions  (36)  and  (37)  produces 
the  trivial  result  C,  =  C:  =  0  and  the  conclusion  that  the 
duct  will  not  buckle  when  the  velocity  is  small. 

In  the  large  velocity  limit  \pV1I0/EI  =  O  (1/e)]  we  pro¬ 
ceed  in  a  similar  manner.  Neglecting  a  and  y  and  introduc¬ 
ing  the  change  of  variables 

_  dy 

y  =  ~T’  (58) 

dx 

Eq.  (34)  becomes 

4 5-+0Ij=/?J*(n  (59) 

dx~ 

K{V)~S+y{l).  (60) 

This  is  equivalent  to  the  horizontal  problem  considered  by 
Lundgren  et  al.9  and  leads  to  the  solution 

y  =  5[cos/#  —  cas(3(x  —  1)],  (61) 

and  the  conclusion  that  for  large  velocities  the  entire  length 
of  the  duct  is  buckled.  Note  that  the  unknown  parameter  5 
serves  only  to  determine  the  amplitude  of  the  solution  and 
does  not  influence  the  buckled  shape. 

Finally,  in  the  moderate  velocity  limit  \pV‘Io/ 
EJ  =  O  ( 1 )]  we  neglect  y,  integrate  once  and  introduce  the 
change  of  variables 

V~  -z/a2li,  (62) 

>'-%■  1631 

These  steps  reduce  (34)  to  the  form  of  a  nonhomogeneous 
Airy  equation 

dr, 2  Vy  a/3  ’  (  ) 

with  the  solution13 

..«  _  ~  S  t  Oi(?)  +  C ,  Ai(7y)  +  C2  Bi(g)  ] 

a/3  [Gi'(b )  +  C,  Ai"(b )  +  C2  Bi"(h )] 

C  _  [Gi'(6  )Bi(a)  —  Gi(a)Bi'(6 )] 

1  [  Ai(a)Bi'(6 )  —  Ai'(6  )Bi(a)]  ’ 
c  =  [Gi(a)Ai'(6)  —  Gi’(6  )Ai(a)] 

2  [  Ai(a)Bi'(6 )  —  Ai'(6  )Bi(a)J 

Gil?;)  =  i-Bil^)  +  ^  fW  V*  -  —  fBiU  )dt,  ,68) 
3  dr]  Jo  dr,  J n 

a  =  -(l  ~a0)/a2,\  (69) 

b  =  -  1/a2'3.  (70) 

The  functions  Ai  and  Bi  are  shown  in  Fig.  5.  The  solu¬ 
tion  165)  will  exhibit  buckling  behavior  when  tj  is  negative 


FIG.  5.  Plot  of  the  Airy  functions  Ai|x|  and  Biui  (Taken  from  Ref.  14,  p 
69.) 

and  exponential  behavior  when  r,  is  positive.  The  transition 
between  the  two  forms  of  behavior  is  defined  by  the  location 
of  the  turning  point  expressed  in  (46). 

III.  LONG  DUCT  EXPERIMENTS 

In  the  previous  section  we  have  shown  that  the  qualita¬ 
tive  behavior  of  the  solution  to  (34)-(38)  agrees  with  the  gen¬ 
eral  observations  described  in  the  introduction.  To  test  this 
agreement  in  quantitative  terms  we  experimented  with  long 
ducts  in  the  fully  buckled  regime  (a0<  1).  After  testing  sev¬ 
eral  different  materials,  we  chose  to  do  our  measurements 
using  surgical  drainage  tubing  made  of  latex  rubber.  We 
used  two  different  tube  sizes:  the  smaller  had  an  inner  diame¬ 
ter  of  9.5  mm  and  a  wall  thickness  of  0.3  mm,  while  the 
corresponding  dimensions  of  the  larger  size  were  12.7  mm 
and  0.34  mm.  For  both  sizes,  the  length  of  the  tube  was  one 
meter.  The  flexible  tube  was  clamped  to  an  inlet  nozzle  con¬ 
nected  to  the  building  water  supply.  The  Reynolds  number 
based  on  duct  diameter  varied  from  104  to  5  X 10*. 

Figure  6  shows  a  schematic  drawing  of  the  experimen¬ 
tal  apparatus.  In  order  to  avoid  the  secondary  effects  asso- 


inlet  nozzle 


loudspeaker 


steel  l  beomi 


rubber  tube 


'load  cell 


FIG.  6,  Schematic  of  enperimertal  apparatus 
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ciated  with  bends  (elbows)  in  the  piping,  the  length  of  the 
inlet  section  was  always  greater  than  50  times  the  rubber 
tube  diameter.  The  inlet  section  had  the  same  inner  diameter 
as  the  rubber  tube,  and  the  inlet  nozzle  was  tapered  to  a  knife 
edge  in  order  to  avoid  flow  separation  at  the  point  where  the 
water  leaves  the  nozzle  and  enters  the  rubber  tube.  The  flow 
was  measured  by  taking  timed  volumetric  samples  of  the 
water  leaving  the  flexible  tube.  At  each  flow  setting,  the 
shape  of  the  buckled  tube  was  recorded  photographically 
(see  Fig.  1). 

As  argued  in  Sec.  II  of  this  paper,  the  bending  stiffness 
of  the  flexible  duct  is  a  pivotal  parameter  in  predicting  the 
buckled  shape  of  the  duct.  The  bending  stiffness  of  the  tubing 
material  was  measured  in  two  different  ways.  In  the  first 
method,  Young's  modulus  was  calculated  using  a  load  cell, 
by  measuring  the  force  required  to  elongate  the  tubing  by  a 
given  amount.  The  experimentally  measured  Young’s  mo¬ 
dulus  was  then  multiplied  by  the  area  moment  of  inertia 
calculated  from  the  physical  dimensions  of  the  tubing.  This 
value  was  double  checked  by  measuring  the  fundamental 
frequency  of  a  short  section  of  the  water-filled  tube.  The 
corresponding  value  for  El  was  then  calculated  by  using  the 
formula  for  a  simple  cantilever  beam.  The  frequency  mea¬ 
surement  was  accomplished  by  shaking  the  tube  with  an  arm 
attached  to  the  cone  of  a  loudspeaker  (see  Fig.  6).  The 
Young’s  modulus  determined  by  these  methods  was  found 
to  be  1 .7  X  106  N/m2  with  the  static  test  giving  values  10%- 
15%  higher  than  those  determined  by  the  dynamic  test. 

A  typical  experimental  run  was  started  with  the  tube 
hanging  limp  (collapsed,  with  elliptical  cross  section).  When 
the  water  was  turned  on  the  duct  began  to  buckle  from  its 
free  end  as  described  in  the  introduction.  The  length  of  the 
buckled  portion  of  the  duct  increased  with  increasing  flow 
velocity  until  a0  reached  the  value  0. 3  5  for  the  larger  diame¬ 
ter  duct  (0.43  for  the  smaller  diameter  duct)  when  self-excit¬ 
ed  vibrations  began  near  the  free  end  of  the  duct.  The  ducts 
were  rotated  end  for  end  during  the  experiment  to  test  for 
systematic  error  associated  with  the  specific  shape  of  the 
duct.  The  only  quantity  that  was  found  to  vary  significantly 
was  the  amplitude  of  the  buckled  shape,  in  agreement  with 
results  (61)  and  (65). 

The  scaling  that  is  appropriate  for  the  ducts  used  during 
the  experiment  differs  from  the  short  ducts  considered  in  the 
previous  section: 


gD/V2  =  0(e),  (71 1 

EI/pV2Iu  =  0{ \/e),  (72) 

L/D  =  0{  1/f),  (73) 

\mf  +  m)/m/  =  0(1),  (74) 

a  =  0(euz),  (75) 

0  =  O(e-'12),  (76) 

Y  =  0  (e).  (77i 

Neglecting  y,  integrating  once,  and  introducing  the  change 
of  variables 

u=‘Q-S^ L,  (78) 

dz  K[V) 

reduces  (32)  to  the  form 

a2^-  +  zu  =  1.  (79) 

dr 

Functions  K  and  z  have  been  previously  defined  in  (60)  and 
(45).  Because  we  expect  a  spatially  periodic  result,  we  will 
seek  a  solution  in  terms  of  an  exponential  WKB  expansion 14 


u~exp^2an*"(2l)- 


The  functions  s„  are  determined  by  substituting  this  expan¬ 
sion  into  the  homogeneous  portion  of  (79) 

*0=  (81) 

r,  =  -  Jlnz.  (82) 

When  expressions  (81)  and  (82)  are  introduced  back  into  Eq. 
(80)  we  find  the  WKB  solution  for  the  homogeneous  equa¬ 
tion  to  be 

“*  -  ^  [c'  c“(  +  c-  ”"( '”’)]■  1831 

The  corresponding  particular  solution  is  found  by  solving 
(79)  in  the  limit  a2  =  0  and  is 

up  =  l/z.  (84) 

Combining  (83)  and  (84),  applying  conditions  (35)— (38)  and 

transforming  back  to  the  nondimensional  physical  coordi¬ 

nate  Jc,  leads  to  the  final  form 


§ " So  I T+^  lc'  “K  to1”’) +  C;  sin( 

£.  _  (l/C3/4)[sin(2/3a)/4  —  cos(2/3g|/q)  —  sin(2C  3;2/3a) 
(cos(2/3a)(C3'2  —  i)/a  +  sin(2/3a)(C3/2  —  11/4] 

C  _  cos(2C  V2/3a)  —  (l/C3/4)[sin(2/3a)/a  +  cos(2/3a)/4] 

J  (cos(2/3a)(CJ/2  -  l)/a  +  sin(2/3a)(C 3,2  -  l)/4]  ’ 

C  =  1  -  cl0. 
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FIG.  7.  Plot  of  An  when  V  -  4.05  m/sec.  a  —  0.01, /?  =  57.4,  Z)  =  12.7  mm, 
L/D  —  72,  £7  =  4.65  x  10~4  Nnr,  and  -r  ml  =  0.895.  (  +  =  ex¬ 

perimental  value,  O  =  theoretical  value.  | 


FIG.  9.  Plot  of  A0  when  V  =  4.56  m/sec.  a  =  0.005.  0  =  75.2.  All  other 
parameters  have  the  same  value  as  Fig.  7.  |  —  =  experimental  value, 
O  =  theoretical  value  ! 


Equation  (85)  was  integrated  numerically  from*  =  0  to 
x  =  1.  During  the  integration  the  distance  between  succes¬ 
sive  crossings  of  an  imaginary  line  extending  from  the  origin 
to  the  free  end  of  the  duct  was  recorded  for  comparison  with 
measurements  taken  from  photographs  of  buckled  ducts. 
We  define  the  distance  between  crossings  to  be  A.0  and  show 
the  results  of  this  comparison  in  Figs.  7-9.  The  experimental 
measurements  are  denoted  by  +  ’s  and  the  numerical  results 
derived  from  ( 85 )  are  denoted  by  O’s.  The  agreement  is  good, 
both  in  terms  of  the  location  of  the  crossings  and  the  magni¬ 
tude  of  A0.  The  total  number  of  crossings  measured  from 
photographs  agree  with  those  predicted  analytically,  except 
for  the  case  shown  in  Fig.  9.  It  is  felt  that  this  is  a  result  of  our 
inability  to  interpret  the  photographs  accurately  near  the 
free  end  of  the  duct. 

IV.  CONCLUSIONS 

We  have  shown  analytically  that  static  buckling  of  a 
highly  flexible  duct  filled  with  flowing  fluid  is  a  result  of  the 
shear  force  created  by  the  unavoidable  initial  curvature  of 
the  duct.  We  found  that  this  shear  force  governs  the  ampli- 


o 
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FIG  8.  Plot  of  A„  when  V  =  4  43  m/sec.  a  =  0  008,  0  =  62.0  All  other 
parameters  have  the  same  value  as  Fig  7  (  +  =  experimental  value. 
O  =  theoretical  value  l 


tude  of  the  buckled  duct,  but  does  not  afTect  its  general 
shape.  In  addition,  we  found  that  a  vertical  duct  exhibits  a 
transition  between  sinuous  and  exponential  deformation. 
The  location  of  this  transition  was  found  to  depend  upon  the 
velocity  of  the  fluid  stream.  We  successfully  tested  our  anal¬ 
ysis  by  measuring  the  buckling  wavelength  of  a  long  duct 
and  comparing  it  with  the  predicted  values. 

It  is  important  to  note  that  the  buckling  property  of  a 
bounded  jet,  which  has  been  explored  in  this  paper,  is  in 
principle  independent  of  the  bending  stiffness  of  the  duct 
which  surrounds  the  jet,  because  of  the  inherent  jet  bending 
stiffness  defined  by  Eq.  (11).  In  the  present  paper  we  have 
examined  the  general  case  where  the  duct  walls  have  finite 
bending  stiffness.  However,  a  similar  buckling  phenomenon 
is  observed  in  jets  bounded  by  ducts  with  zero  bending  stiff¬ 
ness,  for  example,  in  fast  capillary  jets15  and  in  inviscid  jets 
that  discharge  freely  into  an  infinite  reservoir. 16  In  the  latter 
case,  the  hydrodynamic  stability  analysis  demonstrates  that 
the  jet  moves  as  if  it  were  bounded  by  the  annular  shear  layer 
which  surrounds  it,  and  selects  a  longitudinal  wavelength 
that  is  a  certain  multiple  of  the  jet  diameter."’  This  is  consis¬ 
tent  with  the  buckling  wavelength  produced  by  the  EI—J) 
limit  of  the  present  analysis. 1 1  A  compilation  of  photographs 
of  buckling  natural  flows  is  presented  in  Ref.  1 7. 

From  the  point  of  view  of  experimental  methods  in  tur¬ 
bulence  research,  the  present  flexible-duct  experiments  illus¬ 
trate  an  effective  way  to  visualize  the  natural  tendency  of  any 
turbulent  duct  flow  to  engage  in  large-scale,  organized,  mo¬ 
tion.  As  the  beam  stiffness  of  the  hose  decreases,  the  hose 
wall  molds  itself  to  the  large-scale  sinuous  path  preferred  by 
the  turbulent  jet,  and  the  hose  shape  becomes  a  record  of  the 
large-scale  structure  of  turbulent  jet  flow.  As  pointed  out  in 
the  Introduction,  the  search  for  an  effective  jet  flow  visual¬ 
ization  method  is  what  led  to  the  definition  of  the  problem 
addressed  in  this  paper. 
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The  Buckling  of  a  Vertical  Liquid 
Column 

This  paper  documents  the  sinuous  flow  of  a  t  ertical  water  column  uhn  unprrees 
on  a  horizontal  surface.  It  is  shown  experimentally  that  the  solid  ohstaste  hut  kies 
the  column  into  a  plane  sinuous  ( static )  shape  whose  wave/emtih  'tales  with  the 
local  column  diameter.  The  measured  ratio  {buckling  wavelength)  {column 
diameter)  is  of  order  it/ 2,  in  agreement  with  (he  buckling  theory  of  tmiscid  flow. 


1  Introduction 

The  object  of  this  paper  is  to  report  a  series  of  experimental 
observations  concerning  the  buckling  of  a  vertical  water 
column,  which  is  a  phenomenon  commonly  observed  by 
placing  a  flat  obstacle  under  a  faucet.  The  observations  reveal 
the  existence  of  an  interesting  scaling  law,  namely,  the 
proportionality  between  the  diameter  of  the  water  column  and 
its  longitudinal  buckling  wavelength. 

The  flow  of  vertical  liquid  columns  is  a  topic  which  draws 
its  importance  from  its  many  engineering  applications,  for 
example,  gas-to-liquid  chemical  processing  and  condensation 
heat  transfer  in  heat  exchanger  design.  As  shown  below,  this 
topic  has  been  studied  extensively,  however,  with  only  one 
exception,  all  previous  studies  have  been  focused  on  the  free- 
fall  of  the  liquid  column:  in  the  free-fall  regime  the  column 
remains  axisymmetric.  In  the  present  study  we  consider  a 
vertical  column  which  is  interrupted  by  an  obstacle  and,  as  a 
result,  assumes  a  sinuous  (nonaxisymmetric)  shape. 

The  two-dimensional  inviscid  flow  issuing  from  a  very  large 
vessel  through  an  aperture  in  a  solid  wall  is  one  of  the  classical 
problems  in  hydrodynamics,  pioneered  by  Kirchhoff  [1]  and 
Rayleigh  (see  Lamb  (2|)  Milne-Thomson  (3[  considered  the 
same  phenomenon  under  the  effect  of  gravity.  Since,  under 
gravity,  the  "infinite  height  of  the  vessel  leads  to  an  “infinite 
hydrostatic  pressure”  paradox,  Conway  [4]  redid  the  problem 
by  inserting  a  horizontal  lid  above  the  wall  with  the  orifice. 
He  was  able  to  obtain  the  smooth  draw-down  shape  of  the 
column,  employing  the  free  surface  treatment  proposed  by 
Milne-Thomson  13], 

Among  the  numerous  theoretical  studies  reported  in  the 
literature,  the  work  of  Scriven  and  Pigford  (5]  appears  to  be 
the  first  in  which  the  liquid  flow  was  analyzed  based  on  a 
boundary  layer-type  approximation.  Lienhard  [6]  chose  a 
similar  approach  and  included  also  the  effects  of  gravity  and 
surface  tension.  Matovich  and  Pearson  (7]  formulated  a  one¬ 
dimensional  flow  model  and  were  able  to  obtain  several 
closed-form  solutions  for  symmetric  column  shapes  in  various 
parametric  domains.  A  more  rigorous  formulation  was 
presented  by  Duda  and  Vrcntas  (h).  A  one-dimensional 
Cosserat-tvpe  theory  was  developed  by  Green  |4)  and  Bogy 
[I0|.  The  work  of  Bogy  is  concerned  primarily  with  the 

(.  onirihuted  h\  the  Hunts  tnameefinji  Division  lor  puhlnjnoit  in  ihc 
Ioimsal  or  t  n  ids  1-suism  him.  M.it'ii'cnpr  rcvcivcJ  Pv  ihc  !  luiJ- 
l  n*fiiK'criiiki  Dmnon.  (Kriohcf 


condition  of  droplet  formation  at  the  lower  end  of  a  tree¬ 
falling  column.  The  subject  of  axisymmetric  break-up  and 
droplet  formation  has  been  studied  extensively  over  the  past 
100  years;  the  core  of  this  work  was  reviewed  recently  in  the 
context  of  Hydrodvnamic  Stabilit v  Theorv  (Drazin  and  Reid 

(HI). 

To  our  knowledge,  Lienhard  [12]  is  the  only  author  to  have 
analyzed  the  sinuous  shape  of  a  column  impinging  on  an 
obstacle.  Lienhard  viewed  the  phenomenon  as  capillary  waves 
which  originate  at  the  solid  plane  and  travel  upward  at  the 
same  speed  as  the  water  velocity.  It  is  interesting  to  note  that 
the  buckled  shape  of  vertical  water  columns  is  sometimes 
accidentally  photographed  and  published  by  researchers 
investigating  an  unrelated  phenomenon  (see,  for  example, 
Figs.  1(b)  and  1(c)  in  Winston  and  Martin  [13]). 

2  Experiment 

We  carried  out  a  series  of  controlled  experiments  in  which 
we  photographed  the  buckled  shape  of  the  liquid  column.  The 
purpose  of  this  photographic  record  is  to  document  the 
relationship  between  the  geometry  of  the  buckled  column  and 
physical  parameters  of  the  flow  (velocity,  column  diameter, 
column  height). 

The  experimental  setup  is  shown  schematically  in  Fig.  1. 
The  apparatus  consists  of  three  basic  components,  the  water 
supply  reservoir,  the  nozzle,  and  the  round  table  which  in¬ 
terferes  with  the  free-fall  of  the  water  column.  The  water 
reservoir  is  a  cylindrical  vessel  (14  cm  diameter,  60  cm  height) 
which  supplies  water  at  a  known  flowrate  to  the  nozzle  sec¬ 
tion.  The  flowrate  is  known  with  5  percent  accuracy  from 
preliminary  calibration  experiments  which  revealed  the 
relationship  between  flowrate  and  reservoir  water  level.  The 
reservoir  and  the  nozzle  are  connected  through  a  1  m  long 
vinyl  tube  with  an  inside  diameter  of  0.64  cm.  The  nozzle 
holder  was  designed  to  accommodate  four  different  nozzles 
which  produce  water  columns  of  varying  diameters.  \s  show  n 
in  the  detail  drawing  of  Fig.  I .  in  each  noz/ie  the  flow  is  !  ir-t 
lead  through  a  wire-mesh  section  in  order  to  damp  v'ut 
potential  disturbances  originating  from  upstream.  Water 
columns  of  different  heights  ansi  diameters  were  established 
between  the  nozzle  and  the  round  table.  Note  that  the  table 
surface  has  a  slight  spherical  curvature  in  order  to  prevent  the 
irregular  accumulation  of  water  right  under  the  column. 
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Fig.  1  Experimental  apparatus 

The  parametric  domain  covered  by  the  experiments  is 
described  in  terms  of  two  dimensionless  groups, 


(1) 


(2) 


Fig.  2(a)  B  =  0.34,  W  =  1.4,  z/r0  =8.8 


where  p,  g,  u,,,  r.u  and  T  represent  liquid  density, 
gravitational  acceleration,  liquid  velocity  at  the  nozzle  exit, 
nozzle  radius,  and  water-air  surface  tension,  respectively.  We 
varied  the  Weber  number  by  changing  the  flow  rate:  although 
the  Weber  number  attained  in  the  experiment  depends  on 
nozzle  size,  it  falls  consistently  in  the  range  of  0.1-4.  The 
Bond  number  was  varied  from  0.085  to  1 .37  by  changing  the 
nozzle  size  (the  four  nozzle  diameters  were  0.159  cm,  0.138 
cm,  0.476  cm  and  0.635  cm).  The  Reynolds  number  based  on 
conditions  at  the  nozzle  exit,  2r„i/,,,  t\  varied  from  450  to 
1 100.  The  velocity  profile  in  the  water  jet  immediately  below 
the  nozzle  was  measured  recently  by  Yamaguchi  and 
Takahashi  [17]:  their  Laser-Doppler  measurements  showed 
that  the  velocity  profile  becomes  flat  within  a  longitudinal 
distance  of  order  r0  downstream  from  the  nozzle  exit. 
Therefore,  in  discussing  the  present  experiments  it  is  assumed 
that  the  velocity  is  approximately  uniform  in  each  jet  cross- 
section. 

Each  experimental  run  was  conducted  as  follows:  first,  we 
fixed  the  two  parameters  (H’.fl)  by  selecting  a  certain  nozzle 
and  by  keeping  the  flowrate  constant  (note  that  the  water 
reservoir  level  could  be  maintained  by  replenishing  the  water 
supply  during  the  run).  The  table  directly  under  the  nozzle  was 
then  lowered  until  a  new  (additional)  elbow  was  observed  in 
the  buckled  shape  of  the  column.  The  height  of  the  table  and 
the  counted  number  of  elbows  were  recorded  and  a 
photograph  of  the  column  shape  was  taken.  The  procedure 
was  repeated  for  various  Weber  and  Bond  numbers. 

Representative  photographs  of  the  buckled  liquid  column 
are  shown  in  Fig.  2.  The  experimental  observations  are 
summarized  in  Fig.  4,  as  the  total  number  of  elbows  vs.  the 


Fig.  2(6)  B  =  0.34,  W  =  3.65,  zlr0  =7.68 

Fig.  2  Photographs  showing  the  sinuous  shape  ot  buckled  liquid 
columns 


dimensionless  height  of  the  column  z  =  c*/r„.  We  learn  that  in 
all  cases  the  number  of  elbows  increases  almost  linearly  with 
the  height  of  the  column.  This  conclusion  is  of  fundamental 
importance,  because  it  points  towards  the  existence  of  a 
universal  proportionality  between  column  diameter  and 
buckling  wavelength. 


B  =  Bond  number  given  by 
(Pgr,,'-)/  T 

D  -  water  column  diameter 
g  =  gravitational  acceleration 
r  =  nozzle  radius 
Re  =  Reynolds  number,  r 
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T  =  water-air  surface  tension 
u, i  =  liquid  velocity  at  the  nozzle 
evil 

U  =  Weber  number  given  by 
((>!/,, ' r  .)  7 

;  =  hemht  ot  ihe  water  column 


\  =  buckling  wavelength 
u  =  density 
r  =  kinematic  visgositv 

Superscript 

'  dimensional  quantity 

Transactions  ol  the  ASME 


;  .-.V 

ft  -J-  . 


-•ft. _ r  .  -ZW.  -ft  . 


V*  Ajt  v  T_l« 


*3 

.  V 


n 


L4 


zv7-_*v  ; 


.»*«?!.... 


■J»  J  UK  till  MU.  t-UJf.1,'.  I  I.N.'iBH.I.I.  '  I.  I.  I  I  — W  I  1 .1  J  ,4'WIP.  ».  t.  ft 


Fig.  3(a)  8  =  0.34.  W  =  0.59,  z/Fq  =  9.4 


Fig.  3(b)  8  =  0.77,  W  =  1.77,  Zlr0  =  7.5 


3  The  Two-Dimensionality  ol'  the  Buckled  Shape 

An  important  question  arises  in  connection  with  the 
buckled  shape  of  the  liquid  column:  is  this  shape  two- 
dimensional,  i.e..  in  one  plane,  as  the  buckled  shape  of  an 
elastic  rod,  or  is  it  helical?  In  the  course  of  photographing  the 
columns  displayed  as  Figs.  2 (a,  b ),  the  liquid  column 
exhibited  the  tendency  to  buckle  in  a  way  which  made  it 
visible  only  from  a  special  direction.  Indeed,  the  camera 
which  recorded  Figs.  2 [a.  b)  had  to  be  positioned  so  that  the 
buckled  shape  appeared  most  visibly  on  film.  However,  in 
addition  to  plane  buckling,  tall  columns  showed  a  tendency  to 
develop  an  axisymmetric  (varicose)  deformation  at  the  lower 
end.  near  the  round  table  serving  as  obstacle.  It  is  important 
to  establish  whether  this  axisymmetric  deformation  is  an 
integral  part  of  the  plane  buckling  phenomenon  or  simplv  a 
reflection  ot  the  geometnc  symmetry  imposed  by  the  obstacle 
on  the  liquid  column. 

In  order  to  answer  these  questions,  we  conducted  an  ad¬ 
ditional  series  of  experimental  observations  leported  here  via 
figs.  3u-3c.  The  buckled  column  was  photographed 
simultaneously  irom  the  front  and  from  the  side,  bv  using  a 
vertical  mirror  whose  plane  made  a  45  deg  angle  with  the 
camera-column  direction,  in  lugs.  3</-.V  the  direct  (trout) 
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Fig.  3(c)  8  =  0.77.  W  =  1.77,  z/r0  =  6.5 

Fig.  3  The  plane  buckled  shape  ot  the  liquid  column.  Left  side:  the 
direct  view,  right  side:  the  side  view  through  the  mirror 


view  appears  on  the  left  side,  while  the  mirror  (sidei  view 
appears  on  the  right  side  of  the  photographs.  To  focus  on  the 
column  and  its  mirror  image  at  the  same  time,  the  mirror  had 
to  be  positioned  no  farther  than  1  cm  away  from  the  liquid 
column.  As  the  obstacle  we  used  an  asymmetric  object, 
namely,  the  flat  end  of  a  screwdriver.  The  flat  end  was 
positioned  so  that  is  made  approximately  a  10-degree  angle 
with  the  horizontal  plane  (see  Figs.  3<z-3c). 

Figures  3a-3c  show  very  clearly  that  the  liquid  column 
prefers  to  buckle  in  a  plane.  The  view  from  the  direction 
normal  to  this  preferred  plane  reveals  a  nearly  straight 
column. 

In  conclusion,  the  symmetric  (varicose)  or  asymetric 
(sinuous)  deformation  of  the  liquid  column  depends  on  the 
degree  of  asymmetry  of  the  obstacle.  If  the  obstacle  is  from 
the  side,  as  in  Figs.  3o- 3c,  then  the  column  buckles  into  a 
plane  sinuous  shape.  If,  on  the  other  hand,  the  obstacle  is 
symmetric  (i.e.,  the  same  when  seen  from  all  directions)  then 
the  deformation  tends  to  be  symmetric  also. 

4  Discussion  of  Results 

The  total  number  of  elbows  measured  experimentally  is 
shown  in  Figs.  4u-4rf.  Each  figure  shows  the  total  number  of 
elbows  as  a  function  of  column  height  ;  (dimensionless),  for 
various  values  of  B  and  W.  In  each  experiment  the  Bond 
number  is  fixed;  it  was  observed  that  below  a  critical  Weber 
number  the  vertical  column  breaks  up  into  droplets  right 
under  the  nozzle,  and  that  the  buckled  shape  is  most  v  isible  in 
the  range  where  W  is  slightly  higher  than  this  critical  value. 
The  experimental  observations  summarized  in  Figs.  4a-4d 
have  been  performed  in  this  optimum  range  of  Weber 
numbers. 

The  measurements  indicate  that  the  number  of  elbows  of 
the  buckled  column  increases  almost  linearly  with  the  height 
of  the  column,  thus  suggesting  a  local  proportionality  bet¬ 
ween  buckling  wavelength  \  and  column  diameter  D  The 
Weber  number  appears  to  have  only  a  weak  influence  on  the 
measurements  of  Figures  4a-4d.  The  Bond  number  ettect  is 
more  noticeable,  as  the  number  of  elbows  lor  c  =  ss>nstam 
increases  as  the  Bond  number  I.  ncreases:  This  dependence  is 
illustrated  in  Fig.  5. 

I  he  proportionality  between  X  and  D  is  illustrated  turner 
in  lags,  b  and  '.  In  general,  both  P  and  a  dectease  as  the 
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Fig.  4  Total  number  ot  elbows  as  a  (unction  of  column  height 


water  column  is  accelerated  downward  (Figs.  2,  3)  however, 
the  measured  ratio  \/D  appears  to  be  independent  ot  the 
vertical  position  along  the  column  (Fig.  6).  For  this  reason,  in 
Figure  7  we  show  only  the  ,VD  measurements  taken  from  the 
upper  section  of  each  column,  where  the  first  buckling  takes 
place.  The  \/D  ratio  is  largely  insensitive  to  changes  in  both 
Weber  number  and  Bond  number.  Indeed,  in  the  1V-B  range 
documented  by  our  experiments,  the  wavelength/diameter 
ratio  appears  to  be  a  universal  constant  of  order  1 .5.  Figure  7 
'hows  also  Lienhard’s  [12|  theoretical  curve  and"  \/D 
measurements  tor  the  axisymmetric  (varicose)  deformation: 
mite  that  Lienhard  s  data  merge  with  the  present 
measurements  in  the  vicinity  of  «•=()(  1),  and  that  Lienhard’s 
theoretical  wavelength  ot  the  varicose  shape  appears  to  be 
consistent  with  both  sets  ot  data.  However,  Lienhard’s  theory 
does  not  account  for  the  sinuous  (nonaxisymmetric)  shape 
documented  by  the  present  experiments,  nor  does  it  account 
lor  the  apparent  constancy  of  the  V D  ratio  of  buckled 
(sinuous)  shapes. 


We  feel  that  the  universal  proportionality  between  buckling 
wavelength  and  column  diameter  is  of  fundamental  interest 
for  the  following  reasons.  First,  a  similar  \~D  propor¬ 
tionality  is  also  exhibited  by  highly  viscous  fluid  filaments, 
such  as  honey  and  corn  syrup  falling  and  folding  on  a 
horizontal  surface  (see  the  photographs  of  Cruickshank  (14) 
and  Cruickshank  and  Munson  [15],  where  the  folding 
wavelength  is  greater  when  the  filament  is  thicker;  the  same 
scaling  is  revealed  while  experimenting  with  honey  filaments 
ol  increasing  thickness).  The  present  experiments  demonstrate 
that  the  buckling  phenomenon  and  the  \  ~D  scaling  law  are 
not  restricted  to  highly  viscous  columns,  as  originailv 
thought,  rather,  they  are  observed  even  when  "inviscid"  iluid 
columns  impinge  on  an  obstacle.  The  buckling  of  inviscid 
columns  was  linked  theoretically  to  the  origin  of  turbulence 
[lb|.  Thus,  the  present  measurements  confirm  not  onlv  the 
possibility  of  “buckling”  in  inviscid  jets,  predicted 
theoretically  [16],  but  also  the  accuracy  of  ihe  theoretK.il 
buckling  wavelength,  X. D  -  t  2=1.5'  [16], 
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Fig.  5  The  effect  of  Bond  number  on  the  total  number  of  elbows  over  a 
fixed  height 


Finally,  the  agreement  between  the  buckling  theory  of 
reference  [16]  and  the  X/D  data  of  Fig.  7  is  thought- 
provoking,  because  the  theory  neglects  the  effect  of  surface 
tension,  while  the  M''  =  0(1)  data  correspond  to  cases  where  the 
surface  tension  effect  is  not  negligible.  At  this  stage,  the 
agreement  between  theory  and  experiment  suggests  only  that 
surfaces  tension  is  not  a  major  consideration  in  "buckling" 
(i.e.,  in  the  rotational  equilibrium  condition  of  a  nearly 
straight  column  (161).  Whether  or  not  this  suggestion  is 
correct  can  be  decided  by  including  surface  tension  in  the 
buckling  theory  (this  remains  to  be  done)  and  by  extending  the 
present  experiments  to  document  the  Weber  number  domain 
below  W  =  0(1).  Note  that  recent  experiments  with  fast 
capillary  jets  in  the  high  Weber  number  range  10:  -  10J  [18] 
support  the  buckling  theory  of  inviscid  jets  [16,  19], 
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i  Experiments  on  the  Buckling  of 
I  Thin  Fluid  Layers  Undergoing 
K.  R.  Blake  j  End-Compression 

A.  Bejan  This  paper  reports  a  series  of  experiments  concerning  the  buckling  of  a  slender  fluid 

I  layer  in  a  state  of  longitudinal  compression.  The  experiments  consist  of  floating  a 
Department  ot  Mechanical  Engineering,  jj  layer  of  highly  viscous  oil  on  a  pool  of  water  and,  manually,  compressing  the  layer 
University  of  Colorado.  fj  from  the  side.  Photographs  of  the  buckled  layer  show  conclusively  that  the  buckling 
Boulder,  Colo  80309  rj  wavelength  is  largely  insensitive  to  either  the  rate  of  compression  or  the  viscosity  of 
i  the  fluid  layer.  The  observations  suggest  that  the  buckling  wavelength  is  actually  a 
§  characteristic  length  scale  (a  property)  of  the  fluid  layer,  in  contrast  with  the 
i  buckling  theory  of  purely  viscous  layers  ( Buckmaster,  Nachman,  and  Ting.  [  7  ] ) 
I  where  the  buckling  wavelength  remains  to  be  determined  randomly  by  initial 
§  disturbances. 


1  Introduction 

The  concept  of  “fluid  buckling”  is  relatively  new  in  fluid 
mechanics  research:  during  the  past  decade  it  has  been  used 
with  increasing  frequency  to  account  theoretically  for  a 
number  of  flow  phenomena,  the  explanation  of  which  ap¬ 
pears  to  shed  new  light  on  the  theoretical  origins  of  turbulence 
(Cruickshank  [1],  Bejan  [2],  Munson  [3]).  In  the  present 
paper,  the  buckling  concept  is  used  to  describe  a  series  of 
experimental  observations  of  how  a  sheet  of  viscous  fluid 
wrinkles  as  it  is  compressed  from  one  end.  As  shown  in  Fig.  I , 
the  fluid  buckles  and  assumes  a  sinusoidal  shape  with 
characteristic  wavelength  X.  The  objective  of  this  experiment 
is  to  measure  the  buckling  wavelength  and  to  learn  how  this 
wavelength  is  influenced  by  the  geometry  (slenderness)  of  the 
fluid  layer,  the  lluid  properties  and  the  rate  of  compression. 

The  buckling  of  highly  viscous  fluids  was  first  studied  by 
Biot  [4],  In  a  series  of  papers,  Biot  developed  the  equations  to 
describe  the  buckling  of  a  multilayered  viscous  fluid.  The 
equations  were  solved  and  were  shown  to  be  unstable  when 
the  viscous  layers  were  subjected  to  an  arbitrary  finite  strain 
with  a  small  perturbation  superimposed  on  the  initial  state. 
Solutions  were  determined  numerically  and  were  found  to 
agree  with  solutions  given  by  the  theory  of  elasticity  and 
viscoelasticity  when  the  instability  was  of  a  significant 
magnitude.  Biot’s  work  was  motivated  by  applications  to  the 
problems  of  tectonic  folding  of  stratified  geological  struc¬ 
tures. 

The  instability  of  jets,  threads,  and  sheets  of  viscous  fluids 
was  studied  qualitatively  by  Taylor  (5).  In  his  paper,  Taylor 
argues  that  the  instability  created  by  the  compression  ol 
viscous  fluids  is  the  same  as  the  Euler  buckling  of  solid 
columns.  Of  the  many  experiments  he  conducts,  two  are 
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especially  relevant  to  the  work  presented  here.  In  one.  Taylor 
compresses  a  thread  of  an  extremely  viscous  fluid  floating  on 
mercury  and  compares  the  resultant  shape  with  the  shapes  of 
elastica  under  compression  as  calculated  by  Love  [6).  The 
other  experiment  was  designed  to  determine  when  a  sheet  of  a 
viscous  fluid  under  compression  would  become  unstable. 

With  Taylor’s  experimental  results  in  mind,  Buckmaster, 
Nachman,  and  Ting  [7]  considered  theoretically  the  buckling 
of  a  thin  viscous  layer  (the  viscida  problem).  The  problem  was 
that  of  a  two-dimensional  viscida  immersed  in  vacuum,  w  hose 
ends  are  moved  together  slowly  so  that  inertia  terms  could  be 
neglected.  Since  their  analysis  was  not  limited  to  small 
deformations,  it  was  a  generalization  of  Biot’s  work.  They 
derived  a  global  equation  for  the  evolution  of  the  slope  of  the 
centerline  as  a  function  of  time  and  distance  from  one  end,  by- 
integrating  the  momentum  equations  over  the  thickness  of  the 
viscous  layer.  This  equation  was  then  solved  using  asymptotic 
expansions  for  the  case  of  small  centerline  deviation.  Based 
on  these  results  the  solutions  for  large  centerline  deviations 
were  found  numerically.  The  case  where  surface  tension 
effects  cannot  be  neglected  was  incorporated  into  this  theory 
by  Buckmaster  and  Nachman  [8]  in  a  subsequent  paper. 

In  a  more  recent  paper,  Suleiman  and  Munson  (9)  in¬ 
vestigated  the  buckling  of  a  thin  viscous  fluid  layer  subjected 
to  linear  shear.  They  found  that  if  the  dimensionless  shear 
stress  exceeded  a  critical  value,  the  layer  would  buckle  in  a 
manner  similar  to  the  buckling  of  a  thin  elastic  plate.  In 
another  study,  Munson  [10]  examined  experimentally  the 
buckling  of  a  falling  viscous  jet  Bowing  out  of  a  vertical  slit 
orifice  onto  a  horizontal  plate.  In  this  experiment  Munson 
observed  that  the  jet  sometimes  buckles  in  a  manner  similar  to 
the  buckling  of  a  cantilever  beam.  Another  experimental 
study  was  conducted  by  Cruickshank  and  Munson  [1 1|  on  the 
spontaneous  oscillations  of  a  falling  viscous  jet  Bowing  from 
a  horizontal  orifice  onto  a  flat  plate,  flies  determined  the 
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Fig.  t  The  buckling  o(  a  thin  viscous  layer  floating  on  top  of  a  heavier 
inviscld  fluid 


minimum  distance  from  the  jet  orifice  to  the  flat  plate  for 
which  the  jet  would  buckle  (termed  the  “buckling  height”)  as 
a  function  of  fluid  and  flow  variables.  They  found  that  if  the 
flow  exceeded  a  critical  Reynolds  number,  the  jet  does  not 
buckle.  Below  this  critical  Reynolds  number,  surface  tension 
becomes  the  dominant  factor  in  influencing  the  buckling 
height.  The  origin  of  buckling  was  attributed  to  the  jet  flow 
transition  from  tension  in  the  falling  stage,  to  compression  in 
the  deceleration  stage  when  the  jet  strikes  the  plate. 

Commenting  on  Suleiman  and  Munson’s  paper,  Bejan  [2] 
noted  that  the  buckling  phenomenon  is  not  only  a  property  of 
highly  viscous  fluids,  but  also  a  property  of  inviscid  columns 
(streams),  and,  as  such,  explains  the  “meander” 
phenomenon.  The  process  of  inviscid  stream  buckling  was 
treated  analytically  by  Bejan  [12]  who  showed  that  all  inviscid 
fluid  layers  buckle  so  that  the  wavelength  is  proportional  to 
the  layer  thickness  only.  The  proportionality  between 
buckling  wavelength  and  jet  diameter  was  verified  recently  in 
two  separate  experiments,  one  involving  the  meandering  of  air 
streams  driven  by  falling  paper  ribbons  (Bejan  [13])  and  the 
other  focusing  on  the  buckled  shapes  of  fast  capillary  jets 
(Stockman  and  Bejan  [14]). 

At  this  stage  in  our  understanding  of  fluid  buckling,  an 
important  discrepancy  exists  between  the  observed  buckling 
behavior  of  viscous  layers  and  the  behavior  predicted  by  the 
viscida  theory  [7],  The  theoretical  wavelength  of  the  buckled 
shape  is  indeterminate  and,  presumably,  dictated  randomly 
by  the  original  deformation  (disturbance)  of  the  straight 
layer.  Experimental  observations,  on  the  other  hand,  seem  to 
suggest  that  the  buckling  wavelength  is  not  random  (1,  5,  11]. 
Thus,  the  objective  of  the  experimental  work  described  in  this 
paper  is  to  establish  whether  the  buckling  wavelength  of 
viscous  layers  in  indeterminate,  as  in  the  viscida  theory,  or,  in 
fact,  a  "characteristic”  length.  The  experiment  does  not 
correspond  fully  to  the  geometry  of  Taylor  [5]  or  Buckmaster 
et  al.  [7,  8],  however,  it  does  shed  light  on  the  uniqueness  of 
the  observed  buckling  wavelength. 

2  Experiment 

The  experiment  was  designed  to  measure  the  buckling 
wavelength  of  a  viscous  fluid  layer  and  to  determine  quan¬ 
titatively  what  parameters  affect  the  wavelength.  The  ex¬ 
periment  is  shown  schematically  in  Fig.  1.  A  layer  of  viscous 
fluid,  floating  on  a  relatively  inviscid  fluid,  is  compressed  by 
the  partition.  The  partition  extends  into  the  inviscid  fluid  a 
distance  on  the  order  of  the  viscous  layer’s  thickness  and  has  a 
rate  of  compression  u.  The  geometrical  shape  of  the  layer  is 
given  by  the  slenderness  ratio  L/d  where  L  and  d  are  the 
length  and  thickness  of  the  unbuckled  layer.  The  slenderness 
ratio  was  varied  from  about  20  to  70  by  changing  the 
thickness  of  the  layer  and  leaving  the  length  constant . 

The  experiment  was  carried  out  in  an  aluminum  box  with 
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Fig.  2  Construction  details  of  the  experimental  apparatus 


inside  dimensions  measuring  20.3  cm  in  length,  7.6  cm  in 
width  and  10.2  cm  in  depth.  The  sides  of  the  box  serve  as  a 
track  for  a  sliding  mechanism  carrying  a  partition  (gate) 
which  can  be  lowered  into  the  box  (see  Fig.  2).  It  is  by  means 
of  this  partition  that  the  layer  is  compressed.  The  sliding 
mechanism  was  designed  to  operate  without  vibrations  and  to 
hold  the  partition  perpendicular  to  the  viscous  fluid  layer  and 
to  the  walls  on  each  side.  It  also  holds  the  partition  in  contact 
with  the  side  walls  preventing  the  fluid  layer  from  escaping 
behind  the  partition  as  the  viscous  layer  is  compressed. 

The  sliding  mechanism  was  moved  by  hand;  in  order  to 
bring  uniformity  to  the  results,  two  stops  were  positioned  to 
set  the  initial  length  L  and  to  control  the  maximum  excursion 
of  the  partition.  As  the  sliding  mechanism  moved  between  the 
stops,  it  passed  over  a  microswitch  sending  an  electrical  pulse 
to  a  timer  and,  at  the  same  time,  triggering  the  camera.  The 
duration  of  the  pulse  was  used  to  make  a  time  of  flight 
measurement,  and  with  the  assumption  that  the  partition 
undergoes  a  constant  acceleration,  it  was  possible  to  deter¬ 
mine  the  rate  of  compression  at  any  instant.  The  accuracy  of 
the  velocity  measurement  is  determined  by  the  constant  ac¬ 
celeration  assumption.  Since  the  compression  stroke  is  very 
short  (13  to  32  milliseconds)  and  the  force  applied  by  hand 
during  that  interval  is  approximately  constant,  the  ac¬ 
celeration  is  also  approximately  constant. 

Experiments  were  run  to  determine  the  effects  of  the  ex¬ 
cursion  length  and  rate  of  compression  on  the  buckling 
wavelength.  In  varying  the  excursion  length,  it  was  found  that 
there  exists  an  optimum  range  of  lengths  for  the  most  accurate 
portrayal  of  the  buckling  phenomenon.  This  range  was  from  1 
cm  to  2  cm.  If  the  excursion  length  was  less,  there  was  no 
visible  buckling,  and  if  it  was  greater,  the  buckling  became  so 
pronounced  that  the  “waves”  would  collide.  The  effect  of  the 
rate  of  compression  is  discussed  in  greater  detail  in  the  next 
section. 

It  is  important  to  note  that  although  the  buckled  layer  has 
the  appearance  of  a  propagating  gravity  wave,  it  is  not.  This  is 
demonstrated  by  the  observation  that  after  ihe  partition  has 
stopped,  the  buckled  shape  remains  stationary.  In  time,  the 
amplitude  decreases  and  eventually  the  layer  reaches 
equilibrium  again.  In  many  cases  the  buckled  shape  and  its 
wavelength  remained  visible  and  could  be  photographed  for 
up  to  30  minutes  after  the  compression  stroke.  As  shown  in 
section  4,  the  observation  that  the  amplitude  of  the  buckled 
shape  decreases  away  from  the  moving  partition  (Figs.  3,  4i 
can  be  explained  by  the  fact  that  the  layer  is  not  uniformly 
compressed.  The  nonuniform  compression  rate  can  be  due  to 
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Fig.  3  Photograph  showing  the  buckled  layer  as  seen  Irom  above  (the 
partition  moves  from  left  to  right) 

the  no-slip  condition  along  the  side  wall,  or  to  the  inertia  of 
the  fluid  layer  itself. 

The  buckling  wavelength  was  measured  photographically 
using  a  camera  triggered  by  the  trailing  edge  of  the  pulse 
created  by  the  microswitch.  Special  lighting  conditions  were 
needed  in  order  to  make  the  buckling  wavelength  easy  to 
measure  on  the  photographs  (Blake  [15]).  This  system  enabled 
one  side  of  a  “wave  crest”  to  be  illuminated  while  leaving  the 
other  side  in  darkness.  The  photograph  then  contained  a  series 
of  light  and  dark  bands  where  one  light  and  one  dark  band 
together  account  for  one  wavelength.  The  sample  photograph 
shown  in  Fig.  3  has  three  bands,  implying  that  the  buckled 
region  extended  one  and  one  half  wavelengths  into  the  viscous 
fluid  layer.  The  other  photographs  contained  between  two 
and  five  bands  (see  also  Fig.  4). 

In  all  the  experiments  involving  data  acquisition,  Dow 
Corning  200  silicone  oil  was  used  as  the  viscous  fluid  and 
distilled  water  as  the  inviscid  substance.  The  silicone  oil  has 
several  desirable  attributes  such  as  its  small  variation  in 
viscosity  with  temperature,  its  low  surface  tension,  its  im- 
miscibility  with  water,  and  its  availability  over  a  wide  range  of 
high  viscosities.  In  addition,  any  intermediate  viscosity  can  be 
achieved  simply  by  mixing  a  higher  and  lower  viscosity  in  the 
proper  proportions.  Since,  as  shown  in  the  next  section,  the 
buckling  wavelength  is  not  strongly  dependent  on  viscosity,  it 
was  not  necessary  to  know  the  precise  viscosity,  therefore,  the 
oil  viscosity  reported  in  these  experiments  is  the  viscosity 
calculated  based  on  Dow  Corning  information.  Furthermore, 
the  experiments  were  run  at  room  temperature,  which  is 
within  a  few  degrees  of  the  temperature  (75 *F)  at  which  the 
viscosity  was  measured.  The  viscosity  of  the  silicone  oil  used 
in  the  experiments  ranged  from  103  to  10s  centist'okes  (cSt) 
where  one  centistoke  is  approximately  the  viscosity  of  water 
and  the  units  of  stokes  are  centimeters  squared  per  second. 
Lastly,  the  silicone  oil  was  demonstrated  to  be  a  Newtonian 
fluid  by  Suleiman  and  Munson  [9] . 

To  set  up  the  experiment,  a  flat,  uniform  layer  of  the 
viscous  fluid,  of  a  prescribed  thickness,  must  be  placed  on  the 
surface  of  the  water.  The  method  used  was  to  put  the  silicone 
oil  in  a  standard  100  ml  burret  with  the  following 
modifications:  the  tapered  end  of  the  burret  was  ground  off, 
and  the  hole  in  the  valve  was  drilled  out  to  be  the  same  size  as 
the  tube.  The  burret  was  held  about  a  centimeter  above  the 
surface  of  the  water  and  the  oil  was  allowed  to  drain  out  until 
the  desired  thickness  was  obtained.  Depending  on  the 
viscosity,  it  took  from  one  to  five  hours  to  drain  out  the 
required  amount.  This  extremely  slow  pouring  rate  allowed 
the  oil  to  slowly  creep  across  the  surface  of  the  water 
providing  a  flat,  smooth  layer.  Knowing  the  volume  of  oil  and 
the  area  of  the  layer,  the  thickness  d  could  be  calculated  to  an 


Fig.  4 (a)  da .394 ±.003  cm,  L  =  16. 35 ±.03  cm,  u  =  42.3±.2  cm/s, 
p  a  10*  ±400  cSt 
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Fig,  4<6)  d  =  .264 ±.003  cm,  1  =  16.35 ±  03  cm.  u  =  32.6±.2  cm/s, 
»  =  105  ±4000  cSt 

Fig.  4  Two  separate  experimental  runs: 

accuracy  of  0.003  cm  with  most  of  the  error  resulting  from  the 
meniscus  on  the  edges  of  the  oil  layer. 

Since  the  experiment  sometimes  took  days  to  run,  de¬ 
aerated  (distilled)  water  was  needed  for  the  lower  fluid. 
Otherwise  bubbles  of  air  would  form  under  the  layer  of  oil 
creating  little  bumps  in  the  surface. 

4  Results 

Despite  the  fact  that  in  the  present  experiment  the  motion 
of  the  position  is  not  mechanized,  it  is  important  to  document 
the  effect  of  changes  in  the  rate  of  compression  u.  It  is  also 
important  to  report  the  minimum  rate  of  compression 
necessary  for  buckling  the  layer.  Referring  to  buckling  in  a 
purely  two-dimensional  layer  geometry,  Taylor  [5 ]  stated  that 
when 

-4  ptrf-r>o  (i) 

the  viscous  sheet  is  unstable  to  disturbances  of  any 
wavelength,  where  n  and  e  are  the  viscosity  and  the  rate  of 
strain.  The  total  surface  tension  T  is  the  sum  of  the  oil-air 
interface  TUJ  and  the  water-oil  interface  TH„.  For  two  im- 
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Fig.  5  The  weak  dependence  between  buckllr.g  wavelength  (A)  and 
rate  of  compression  (u).  (Uncertainty  in  A/d=  *0.08,  in  4(p»u/7V 
(L/d>  =  ±0.06  at  20:1  odds) 


miscible  liquids,  the  value  of  T„0  is  approximately  Twa  -  T0j, 
where  Twa  is  the  surface  tension  of  the  water-air  interface 
(Suleiman  and  Munson  19),  Davies  and  Rideal  [16]). 
Therefore,  T  =  Tna  +  (T^-T*)  =  r„„,  that  is,  the  total 
surface  tension  T  is  just  the  surface  tension  of  the  water-air 
interface.  Upon  substituting  - u/L  for  e  in  Taylor’s  equation 
and  rearranging,  it  is  found  that  the  layer  buckles  when 

4(pvu/T)/(L/d)>  l  (2) 

It  should  be  noted  that  Taylor’s  conclusions  do  not  apply 
exactly  to  the  present  experiment.  As  shown  i..  Figs.  3  and  4, 
the  no-slip  condition  along  the  side  walls  gives  rise  to  a  three- 
dimensional  flow  such  that  only  the  middle  portion  of  the 
buckled  layer  shows  parallel  waves  and  can  be  regarded 
approximately  as  two-dimensional.  Therefore,  the  sub¬ 
stitution  «  =  -u/L  used  to  derive  equation  (2)  can  only  be 
approximately  valid.  That  —u/L  is  not  the  rate  ot  strain 
everywhere  in  the  layer  is  demonstrated  by  the  fact  that 
buckling  does  not  occur  throughout  the  layer  but  only  near 
the  compressed  end.  Figures  3  and  4  show  that  sufficiently 
close  to  the  partition,  the  bands  (waves)  are  relatively  constant 
size,  hence,  it  is  reasonable  to  assume  that  in  that  region  the 
rate  of  strain  is  relatively  constant.  In  cases  where  the 
photographed  bands  are  unequal  in  size,  the  reported 
wavelength  was  calculated  by  averaging  the  first  two  bands. 

For  this  experiment  the  viscous  fluid  layer  had  a  slenderness 
ratio  (L/d)  of  25,  a  viscosity  (v)  of  3  x  104  cSt,  a  density  (p)  of 
0.975  g/crn',  and  a  surface  tension  ( T)  of  72.8  dyne/ cm.  The 
results  of  the  experiment  are  given  in  terms  of  the  buckling 
wavelength  to  thickness  ratio  X/e/  versus  4(pvu/T)/(L/d)  and 
graphed  as  shown  in  Fig.  5.  The  data  fall  into  two  categories: 
points  where  the  fluid  buckled  and  points  where  the  fluid  did 
not.  The  cases  in  which  buckling  did  not  occur  have  no  \/d 
value  and  are  indicated  with  short  vertical  lines  at  the  bottom 
of  the  graph.  The  large  vertical  line  indicates  the  approximate 
minimum  value  of  4(pvu/T)'{L/d)  that  will  buckle  the  layer. 
This  value  is  approximately  7,  i.e..  greater  than  the  value  1 
appearing  in  equation  (2'  derived  trom  Taylor’s  buckling 
criterion  (1).  This  discrepancy  is  not  surprising  in  view  of  the 
three-dimensional  effects  and  nonuniform  t  that  distinguish 
the  present  experiment  from  Taylor’s  two-dimensional 
geometry. 

An  important  result  of  the  present  experiment  is  the  in¬ 
dependence  of  \id  on  the  rate  of  compression,  as  witnessed 
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Fig.  6  The  relationship  between  buckling  wavelength  (A)  and  both 
length  to  thickness  ratio  (L,d)  and  viscosity  (>).  (Uncertainty  in 
\/d=  ±0.06,  in  L/d  =  ±0.01,  In  »  =  ± 0.04  at  20:1  odds.) 

by  the  nearly  horizontal  distribution  of  the  data  on  Fig.  5.  The 
data  was  limited  on  the  right  by  the  inability  to  compress  the 
layer  at  a  faster  rate.  This  conclusion  is  further  supported  by 
the  results  given  in  Fig.  6  where  each  datum  is  an  average 
obtained  from  three  different  rates  of  compression.  In  each 
set  of  three,  the  group  4 (pvu/T)/(L/d)  varies  with  an  average 
55  percent  increase  from  the  smallest  to  the  largest  value. 
However,  the  corresponding  values  of  \/d  only  have  an 
average  5  percent  increase,  as  indicated  approximately  by  the 
size  of  the  circles.  Thus  for  a  variety  of  L/d  ratios  and 
viscosities  it  can  be  seen  that  the  sizable  change  in 
4  (pvu/  T)  /(L/d)  has  only  a  marginal,  if  any,  effect  on  \/d. 

To  determine  the  relationship  between  the  buckling 
wavelength  and  the  slenderness  ratio,  a  series  of  experiments 
were  conducted  varying  the  thickness  of  the  layer  and  using 
different  viscosity  oils.  The  data  were  obtained  in  the 
following  manner.  The  partition  was  placed  at  the  maximum 
distance  from  the  far  end  so  that  L  =  16.4  cm.  The  box  was 
filled  with  water  and  then  a  thin  layer  of  oil  was  poured  onto 
the  water  surface  from  the  burret.  To  obtain  different 
slenderness  ratios  more  oil  was  added  from  the  burret,  and 
allowed  to  reach  equilibrium.  For  each  slenderness  ratio  the 
layer  was  usually  compressed  three  times.  Between  each 
compression  the  partition  was  moved  back  to  its  initial 
position  and  the  oil  layer  was  given  enough  time  to  reach 
equilibrium  again.  For  each  change  in  viscosity,  the  burret 
was  drained  of  the  previous  oil  and  the  new  viscosity  oil  was 
poured  in. 

The  results  of  these  experiments  are  graphed  as  \/d  versus 
L/d  in  Fig.  6.  The  plot  contains  three  sets  of  points  indicated 
by  an  open  circle  and  right  and  left  hall -filled  circles  tor 
viscosities  of  10!  cSt,  104cSi,  and  105cSt,  respectively.  It  is 
clear  from  Fig.  4  that  \/d  is  dependent  on  L/d.  For  L/d  less 
than  40  the  slopes  of  the  curves  are  almost  unity,  indicating 
that  X  is  proportional  to  L.  Due  to  the  construction  of  the 
experimental  apparatus,  the  value  of  L/d  could  not  be  in¬ 
creased  beyond  65. 

The  effect  of  viscosity  on  the  buckling  wavelength  can  also 
be  interpreted  from  the  data  in  Fig.  6.  The  value  ot  X  </ 
decreases  by  about  one-third  when  changing  the  viscosity 
from  10'  cSt  to  10'  cSt  and  keeping  L/d  constant.  Thi^ 
change  is  fairly  consistent  throughout  the  range  of  values  oi 
L/d.  Within  this  range,  then,  the  effect  of  viscosity  on  \/J  is 
seen  to  be  very  slight. 

The  error  in  the  results  presented  here  is  governed  primarily 
by  the  measured  value  ot  X.  The  variables  d,  L,  and  u  have  all 
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been  calculated  to  within  1  percent.  The  other  variables  in¬ 
fluencing  the  results,  although  unknown  precisely,  remain 
constant  throughout  each  experiment  and  therefore  do  not 
effect  the  overall  trends  shown  in  the  graphs. 

5  Conclusions 

The  object  of  this  experimental  report  has  been  to  present  a 
series  of  observations  concerning  the  buckling  wavelength  of 
a  highly  viscous  fluid  layer  undergoing  end-compression.  The 
motivation  for  designing  and  running  these  experiments 
stemmed  from  the  unresolved  issue  of  whether  the  buckling 
wavelength  is  arbitrary  (as  in  the  viscida  theory,  Buckmaster 
et  al.  [7]),  or  a  characteristic  length  of  the  layer.  Based  on  the 
measurement  produced  by  the  present  study,  the  following 
conclusions  may  be  drawn: 

1  In  a  fluid  layer  of  fixed  geometry  and  viscosity,  the 
buckling  wavelength  is  practically  independent  of  the  rate  of 
compression,  that  is,  independent  of  the  layer  velocity  relative 
to  the  fluid  ambient  (Fig.  5). 

2  The  viscosity  of  the  buckled  layer  has  only  a  minor 
impact  on  the  buckling  wavelength  (Fig.  6). 

3  In  the  range  L/d<  100  the  buckling  wavelength  X  scales 
with  the  length  of  the  compressed  layer,  L  (Fig.  6). 

In  view  of  these  findings,  the  wavelength  of  a  buckled 
viscous  layer  emerges  as  a  property  of  the  layer  (a  charac¬ 
teristic  length),  as  opposed  to  the  undetermined  length  dic¬ 
tated  by  random  initial  disturbances  assumed  in  the  viscida 
theory. 
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